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THE HYDRAULIC JUMP AT AN ABRUPT DROP! 


Walter L. Moore,* A.M., ASCE and Carl W. Morgan,** A.M., ASCE 
(Proc. Paper 1449) 


SYNOPSIS 


The hydraulic jump may form at various locations relative to a low abrupt 
Op in a rectangular channel. The role of the drop in determining the form 

the jump and in stabilizing its position is clarified by analysis and experi- 
sent. An example illustrates the application of the results to the analysis of 
stilling basin. 


INTRODUCTION 


In many types of hydraulic structures a considerable portion of the kinetic 
ergy of the flowing water must be dissipated to prevent erosion of the chan- 
| downstream from the structure. Various methods of energy dissipation 
ve been uSed to achieve tranquil flow conditions as the flow enters the down- 
‘eam channel. The primary purpose of all such methods is to convert as 
ich as possible of the kinetic energy of the flow into turbulent energy and 
imately into heat. The required energy dissipation is frequently most ef- 
tively accomplished by means of a hydraulic jump which may assume 

yeral different forms depending on the geometry of the channel boundaries. 
an effort to reduce the length of stilling basins and to insure that the jump 
1 not be swept from the stilling basins and to insure that the jump will not 
Swept from the stilling basin under varying flow conditions, an abrupt drop 
rise in the channel bottom may be introduced. 

For economic reasons, there has been a tendency in practical design to 

ke the length of the stilling pool considerably less than the length of the 

ap. This paper presents information concerning the characteristics of the 


é: Discussion open until May 1, 1958. Paper 1449 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers ,Vol. 
83, No. HY 6, December, 1957. 

Paper presented at a meeting of the American Society of Civil Engineers in 
Knoxville, Tenn., June, 1956. 

Prof., Dept. of Civ. Eng., Univ. of Texas, Austin, Tex. 

Asst. Prof., Dept. of Civ. Eng., Univ. of Texas, Austin, Tex. 
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jump at an abrupt drop in the channel bottom with particular emphasis on th 
characteristics when the drop is located within the length of the jump. 


Previous Study 


Basic information concerning the characteristics of the hydraulic jump 0 
a level surface is published in many books on hydraulics and fluid mechanic 
Mr. Edward A. Elevatorski(1) has compiled a bibliography of published infor 
mation concerning the hydraulic jump under various channel and flow condi- 
tions. This publication lists some 504 references arranged by years startin 
in 1819 and extending to 1954. 

Information concerning flow problems closely allied with those discussed 
in this article have been published. Moore(2) made an analytical and experi: 
mental study of a free overfall and determined the energy, velocity, and sur: 
face profiles for various conditions of flow. Forster and Skrinde(3) studied 
both theoretically and experimentally the stabilizing effect of an abrupt rise 
at the end of a jump and developed a composite graphical summary of result 
by which the downstream depth after the rise could be computed. Valuable 
contributions to the flow analysis were also made by the several persons dis 
cussing this paper. Rouse, Bhoota, and Hsu(4) investigated the stabilizing 
characteristics of an abrupt drop in a channel bottom. They found that two 
different types of flow may occur at a given drop, the transition from one to 
the other being characterized by an undular wave. The type of flow which w: 
form is dependent upon whether the downstream depth is above or below tha 
which produces the undular standing wave. The conditions which produce th 
undular stage (critical transition zone) cannot be foretold by equations utiliz 
ing only the momentum and continuity relationships but must be obtained by 
experimental measurements. Rouse, Bhoota, and Hsu presented such data 
which verified their analysis and indicated a systematic transition between 
the two regimes of flow depending on the Froude number. To simplify the 
analysis for the two types of jumps, the pressure on the face of the drop wa: 
assumed to be hydrostatically distributed and assumed to be dependent on tt 
downstream depth if the jump is upstream from the drop and dependent on tl 
upstream depth if the jump is downstream from the drop. Experimental dat 
was not presented for intermediate positions of the jump which were not suk 
ject to a simple analysis based on the two assumed hydrostatic pressure dis 
tributions on the face of the drop. q 


. 


Ml Analysis : 
The different forms of a hydraulic jump which will occur at an abrupt dr 
in a rectangular channel are shown in the definition sketch, Fig. 1. For con 
venience the different forms are identified as A, W, and B as indicated on qd 
Y 
sketch. The relative height of a stationary jump = will be dependent upon 
1 
AZ 


; 

the dimensionless parameters 2 . and the Froude tee of tk 
a6 vo lov 4 

1 1 y sY, 

entering flow. As the jump changes from one form to the other, the te q 
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y,=Depth upstream from drop 
Yo= Depth downstream from drop 
4z= Height of drop 

Np= Piezometric head on drop 


FIG. |. FORMS OF HYDRAULIC JUMP AT AN ABRUPT DROP. 
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the hydrostatic head, hp, will change causing the pressure on the face of th 
drop to change from a maximum value of yhp = y Y2 when the jump forms 
stream from the drop to a minimum value of yhp< y (AZo + Y1) when the 
jump forms downstream from the drop. Photographs illustrating various 
forms of the jump are shown as Fig. 2. 

If the pressure distribution on the face of the drop is assumed to be hyd 
static under a head of hp, if the shear stress along the solid boundaries be 
tween section 1 and 2 is neglected, and if the momentum coefficients B; an 
Bo are assumed as unity, the momentum equation for sections (1) and (2) is 


2 


2 
Yq ae! 12 
soe ees = 
+ v(oy ; AZ, * av, V,) 


Zz 


in which ¥ is the specific weight of the fluid and q is the discharge per unit 
width of the channel. If the continuity equation, 


q=V,Y, = V.Y ( 


1h | 22 


and Eq. 1 are combined and reorganized, the following general equation rel 
ing the several variables is obtained. 


a 


2 
s AZ, AZo fhe AZ, 
ye me ny sews Y cai ice 
2 1 1 1 1 Y) ( 


Plots of Eq. 3 for different values of the variables are shown as Figs. 3, 
4, and 5. 


It can be shown that for the condition hp = Y2 (Fig. 1A) Eq. 3 reduces to 


ce Kee he 
weet 
2 Y Y - 1 


1 i 


04 
om he (eee 
Y2 ] 
A plot of this equation gives the upper limiting curve (plotted as a dashe 
line) for Figs. 3, 4, and 5. 


| a ae 
cut 4 : One. 
For the condition hy Tee ae 25) that is yy, ae vy, (Fig. 1B) the 
general equation reduces to 


E MOORE - MORGAN 1449-5 


| MIN. JUMP B 


FIG. 2. PHOTOGRAPHS SHOWING FORMS OF THE 
HYDRAULIC JUMP AT AN ABRUPT DROP 
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FIG.3 RELATIVE DOWNSTREAM DEPTH AS A FUNCTION OF 
FROUDE NO. AND PIEZOMETRIC HEAD AT THE  OROP. 
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FIG.4 RELATIVE DOWNSTREAM DEPTH AS A FUNCTION OF 
FROUDE NO. AND PIEZOMETRIC HEAD AT THE DROP. 
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Solid pass are plots 
of Equation (3) 
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FROUDE NO. AND PIEZOMETRIGC HEAD AT THE DROP 
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Equation 4 and 5 are equivalent to the equations given by Hsu, et al, (5) for 
the two values of hp which were assumed to determine the characteristics 
analytically. 


It is noted that a reduction of y_ causes a shift to the right in the curves of 
1 
Figs. 3, 4.and 5. That is, for a given condition of the entering flow, the jump 
can be held in place with a reduced downstream depth Y9. The converse is, 


of course, equally true. 


Experimental Apparatus 


Experimental observations were made on two drops which were con- 
structed of 3/4 inch waterproof plywood and each in turn mounted in a one- 
foot wide glass walled flume. One drop was 0.400 ft. deep and the other was 

0.200 ft. deep. Each was 3 ft. long with an elliptical section at the upstream 
end and a vertical section at the downstream end. A brass sluice gate adjust- 
ee by an attached screw device could be positioned to 0.001 ft. The sluice 

te had a carefully machined 1/2 inch rounded edge which gave only a slight 

Be acticn in the jet. 

_ Water was supplied from a constant head tank through an 8 in. pipe to the 
-forebay of the flume. The discharge was regulated with an 8 in. gate valve, 

and the flow rate was measured by a 5 1/4 x 8 in. Venturi-Orifice Meter with 

a 4 3/8 in. sharp edged orifice plate inserted in the throat section. The dif- 

ferential head on the meter could be read to t .001 ft. which made it possible 
0 set the rate of flow with an error less than 1%. 

Surface elevations were determined with an electric point gage mounted 

“upon an instrument carriage that could slide along rails above the channel. 
‘The vertical position of the gage was read with a direct vernier reading to 
0.001 ft. The depth of water in the downstream channel was controlled by a 
Wicket type gate, which caused a minimum disturbance to the flow pattern up- 
‘Stream. It was located 17.5 feet downstream from the last section at which 
‘Measurements were taken. 
Velocities were measured with a standard Prandtl pitot-static tube that 
had a coefficient of 1.00. The pitot tube was fitted to the same instrument 
Carriage which held the point gage and could be positioned to 0.001 ft. Velo- 
‘Cities were measured near the bottom of the channel at several positions 
downstream from the drop. Distances were measured on a section of engi- 
Neers steel tape with an error of no more than t 0.002 ft. 

Values of hp were measured by means of piezometers attached to 

piezometric holes drilled in a 1/2 in. x 1 in. brass strip set into the face of 
the drop so as to be flush with the vertical surface. The openings were set at 

0.073 ft., 0.200 ft. and 0.327 ft. from the bottom of the 0.400 ft. drop and at 

0.100 ft. from the bottom of the 0.200 ft. drop. 
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Z 4 eX E Type of Jump 
1 Y, 1 
0.400 0.100 4 2,4, 6,8 B,W,A 
0.400 0. 133 3 2,4 B,W,A 
0.200 0.067 3 6,8 B,W,A 
0.200 0.100 2 2,4,6 B,W,A 


Values of the depth Yo, the piezometric head at each piezometer, and 


The following flow conditions were chosen for observations. 
velocity along the bottom were read for each of the flow conditions listed. : 


Results 


It was found that the pressure on the face of the drop was hydrostatically 
distributed under a piezometric head, hp, which varied with the location of 


the pump relative to the drop. The measured values of Y are shown on 
1 

Yo 

Figs. 3, 4, and 5 at the values of YY and Fy at which they were observed. 
1 

The decimal point indicates the plotted point. These observed values confirn 
the analytical calculations within the limits of accuracy of the observation an 
clearly indicate the variation of hp with the change in the flow conditions. 


The limiting values of < within which each type of jump would form at 
1 
given entering flow conditions are shown in Figs. 6, 7, and 8. The conditions 
under which Jump A will form are shown as a very narrow band, but actually 
ny. 
Jump A may form at values of = greater than the upper bound being limited 
al 
only by the free distance upstream and the slope and resistance of the chan- 
nels. The undular wave will form at values of the variables which fall within 
the middle band. The height of this wave may exceed the value of Yo by as 
much as 50 per cent and might cause considerable damage unless adequate 
training walls were provided. Jump B will form when values of the variable: 
are those within the lower band (dotted). The narrow gap between the inter- 
mediate and lower bands in Figs. 6, 7, and 8 represents a transition zone 
where there is uncertainty as to whether the wave or Jump B will form. For 
Y 
values of < and F, which would lie below the lower band, the jump would be 
1 S , 
swept downstream from the drop. Perhaps the most important characteristi 
ay. % 
presented is that the jumps can form over a wide range of = for the entire 
1 
range of entering Froude numbers covered by the experiments. This might 
be anticipated from the previous analysis because the hydrostatic head, hp, 
should vary continuously between its limiting values. Visual observation of 
the gradual change in curvature of the main stream of the flow as it crossed 
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FIG.6 FORMS OF THE HYDRAULIC JUMP AS A FUNCTION 


OF FROUDE NO. AND RELATIVE DOWNSTREAM DEPTH. 
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FIG.7 


OF FROUDE NO. AND RELATIVE DOWNSTREAM DEPTH. 
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FIG. 8 FORMS OF THE HYDRAULIC JUMP AS A FUNCTION 


OF FROUDE NO. AND RELATIVE DOWNSTREAM DEPTH. 
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the drop also indicated a continuous variation of hp. As indicated by Eq. 3, 


; 2 : 
continuous variation of hp implies a continuous variation of Yy between its 
1 
extreme values. 
To demonstrate the effectiveness of the different types of jumps in reduc 
ing the velocity along the channel bottom, the velocity near the channel bott 
is plotted against distance from the drop in Figs. 9, 10, and 11. Of particul 


interest are the high velocities produced by Jump B which may give ratios i 


of 5 or 6 at = of approximately 2. The wave and Jump A give values of 
2 

Vy 
Vv. from 2 to 3 in the upstream direction near the drop. The return flow re 

2 
sults from the reverse roller produced by both the wave and the Jump A. 
Jump A gives the smallest values of velocities along the bottom in the down: 
stream direction. This condition would be expected since the greater pres- 
sure on the face of the drop causes the entering flow to be directed slightly 
upward and the violent eddy motion to occur in the upper levels of the down- 
stream channel. 


Application of Results 


The investigation has clarified the action of an abrupt drop in stabilizing 
the position of a hydraulic jump. The results obtained are directly applicah 
to the design of a stilling basin incorporating an abrupt drop. 

If the downstream depth, Yo, is required for given entering flow conditio 
it may be read from Fig. 6, 7, or 8 using the chart with the applicable relat 


AZ 
drop height. For example, given Y, = 2 ft., Froude No. = 4 and r = 3. 
1 
Then, from Fig. 7, the following information could be obtained: 
rs 
Type of Jump Range of —— Range of Y 
Yy 2 
Jump B 5.1 - 6.8 ; 10.72) =" 13 6.fts 
Wave 7.0 - 7.6 14.0 - 15.2 ft. 
Jump A 7.6 -7.8 L5ecr— (57 the 


The results may be applied to analyze the performance of a stilling basi 
for a range of discharges and tailwater elevations. For example, assume a 
spillway with a crest 100 ft. long at elevation 460 and a coefficient of 4.0. | 
the apron be at elevation 400 with a drop of 8 ft. to elevation 392. Fora 
series of values for the head on the spillway, discharge may be determined 
well as the depth and Froude number at the base of the spillway; i.e., Yq ar 
Fj. (For convenience, the effects of resistance and air entrainment were r 
lected in the example.) By use of Figs. 6, 7, and 8, and by interpolating for 
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FIG. 9. VELOCITY ALONG BOTTOM Vs. DISTANCE FROM DROP 
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FIG. 12 TAILWATER ELEVATIONS AND JUMP 
CHACTERISTICS FOR DOWNSTREAM CHANNEL I 
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AZ 
ntermediate values of on , the range of values of Yo for each type of jump 
an be determined. These results are shown in Fig. 12 in which the depths 
'9 have been converted into elevations. Also shown in Fig. 12 are curves of 
he tailwater elevations at sections 2 and 3 for a rectangular downstream 
hannel 100 ft. wide with its bottom at elevation 398.5, a slope of 0.0016, and 
‘Manning n of 0.03. This is referred to as Channel I. In making the compu- 
ations it was assumed that the depressed apron (elevation 392) was 125 ft. 
ong which would allow tranquil flow to develop at section 2. If the depressed 
pron were shorter so that high velocity flow would impinge on the step at the 
nd of the apron, dynamic pressures would develop on the face of the step 
dding an additional upstream force component tending to stabilize the jump. 
ince it is not yet possible to evaluate the effect of this force, it was elimi- 
ated from consideration by choosing a long apron. 

From Fig. 12, it is apparent that the type A jump will occur at discharges 
bove about 20,000 cfs with the form changing to the wave at lower rates of 
ischarge. There is no danger of the jump washing downstream for these 
Onditions and from Figs. 9, 10, and 11, it is apparent that the velocity near 
he bottom would be negligible while the type A jump forms; i.e., Q > 20,000 
fs. At lower discharges when the wave forms the velocity near the bottom 
t the end of the depressed apron would not exceed 3V9 or about 3 x 7 = 21 
./sec. 

Figure 13 shows what would happen if the same spillway discharged into a 
hannel 100 ft. wide with its bottom at the same elevation as the depressed 
pron (392) and with a slope of 0.0010 and a Manning n of 0.035. This is re- 
srred to as Channel Il. From Fig. 13, it is apparent that type A jump would 
rm at the maximum discharge with the form changing successively through 
le wave to type B as the discharge decreases. As the form of the jump 
hanges, it would move progressively down over the drop but, even at the 
west tailwater, it would remain close to the location of the drop. 

These examples illustrate the type of analysis made possible by the re- 
ults of these experiments. By making it possible to determine the jump 
haracteristics for a wide range of Froude number and tailwater conditions, 
ie effect of both upstream and downstream geometry on the jump character- 
stics can be determined. 


. CONCLUSION 


r S A 

1. An abrupt drop in the bottom of a rectangular channel is effective in 
fabilizing the hydraulic jump over a broad and continuous range of values of 
relative downstream depth. 


2. As the relative downstream depth is reduced from its upper limit, the 

p shifts progressively downstream resulting in a corresponding continuous 
rease in the pressure on the face of the drop as the jump passes through 

e different forms. 


3. A momentum analysis which includes a term for the pressure force on 
vertical face of the drop gives results which are consistent with the ex- 
imental observations. 
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4. The maximum velocity attained near the channel bottom, which may be 
considered as an index of the scouring power of the flow leaving the jump, 
varies over a wide range and depends mainly on the form of the jump. For 
the type B jump (Fig. 1) this velocity may attain a value several times greater 
than the downstream velocity, but for the type A jump, the maximum velocity 
near the bottom exceeds only slightly the average velocity in the downstream 
channel. 
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5 TRANSITION FROM LAMINAR TO TURBULENT FLOW IN A PIPE 
J 
. M. R. Carstens! 
4 (Proc. Paper 1450) 
ABSTRACT 


a 


_ Transition from laminar to turbulent flow was experimentally observed 
du ing flow establishment in a smooth straight pipe. Turbulence originated 
spots. The downstream face of the turbulent spots was transported with 
the centerline velocity of the preceding laminar flow. Turbulent spot forma- 
* was in qualitative agreement with the Tollmein-Schlichting theory of 
§mall disturbances. 


INTRODUCTION 


A study of flow establishment in a pipe was performed in the Hydraulics 
aboratory, School of Civil Engineering, Georgia Institute of Technology. 
he experiments were performed with a smooth, straight, circular pipe which 
horizontally aligned. The upstream extremity of the pipe was a well- 
ded inlet which was located in a large reservoir. The downstream ex- 
emity of the pipe was unobstructed so that a liquid jet was formed in the at- 
osphere downstream from the pipe. The flow was established by rapid 
ffectively instantaneous) release of a disk valve placed against the square- 
ged downstream end of the pipe. Velocity-time data were obtained from a 
ion picture record of the jet. Pressure-time data were obtained at 
lected stations along the pipe. These data were recorded as an oscillo- 
raph record of the output signals of pressure transducers. 
The experiments were performed with a systematic variation of the inde- 
endent dimensionless variables in the equations: 


@: Discussion open until May 1, 1958. Paper 1450 is part of the copyrighted 

Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 

83, No. HY 6, December, 1957. 

Assoc. Prof., School of Civ. Eng., Georgia Inst. of Technology, Atlanta, 
Ga. 
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and 


The letter symbols are defined as follows: 


D - inside diameter of pipe; 

g - acceleration of gravity; 

h - piezometric head in the pipe, referenced to the pipe centerline; 

ho - piezometric head in the reservoir, referenced to the pipe center 
line; 

L - length of pipe; 

v - kinematic viscosity of the fluid; 

t - time since opening of the valve; 

V - mean velocity in the pipe; and, 

x - axial distance from the inlet. 


Results obtained from ten experimental runs (Runs 20-29, inclusive) are re 
ported in this paper. The value of L/D was 380.3 in Runs 20-24 and was 
475.3 in Runs 25-29. The other externally controlled variable, ghoD3/Lv2, 
was established at values of 6(106), 12(106), 18(106), 24(106), and 30(106). 

In this paper only the phenomenon of the transition from laminar to turbi 
lent flow is considered. The specific objectives of this study are (1) to dete 
mine the position of turbulence formation, (2) to determine the manner of e: 
pansion of the turbulence, and (3) to determine the requisite flow conditions 
for the formation of the turbulence. 


Experimental Equipment and Technique 


Figure 1 is a schematic drawing of the equipment. The description of th 
equipment is presented from the upstream section to the downstream. 

The reservoir was a sealed cylindrical steel tank with a diameter of 42 
inches. The air pressure in the reservoir above the water was controlled t 
any desired magnitude by means of an air overflow line. The air overflow 
was allowed to bubble into the bottom of a standpipe. Thus by controlling tl 
water level in the standpipe, the air pressure in the reservoir could be ma: 
tained so as to obtain the desired magnitude of the reservoir piezometric 
head, hy. No water was added to the reservoir during a run, and no run we 
made without allowing the reservoir to be quiescent for at least four hours 
Since the tank cross section was large and the volume of air above the wate 
was large, the piezometric-head change was negligible (A ho/ho < 0.06%) 
during the short time required for a run. With this arrangement, the value 
ghoD3/Lv2 was determined within an estimated accuracy of + 0.4% and the 
magnitude was controlled within + 1.0% of the desired value. 

The pipe leading from the reservoir consisted of coupled sections of ext 
heavy smooth brass pipe. In order to vary the pipe length easily and 
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systematically, five sections approximately 95 diameters in length were ma- 
chined and fitted to insure close abutment and proper alignment when forced 
into the machined-brass sleeve couplings. Each pipe length was carefully 
aligned along a straight horizontal axis. 

The extruded brass pipe proved to be extremely uniform in diameter. 

_ Micrometer measurements were taken on orthogonal diameters at each end 
section. The mean diameter of these twenty measurements was 0.544 inch 
with all readings being within t 0.001 inch of this value. The average diame- 
ter of each pipe length was also checked by weighing the water that would fill 

each length. 

_ The rounded pipe inlet was located in the reservoir. The inlet piece was 

_ attached to the first pipe length by a machined connection similar to the 
Sleeve coupling. The machined brass inlet piece was made with a radius of 

rounding equal to 1.13 pipe diameters. A band of non-uniform sand roughness 
approximately 1/8-inch wide was cemented to the inlet piece slightly upstream 
from the tangent point of the inlet roundness. 

Wall piezometer holes were placed 1/2-inch from the downstream end of 
@ach pipe length. Four piezometer openings were drilled on orthogonal diame- 
ters with a No. 45 (0.0820 inch) drill. The holes were carefully inspected, 
and any burrs caused by drilling were carefully removed. The four piezome- 
ter openings were connected by a piezometer ring which was machined in the 
Pipe sleeve coupling. 

_ Pressure-time measurements were obtained by means of pressure trans- 
ducers, the signal of which was electrically amplified and recorded. These 
differential-pressure transducers were mounted on a transducer collar in or- 
der that the water passage from the pressure transducer to the piezometer 

Ting was large in cross section and short in length. The pressure transducers 

vere of the unbounded-resistance-wire type. The transducer impulses were 

mplified and recorded on an oscillograph. The oscillograph contained an in- 
lependent timing mechanism with the result that the oscillograph record con- 
sisted of the timing pips as well as the output record of the two pressure 

ansducers. Because of the cost of the transducers and the recording equip- 
nt, only two transducers with a two-channel recorder were utilized. Runs 

h four and five pipe sections (Runs 20-29, inclusive) were duplicated in or- 

er to obtain pressure-time records at all piezometer stations. 

In order to simplify the boundary conditions associated with the flow es- 

tablishment, a downstream valve was required which could be instantaneously 

pened and which in the opening process would not influence the resulting 

uid motion. The physical requirements of such a valve are (a) in a negligi- 

le time after valve release the piezometric head, h, at the pipe outlet must 

irop from the value, ho, to zero, and (b) the valve must be moved in the same 
irection and with a faster motion than the fluid particles issuing from the 
ipe outlet. The type of valve selected was a flat disk, which was positioned 
lush against the downstream end of the pipe. The force of a stressed tension 
pring was utilized to remove rapidly the valve disk from the oncoming fluid. 
Velocity determinations were made from individual frames of a 35 mm 
notion-picture camera record of the jet. The exposure time for each frame 

fas 0.0015 second. The camera was run at a speed of approximately 36 

frames per second. Rear lighting was used in order to eliminate the shadow 

of the jet. Also, a clock image was recorded in each motion picture frame. 

‘ime could be determined to the nearest 0.001 second. A horizontal and a 

ertical grid scale were also photographed. These scales were positioned in 


he plane of the jet to avoid parallax. 
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Transformation of Data into a Form Suitable for Analysis 


With this experimental equipment, the time-dependent data were available 
in the form of a 35 mm motion-picture film strip of the jet and the oscillo- 
graph records from the pressure-transducer circuit. The reduction of the 
pressure data into a useable form was straightforward. The reduction of the 
film strip data into a useable form was not straightforward. In order to ex- 
plain the process of data transformation, the data of Run 29 is used through- 
out the following description. In Run 29 all five pipe lengths were utilized, 
that is, the value of L/D was 475. The other independent variable, 
gh D3/Ly2, was established at the value, 30(106). : 


Pressure Transducer Data 


In each run the pressure data as recorded by the oscillograph was re- 
plotted in the form of the dimensionless piezometric head ratio, h/ho, as a 
function of dimensionless time, t Vgho/L. Time, t, was measured from the 
instant of opening of the valve. The very rapid change in piezometric head, 
h, from the initial value, hp, at the instant the valve was opened could be 
easily discerned on the oscillograph record. Time was measured from this 
point on the record by the pips made by the independent timing stylus. 
Figure 2 is a graph of h/hg at each of the five measuring stations as a func- 
tion of t¥gho/L for Run 29. The position of measurement is given by the 
dimensionless distance from the inlet, x/D. Runs 29 and 29a were duplicate 
runs in order to obtain pressure measurements at four stations with the two 
pressure transducers. The pressure at the outlet, x/D = 475, is zero for t>{ 


Film Strip Data 


Six frames of the film strip of Run 29 are shown in Fig. 3. 

Two assumptions were used to determine fluid velocity from the photo- 
graphically recorded jet. First, air resistance was considered negligible. 
Second, the jet was considered to be a succession of fluid masses composed 
of individual fluid particles each having the velocity characteristics of this 
larger group. If the velocity characteristics of the individual particles of a 
fluid mass are identical, the motion can be completely described by the mo- 
tion of the mass center of that fluid mass. In Fig. 4 the dimensionless “ob- 
served velocity,” V'/V 2ghg, is plotted as a function of t ¥gho/L. 

The magnitude of the error involved in assuming negligible air resistance 
can be determined by observing the same fluid mass on successive film stri 
frames. The relationships between V' and t determined at + 0.5 ft. and 
+ 2.50 ft. below the pipe were identical, indicating negligible air resistance. 

The correctness of the assumption that the jet can be visualized as a 
series of individual fluid masses is confirmed by the appearance of the jet 
(Fig. 3). In all the photographs, discrete masses of fluid appear to be formir 
in the jet. A logical explanation of why the jet can be considered to be a 
series of fluid masses is as follows. During the time the fluid is in the pipe, 
the external force of boundary shear acts on the fluid at the pipe wall. This 
boundary-shear force results in a system of internal shear forces within the 
fluid and results in a radial variation of linear momentum. The fluid par- 
ticles close to the wall are deficient in momentum. As the fluid issues from 
the pipe into the atmosphere, the external boundary-shear force becomes 
zero. However the internal shear forces are reduced to zero only if the 
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linear momentum is constant in a cross section. The internal shear forces 
persist until the jet has constant linear momentum in the cross section. Thus 
the internal shear forces are the implement for effecting the radial transfer 
of linear momentum. The result of this process is a jet with uniform linear 
momentum in a cross section which is, of course, the assumption utilized in 
determining the “observed velocity,” V'. 

However, the “observed velocity,” V', is not the mean velocity, V, of the 
Pipe flow. Because of the boundary-shear force between the fluid and the pipe 
wall the velocity distribution is not uniform within the pipe. Thus the fluid 
upon leaving the pipe has a greater linear momentum than if the velocity of 
all particles were equal to the mean velocity, V. Since the linear momentum 
is constant in a fluid mass after this fluid mass leaves the pipe, the fluid 
mass attains the velocity, V' , which is greater than V. The two velocities 
are related by the relationship, V = V' /Cm; in which Cy is the momentum 
correction coefficient. Cy) is a function only of velocity distribution and is 


defined by the relationship, C,, = (1/A) ap (v/v)? dA, in which v is the point 
A 


velocity and A is the cross sectional area of the pipe. Therefore in order to 
utilize quantitatively the data of Fig. 4 in a study of pipe flow, the velocity 
distribution had to be determined as a function of time. 

The velocity distribution, or specifically the value of Cy, in unsteady 
laminar flow was determined by a combination of theoretical analysis and 
empirical results. 

A solution of the Navier-Stokes equations was obtained in which the experi- 
mentally determined data relating piezometric head and time were introduced. 
If the pipe flow is uniform, that is, the velocity is axial, two of the Navier- 
Stokes equations are indicative that (a) the piezometric head is constant in a 
fluid cross section which is normal to the flow and (b) the axial piezometric- 
pressure gradient 0 p*/@x is only a function of time. The remaining equation 
is 


2 % 
cline Onve Yay pees) 
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in which v is the point velocity and P is the mass density. The boundary con- 
ditions are that the velocity, v, is zero at the boundary and that the velocity 
gradient, Ov/dr, is zero at the pipe centerline. The initial condition is that 
the velocity, v, is zero when time, t, is zero. 

_ The analogous equation with analogous boundary and initial conditions has 
been solved for temperature distribution in a cylindrical rod.(1) The solution 
for the pipe flow problem is 


: 
28 t 
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The various terms are as follows: J, and Jj are Bessel functions of the fir 
kind of orders zero and one, respectively; the positive numbers a, are the 
zero roots of Jo; ro is the value of the radius, r, at the pipe wall; f is the 
dimensionless radial coordinate, r/rp; Tis a dummy variable of integration; 
e(T) is - @p*/ox; and Cis yt/ro2. 

In order to apply Eq. (2) in the analysis of the experimental results, the 
existence and location of a uniform laminar flow region must be determined 
and, then, the time variation of the piezometric-pressure gradient, @p*/dx, 
in this region must be determined from the pressure-time data of Fig. 2. : 

A physical interpretation of the pressure-time data (Fig. 2) is requisite t 
locating a laminar flow region. Since the fluid starts from rest, the flow is — 
laminar in the entire pipe during the initial period after the downstream valv 
is opened. The first external indication of a change in flow regime is the 
abrupt pressure increase at the piezometers. For example, as shown in 
Fig. 2, the first of these abrupt pressure increases occurred at the location 
x/D of 189 at the time t(gho/L of 0.616. The interpretation of these pressur 
increases is that an interface between laminar flow and turbulent flow is pas 
ing the piezometer at that instant. The reason for the rapid change is that 
the mean kinetic energy of the cross section is greater in the laminar flow 
than in the turbulent flow even though the mean velocity, V, is identical on 
both sides of the interface. The difference in these kinetic energies is due t 
the different velocity distributions. Thus at x/D of 189 the flow was laminar 
prior to tVgho/ L of 0.616 and was turbulent during the remainder of the run 

The time when turbulence passes through the outlet can also be determine 
from the pressure record. In Figs. 2 and 4, the time coincidence of the rapi 
pressure drops and the rapid drop in the “observed velocity” V' is apparent 
Both occurrences are a manifestation that the laminar-turbulent interface is 
passing the pipe outlet at x/D of 475. Prior to tVgho/L of 1.38 the flow is 
laminar at the outlet. During this time a laminar-turbulent interface is mov 
ing toward the outlet. The flow is turbulent upstream from the interface. A 
explained in the preceding paragraph, the piezometric head is greater in the 
turbulent flow immediately behind the interface than in the laminar flow im- 
mediately preceding the interface. The outlet of the pipe is a station of fixe 
piezometric head since the liquid jet is projected into the atmosphere. Thus 
as the laminar-turbulent interface passes the outlet, the piezometric-head 
line is rapidly lowered just upstream from the outlet. The effect is a pivotil 
of the piezometric-head line. This rapid pivoting of the piezometric-head 
line results in a piezometric-head drop at every measuring station. The 
magnitude of this drop is directly proportional to the distance from the inlet 
This effect is clearly illustrated in Fig. 2 at the time, tVgho/ L, of 1.38. 

With this preliminary knowledge concerning the movement of turbulence, 
regions can be definitely located in which the flow is laminar. For example 
in Fig. 2 between the time, tYgho/ L, of 0.715 and 1.06 the flow is laminar be 
tween the x/D stations 379 and 475, and prior tot Vgho/L of 0.715 the flow 
is laminar between the x/D stations 284 and 475. 

The piezometric-head lines are indicative that the laminar flow regions 
are uniform. At constant values of time, t ¥gh,/L, the measured values of 
piezometric head were plotted as a function of x/D. Since each plotted 
piezometric-head line was straight at a given time in this laminar flow re- 
gion, the conclusion was that the region was also a uniform flow region as 
visualized in the derivation of Eq. (2). A contributing factor in the existencs 
of this uniform downstream laminar flow region is the presence of upstrean 
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turbulent flow regions, since an upstream turbulent flow region has the effect 
of isolating the inlet. Thus the downstream laminar flow region is similar to 
flow in an infinitely long pipe since the boundary-layer thickness, 
piezometric-head gradient, and developing velocity profile are functions only 
of time. 

The slopes of the piezometric-head lines of the uniform laminar flow were 
determined and plotted. These slope values are shown as points on Fig. 5 for 
the run being discussed. In order to incorporate this experimental data into 
Eq. (2), an analytic function for g(T) must be selected that approximates the 
empirical results. A very simple function was selected, that is, 0 p*/Ox is 
proportional to cos(™t/2t*) in which t*/gho/L is empirically defined as shown 
on Fig. 5. The value of t*\/ghj/L is dependent upon the experimental parame- 
lers as well as the time and position of turbulence inception. The function for 
B(T) was introduced into Eq. (2) and the integration performed. The resulting 
expression for the time-developing velocity profile is 


an? ot (+/t*) 
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All terms of Eq. (3) have been previously defined except a * which is t*/ro2. 

Values of v/V were computed in systematic steps so that by interpolation 
of the computed results, these results could be applied to all experimental 
runs. The first nine zero-roots of the Bessel function were used in the nu- 
merical computations. The range of the other variables were as follows: 
N<B<1, 0.01< o *<0.05, and 0.2<t/t*< 1.0. Values of the linear momentum 
correction coefficient, C;,, were then computed using numerical-graphical 
methods. Hence to determine the temporal variation of Cy for any experi- 
mental run only the value of t*Vgho/ L was required. 

The analytical solution does not cover the complete temporal range of 
laminar flow. An empirical method was used to extrapolate the analytical 
solution. The method is based upon the ‘observed velocity” data. At the time 
vf transition from laminar flow to turbulent flow at the pipe outlet, there is a 
sharp discontinuity in the “observed velocity,” as shown in Fig. 4. However 
here can be no discontinuity in the mean velocity. Thus the observed dis - 
sontinuity must be due to the difference in Cy, of the laminar flow and turbu- 
ent flow. Hence at the time of transition 


', — 
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ince (Cm)turbulent is equal to 1.02. The value of Cy of 1.02 for turbulent 
w is based upon the seventh-root law of velocity distribution. 
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The mean velocity, V, was then computed from the “observed velocity,” 
V'. The displacement boundary-layer thickness was also computed using t 
following relationship 


b 3 
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in which v¢ is the centerline velocity and 6* is the displacement boundary- 
layer thickness. The values of V/V 2gho, 5*/r,g, and v,/V of Run 29 are p 
sented in Fig. 6 as a function of time. 


Analysis of Results 


The specific objectives of this paper pertain to the determination of the 
manner of turbulence inception and the manner in which the turbulent flow 
zones are expanded until the pipe is filled with turbulently flowing fluid. 


Downstream Movement of Turbulence 


The pressure-time records, of which Fig. 2 is an example, of all the run 
are indicative that the turbulent flow region is moved toward the outlet in a 
systematic manner. As previously explained, the abrupt pressure increase 
at the measuring station is an indication that the laminar-turbulent interface 
is passing the piezometer, and the simultaneous pressure decrease at all 
measuring stations is an indication that the laminar-turbulent interface is 
passing the pipe outlet. A systematic movement of the turbulence is indicat 
in Fig. 2 between stations at which x/D is 284 and 379 and between stations” 
379 and 475. The mean transport velocity, V*, (between the measuring sta- 
tions) of the laminar-turbulent interface was computed. The results are pre 
sented in Table I. These results are indicative that the downstream move- 
ment of turbulence is a physical transport phenomenon in which the turbuler 
eddies are transported downstream with a velocity equal to the centerline 
velocity of the preceding laminar flow. 


Nature of the Downstream Face of the Turbulent Zones 


A qualitative picture of the nature of the downstream face of a turbulent 
zone can be constructed from the occurrences recorded as the interface is 
moved past the pipe outlet. The appearance of a typical jet shortly after the 
interface has passed the outlet is shown on Fig. 3 (e). The simultaneous 
pressure decrease at all measuring stations as shown in Fig. 2 occurs duri 
the time the interface is passing the outlet. The rapid decrease in the “ob- 
served velocity,” V', as shown in Fig. 4 also occurs when the laminar- 
turbulent interface passes the outlet. 

In order to establish the chronological order of events, the times of vari 
ous occurrences are presented in Table II. Columns 2, 3, and 4 of Table II 
are based upon motion picture data. In columns 2 and 3 are recorded times 
of passage through this outlet of the last laminar flow and of the first turbu- 
lent flow, respectively. The fluid that passed the outlet in the interveing tir 
is presumably in the interface region. Based upon this time increment the 
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computed lengths of the interface region are presented in column 4. These 
computed lengths are useful only to establish the order of magnitude of the 
length of the interface region because the time interval between the last lami- 
nar flow and the first turbulent flow is too short to be determined with pre- 
cision. The time of passage of the last laminar flow through the outlet as in- 
dicated by the pressure record is presented in column 5. The time of passage 
of the first turbulent flow through the outlet as indicated by the pressure 
record was not tabulated because the rate of pressure decrease as indicated 
on the record was of the same order of magnitude as the capability of the 
oscillograph to follow rapid changes. The interpretation was that the rapid 
pressure drop was at least as rapid as indicated and was probably more 

rapid than indicated. 

Based upon the previously determined transport velocity and the recorded 
jata (Table II) at the time the laminar-turbulent interface passed through the 
pipe outlet, a possible configuration of the interface can be conceived which is 
in basic agreement with the data. Fig. 7 is a schematic drawing of the inter- 
face region. The flow is from left to right with laminar flow preceding the 
urbulent flow. The turbulence in region (a) of Fig. 7 is transported down- 
stream at a velocity, V*, equal to the centerline velocity, vce, of the laminar 
low. The turbulence diffuses radially outward from the central region (a) in- 
0 the boundary layer (b). In this advanced stage of laminar boundary-layer 
levelopment, the boundary layer is unstable for a wide range of disturbance 
Trequencies. As the turbulence is diffused into the laminar boundary layer, a 
joint is reached at which the entire laminar boundary layer (c) suddenly be- 
-omes turbulent. The last indication of laminar flow from the pressure re- 
ord (Table II, Column 5), most closely corresponds with the first indication 
4 turbulent flow from the motion picture data (Table II, Column 3). The wall 
iezZometric head, h, probably does not rapidly change (d) until the entire 
aminar boundary layer becomes turbulent (c). The deceleration of the 
locity in the central region as the turbulence penetrates the boundary layer 
b) is the cause of the radial piezometric head gradient, (e) to (d), necessary 
0 effect the radial mass transfer involved in changing from the laminar 
elocity distribution to the turbulent velocity distribution. 


Yelineation of the Turbulent Flow Regions 


5 The spatial and temporal location of the turbulent regions within the pipe 
an be determined from the piezometric-head lines. The details of this 
rocess are illustrated with the data of Run 29 (Figs. 2 and 6). The values of 
he piezometric-head ratio, h/ho, were determined from Fig. 2 at each of the 
aeasuring stations at constant values of the time parameter, ty gho/L. Since 
he time that a laminar-turbulent interface passed a measuring station could 
e determined from the pressure record and Since this interface was moved 
‘ith the centerline velocity of the laminar flow, the position of the downstream 
xtremity of a turbulent region could be calculated at the constant values of 
ie time parameter for which the piezometric-head lines were to be plotted. 
‘he piezometric-head changes at the laminar-turbulent interfaces were cal- 
ulated as being the difference in the mean kinetic energy of the laminar and 
irbulent flow. The piezometric-head lines were then plotted as shown in 
ig. 8 A. 

A rather complex plotting procedure was adopted in Fig. 8 in order to show 
oth the temporal and spatial location of the turbulent flow regions. The 
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constant value of t V¥ gh,/L for each piezometric-head line is shown at the | 
right hand margin. The position of the zero value of h/hy was lowered wi 
increasing values of t V gho/L in order to avoid a superposition of many 
piezometric-head lines. Since the plotted points at the largest values of x 
are for the pipe outlet, the zero value of each piezometric-head line is at 
this level. With the known positions of the downstream extremity of the t 
bulent regions and with the slopes of the piezometric-head lines of the la 
nar and turbulent flow, the upstream extremity of the turbulent regions c 
be graphically determined. The regions of turbulent flow are clearly deli 
ated in both time and space on Fig. 8. The turbulent flow regions are bel 
the dashed line and the laminar flow regions are above the dashed line. A 
matter of incidental interest, the last laminar flow of Run 29 (Fig. 8A) at- 
tained a pipe Reynolds number, VD/V, of 5.1 (10)4. 


Inception of Turbulence 


With the turbulent flow regions located both spatially and temporally as 
shown in Fig. 8, the flow conditions at the time of origin of the turbulent s; 
can be qualitatively analyzed with reference to the stability of small distuz 
ances.(2) Basically, the method is to introduce a two-dimensional periodi 
disturbance into the Navier-Stokes equations and to determine whether the 
disturbance is amplified or diminished as the disturbance is moved down- 
stream. The neutral curve which separates the regions of amplification a 
decay is presented on Fig. 9 for Blasius flow over a flat plate.(3) The ord 
nate value is a ratio of the displacement boundary-layer thickness §6* to tl 
disturbance wave length X. The abscissa value of Fig. 9 is the boundary- 
layer Reynolds number R §*. 

The salient features of turbulence inception are apparent in Fig. 8, and 
pertinent flow parameters are presented in Table II]. The most prominen 
characteristic is that turbulence started as spots in laminar flow. The 
downstream face of these spots was transported downstream in a systema 
manner as shown previously. The upstream face of the turbulent spots di 
not indicate a systematic movement along the pipe; but even for the spots 
which the upstream face moved downstream, this movement was at a less 
velocity than for the downstream face. The result is a joining together of 
spots as shown clearly in Fig. 8A. The pattern of turbulent spot generatic 
was approximately the same for equal values of the independently controll 
parameter, ghoD’/Lv2. That is, in Runs 29 and 24 the value of ghoD3/Lv 
was 30(10)6, in Runs 28 and 23 the value was 24(10)6, in Runs 27 and 22 th 
value was 18(10)6 in Runs 26 and 21 the value was 12(10)6, and in Runs 25 
and 20 the value was 6(10)®. 

Since the flow throughout the pipe was laminar prior to turbulent spot 
generation, the initial transition from a laminar boundary layer to a turbu 
boundary layer must have been aided by an instability of the laminar boun 
layer rather than the external source of free-stream turbulence. The ma 
tude of R§*, v.6*/y, at turbulence inception (Table Il) in comparison wit 
the magnitude of R§* required for boundary-layer instability (Fig. 9) is ix 
dicative that turbulence exists soon after the boundary layer is unstable. 
Direct numerical comparison between the values of R§* in Table II and 
Fig. 9 is not possible for two reasons. Firstly, the curve of neutral stabi 
of Fig. 9 is for Blasius flow over a flat plate rather than for a temporally 
growing boundary layer in a circular pipe. Secondly, the experimentally 
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determined values of R§* of Table II are the values at the time of turbulence 
inception; whereas, the values of R §6* of Fig. 9 are the values at which dis- 
turbances are either amplified or diminished. Since a disturbance undoubted- 
ly must undergo some amplification before the disturbance is transformed in- 
to turbulence, turbulence inception occurs at a later time and at a farther 
dewnstream position than the time and place of disturbance generation. 

The location of the points of turbulence inception (Fig. 8 and Table III) is 
indicative of the source of the disturbances. 

Turbulent spot number 1 occurred in all runs and in the position range of 
0< x/D< 62. The source of disturbances for turbulent spot number 1 was 
ithe band of sand grain roughness which was cemented to the bellmouth inlet. 
Since the source of disturbances was at the inlet, the value of x/D at the posi- 
tion of turbulence generation is evidence of the distance required to amplify 
the disturbance before turbulence is actually generated. The flow pattern is 
complex in the region just downstream from the inlet since the boundary layer 
is expanding both with increasing time and increasing distance from the inlet. 
Because the disturbance is generated in a region of small boundary-layer 
thickness, the disturbance may first diminish in amplitude due to the small 
and hence stable value of R§*. Since the boundary layer increases in thick- 
ness both with distance and time, the disturbance which was initially diminish- 
ing may pass into the region of instability or amplification. The disturbance 
is moved further downstream during amplification. At some point the disturb- 
ance attains a sufficient amplitude to break into turbulence. The heterogene- 
ous mixture of disturbance wavelengths from the nonuniform sand grain 
roughness as well as the spatially varying laminar boundary-layer thickness 
precluded the possibility of avoiding the unstable region (Fig. 9) with the re- 
sult that turbulent spot number 1 existed in all runs. 

Turbulent spot number 2 occurred in all runs except Run 25 and in the 
~0Sition range 99< x/D < 134. This position is indicative that the source of 
listurbances was the pipe junction at x/D of 95. In spite of the care of pre- 
Jaring the pipe junctions, the pipe edges could have been mismatched as much 
is 0.001 inch because of diameter variations between the ends of each pipe 
ength. With the pipe junction as the source of disturbances for turbulent spot 
lumber 2, the distance required for amplification of the disturbance was from 
{to 39 pipe diameters. More significantly, this distance for amplification 
aries between 1306* and 6506*. In order to determine the minimum possi- 
le wavelength the maximum value of 277§*/) of 0.409 is utilized from Fig. 9. 
thus the minimum wave length,A, is approximately 156 *. With this minimum 
vavelength and the distance for Bmplifieation, the maximum number of dis- 
urbance waves between the disturbance source and the point of turbulence in- 
eption would be approximately 40. 

_ The existence of turbulent spot number 2 (Table Ill) in all runs except Run 
5 can be qualitatively explained by reference to the neutral stability curve. 
the dimensionless boundary-layer thickness, 6 */ro, and the boundary-layer 
teynolds number, R6*, were computed for each run. In Fig. 10 the value of 
*/r 0 is presented as a function of R§*, vc 6*/V for the early portion of all 
uns. The boundary-layer parameters as shown in Fig. 10 are very similar 
or Runs 29 and 24, Runs 28 and 23, and so forth. This similarity of boundary- 
ayer parameters is the reason for ie similarity of the turbulent spot pat- 
rns which was mentioned previously. With the hypothesis that the disturb- 
nce wave length A is proportional to the roughness height at the pipe junc- 
ion, a possible neutral stability curve for turbulent spot number 2 is shown 
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on Fig. 10. This neutral stability curve is for a disturbance wavelength, A, ‘ 
approximately 2.1 rg. This neutral stability curve is for the longest disturb. 
ance wavelength which originates at the pipe junction. Thus, the neutral 
stability curves for shorter wavelengths would be below the curve shown in 
Fig. 10. With this concept there is no conflict in that some of the points of 
turbulence inception (indicated by the circular symbols on Fig. 10) are below 
the lower branch of the curve shown. The curve of neutral stability is a line 
of zero amplification. Other lines of constant amplification are similar in ~ 
shape to the neutral curve and resemble the contour lines of the end of a : 
ridge on a topographic map. Since the boundary layer is undergoing essen- 
tially temporal growth, the boundary-layer curves of Fig. 10 can be visua- 
lized as pertaining to a single wave or disturbance. Thus a disturbance of 
Runs 20 and 25 (Fig. 10) would be subjected to lesser amplification in cross- 
ing the unstable zone than for the other runs. The disturbance of Run 20 re- 
ceived sufficient amplification to break into turbulence, whereas the disturb- 
ance of Run 25 passed into the upper stable zone before breaking into 
turbulence. Once a disturbance is in the upper stable zone, this disturbance 
will diminish without breaking into turbulence. 

Turbulent spot number 3 occurred in only Runs 24 and 29 at the x/D posi- 
tion of 200 and 213. This position is indicative that the source of the disturb 
ances was the pipe junction at x/D of 190. This axial distance is equivalent t 
3106 * and 6006 *. The maximum number of disturbance waves between the 
source and the point of turbulence would be approximately 40. A possible 
neutral stability curve for the longest wavelength, equal to 1.5 ro, is shown i 
Fig. 10 for turbulent spot number 3. The interpretation is the same as befor 
except that only in Runs 24 and 29 was a disturbance amplified sufficiently tc 
break into turbulence. In all other runs a disturbance was not sufficiently 
amplified in crossing the unstable zone to break into turbulence. 


Upstream Face of a Turbulent Spot 


The experimental data are of little use in defining the nature of the up- 
stream face of the turbulent zones. Since every individual turbulent eddy is 
transported downstream, new turbulence would have to be generated even fo! 
the upstream face of the turbulent zone to remain in a fixed position. Since 
the position of the upstream face of the turbulent zones depends upon the 
generation of new turbulence, the probability is that this face has an erratic 
motion. If two turbulent spots are joined before the downstream face of the 
upstream spot passed a pressure transducer station, the presence of the up- 
stream spot could not be detected from the data. In fact, it is probable that 
turbulent spot number 1 of Fig. 8B and Fig. 8D is formed from two or more 
spots which are soon joined. The upstream face of the turbulent zone would 
move downstream if no new turbulence were generated. In a temporally 
growing boundary layer, a single burst of turbulence might be generated, sul 
sequent disturbances from a fixed source of disturbance being insufficiently 
amplified to break into turbulence. Referring to Fig. 10, this occurrence ap 
pears to be probable for turbulent spot number 2 in Run ’20 and for turbulent 
spot number 3 in Runs 24 and 29. The definite downstream movement of the 
upstream face of these turbulent spots, as shown in Figs. 8A, 8B, and 8J, is 
indicative of this occurrence. 
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CONCLUSIONS 


Based upon an experimental study of unsteady flow, the following conclu- 
sions concerning transition from laminar to turbulent flow ina pipe are con- 
sidered to be significant: 


_ 1. In an experimental investigation of laminar boundary-layer transition 
he effects of free-stream turbulence can be avoided entirely with a temporal- 
y expanding boundary layer. 


2. Transition from laminar to turbulent flow in a pipe occurs in spots. 


3. The magnitudes of the laminar boundary layer Reynolds number at 
urbulence inception was in reasonable agreement with the theoretical values 
letermined for Blasius flow over a flat plate. 


4. The downstream movement of a turbulent zone into a laminar zone is 
hat of physical transport. The transport velocity is that of the centerline 
‘locity of the laminar flow. In the region of this interface the free stream 
urbulence is the disturbance source for laminar boundary layer transition. 


5. The regions of laminar and turbulent flow in a pipe can be defined from 
he piezometric-head lines. 


6. From the standpoint of boundary shear, the pipe utilized in the experi- 
nent would be classified as smooth. However, from the standpoint of laminar 
oundary-layer transition, disturbance sources were found at the pipe junction 
S well as at the sand-roughened inlet. 


7. With a systematic variation of boundary-layer parameters, the presence 
r absence of boundary-layer transition at the pipe junctions was found to be 
n qualitative agreement with the laminar boundary-layer stability theory. 
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Figure 8, Piezometric Head Lines of Runs 20-29, Inclusive. (Cont inued) 
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Figure 8. Piezometric Head Lines of Runs 20-29, Inclusive. (Continued) 
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Figure 8. Piezometric Head Lines of Runs 20-29. Inclusive. (Continued) 
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Table I 


% 
Downstream Transport Velocity of Turbulence 


aa west le at wescay (elt 
From To Velocity v/v, 
2la 189 28k 0.97 
22 284 380 0.99 
22a 189 380 0.96 
23 28), 380 1.00 
23a 189 380 0.98 
23b 28) 380 1.01 
2h 28) 380 0.97 
27 189 379 0.97 
27b 189 379 1.01 
27c 189 379 1.00 
27d | 189 379 1.02 
27f 189 9) 1.03 
28a 189 379 1.00 
29 28) 475 1.01 
29a 319 475 1.02 


“Data from Run Numbers 20, 25 and 26 as well as a portion of the 
data from Run Numbers 21, 27, and 28 have been omitted from this tabula- 
tion. The omitted data occurred during the time the jet was separating 
into individual fluid masses having different values of linear momentum. 
In other words, the radial transfer of linear momentum by the internal 
shear forces was incomplete and some fluid masses fell from the jet 
with a lesser momentum than the main body of the jet, The calculated 


values of the centerline velocity Vv, were questionable during this type . 
of jet separation. 
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Table II 


Laminar - Turbulent Interface at the Time of 
Passage Through the Pipe Outlet 


x. 4 me ae cow one: Bate Data Hoe Ven Record 
Interface 
Last Laminar First Turbulent Last Laminar 
20 1.532 1.566 9D 1.587 
21 1.447 1.493 11D 1.501 
22 1.386 1.479 23D 1.468 
23 1.433 1.462 7D 1.460 
2),* ay ie as 1.214 
25 1.781 1.835 17D 1.844 
26 Teeiesp 1.553 7D 1.583 
27 1.526 1.567 13D 1.579 
28 1.530 1.567 11D 1.584 
29 1.333 1.365 9D 1.386 


— 


¥ 
' *the length of the interface region was very short with the result 
pret meaningful determinations of the time of last laminar flow passage 
and of the time of first turbulent flow passage were impossible. 


: 
rears! OFan tO See ee 


f 


1450-30 | ) HY 6 December, 1957 | 


Table III 


Flow Parameters at the Inception of Turbulence 


Run No. Turbulent Position of Boundary Boundary 
Spot Turbulence Lay! r 


Number poceene ee Thickness Reynolds Number 
in Diameters S /r, Re aes Ye oie 4 
20 112 0.100 1140 ‘ 
va 13h 0.120 - 1810 + 
21 1 te) 0.078 1030 
2 99 0.062 Lio “a 
22 1 62 ~ 0.070 1100 me 
2 105 0.062 710° 
23 1 48 0.068 1020 
2h 1 37 
Bans 2 100 
ae OO iia 
25 eee ue 
of vve.e ¥ a ‘ ¥ 
: ae 
se “Says? ae viet nis 
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ABSTRACT 


The feasibility of using rainfall data to obtain flood frequencies has been 
demonstrated before. Most investigations, however, utilized rainfall data and 
rainfall-runoff relations specific to the basins for which the flood-frequency 
data were desired. The present study, involving 36 basins, shows that the 
rainfall data need not be from the problem basin and that the data from one 
precipitation network can be used to synthesize streamflow data for a number 
of basins. Similarly, rainfall-runoff relations already derived for other 
basins may be easily modified to apply to the problem basin. The transpo- 
Sition of rainfall data and rainfall-runoff relations not only make it possible 
to extend the streamflow record of basins not having rainfall data, but greatly 
reduces the cost of synthesizing streamflow data. 


INTRODUCTION 


Among the procedures for extending streamflow records for design purpos- 
es is that reported by Kresge and Nordenson, 1) who demonstrated that 
streamflow data, specifically annual flood peaks, could be synthesized from 
rainfall data by means of a rainfall-runoff relation and unit hydrograph. Us- 
ing the Hocking River at Athens, Ohio, area 944 sq. mi., which had a 4-station 
recipitation network and a 37-year period of concurrent precipitation and 
streamflow records, they showed that the flood-frequency characteristics of 
he basin could be reproduced from rainfall data. The procedure requires a 
reat deal of work and the cost of application on a broad scale might be con- 
sidered prohibitive, especially in connection with the design of low-cost hy- 
lraulic structures, such as highway culverts. The results were so encourag- 
ng, however, that in 1954-56 the U. S. Weather Bureau and the Bureau of 


jote: Discussion open until May 1, 1958. Paper 1451 is part of the copyrighted 
_ Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 6, December, 1957. 
. Chief, Cooperative Studies Section, U. S. Weather Bureau, eee DIC: 


, U. S. Weather Bureau, Washington, D. C. 
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Public Roads jointly undertook further investigations along these lines to de-| 
termine if the procedure could be modified to lower the cost of application 
when used on a broad scale. 

A preliminary study based on three basins: Delaware River at Valley 
Falls, Kans. (922 sq. mi.); Meramec River near Steelville, Mo. (781 sq. mi.); 
and Ralston Creek at Iowa City, Iowa (3 sq. mi.), showed that their flood- 
frequency characteristics could be satisfactorily reproduced from rainfall 
data outside the basins. The tests were made by taking rainfall data from a 
number of networks around each basin and comparing the results of the 
frequency analysis of the 30 years of annual flood peaks at each basin with 
those of the annual flood peaks synthesized from the rainfall data for the 
corresponding networks for the same period as the streamflow record. It 
was found that a precipitation network would yield satisfactory results provid 
ed (1) its area roughly approximated the size of the problem basin, (2) its 
rainfall-intensity frequency characteristics were similar to those of the 
basin, and (3) the seasonal distribution of rainfall was the same as that of the 
basin. 

These preliminary studies thus led to the large-scale undertaking of synth 
sizing streamflow data for 36 basins (Table 1). These 36 basins, ranging in 
size from 2.2 to 521 sq. mi., are scattered in a long, narrow region having 
homogeneous physiographic characteristics* and extending from northern 
New Jersey to northern Georgia (Fig. 1). Each basin has at least 15 years of 
streamflow record. 


Procedure 


Selection of Precipitation Networks 


As it had been decided to synthesize streamflow data for no less than a 50 
year period for each problem basin, the basic precipitation networks were 
made up of stations with at least 50 years of concurrent rainfall data. This 
restriction severely limited the number of networks available. Because of 
the relatively wide scatter of stations with 50-year rainfall records, the indi: 
vidual networks were limited to no more than 3 stations, since a larger 
number would have made the network area incompatible in size with the 
problem basins. 

The effect of the size of the network area on synthesized peak discharges 
was investigated. For this purpose the West Fork River Basin at Clarksburg 
W. Va., 384 sq. mi., was used. This basin contains a fairly evenly distribute 
precipitation network of 8 stations, and is surrounded by other stations not 
far outside its boundaries. Fig. 2 shows the equation used to obtain an index 
of the effective area of the 1-, 2- and 3-station networks selected for this 
investigation and the effect of network area on synthesized 2-year peak dis- 
charges. The 2-year peak discharge synthesized from the storm-precipitati 
averages of the basin’s 8-station network served as the control. The curve 
of Fig. 2 indicates that for network areas of about 40 to 600 sq. mi., a differ: 
ence up to 350 sq. mi. between network and basin areas affects the results by 
no more than 5 percent. The range of network area from 0 to 40 sq. mi, 


*As indicated by the map of Problem Areas in Soil Conservation, published 
by the U. S. Soil Conservation Service in July 1950. Most of the problem 
basins were in the area designated as the Piedmont Plateau, or B-10, regio 
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Table 1 


PROBLEM BASINS 


STREAM AND GAGING STATION 


Whippany River at Morristown, N. J. 


North Branch Raritan River nr Far Hills, N. J. 


Walnut Brook nr Flemington, N. J. 
Neshanic River nr Reaville, N. J. 
Tohickon Creek nr Pipersville, Pa. 
Perkiomen Creek at Graterford, Pa. 
Chester Creek nr Chester, Pa. 
Codorus Creek at Spring Grove, Pa. 
Octoraro Creek nr Rising Sun, Md. 
Piney Run nr Sykesville, Md. 
Linganore Creek nr Frederick, Md. 
Seneca Creek at Dawsonville, Md. 


Northwest Branch Anacostia River nr Colesville, Md. 


Goose Creek nr Leesburg, Va. 

Slate River nr Arvonia, Va. 

Tye River nr Lovingston, Va. 
Appomattox River at Farmville, Va. 
Goose Creek nr Huddleston, Va. 

Sandy River nr Danville, Va. 

Mayo River nr Price, N. C. 

Fishing Creek nr Enfield, N.C. 

Dial Creek nr Bahama, N. C. 

Eno River at Hillsboro, N. C. 

Horsepen Creek at Battleground, N. C. 
East Fork Deep River nr High Point, N. C. 
Deep River nr Randleman, N. C. 

Fisher River nr Copeland, N. C. 

South Yadkin River nr Mocksville, N. C. 
Uwharrie River nr Eldorado, N. C. 
Brown Creek nr Polkton, Ne C. 

South Tyger River nr Woodruff, S.C. 
South Tyger River nr Reidville, S.C. 
Middle Tyger River at Lyman, S. C. 
Apalachee River nr Buckhead, Ga. 
Little Tallapoosa River at Carrollton, Ga. 
Tobesofkee Creek nr Macon, Ga. 
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however, encompasses a rapid change in selectivity, suggesting that synthe- 
sized peak discharges should be adjusted when there is any appreciable differ 
ence between basin and network area within this range. The data of Fig. 2 
also suggest that 2-station networks yield results as satisfactory as those 
obtained from 3-station networks. Use of the 2-station instead of the 3-statio 
network results in a great saving of time in making storm surveys and de- 
termining time distribution of rainfall. Best results for very small basins, 
i.e., under about 15 sq. mi., were obtained from single stations. This is true, 
of course, only because closely-spaced long-record precipitation stations 
were not available. Another important consideration in the selection of pre- : 
cipitation networks was their geographic proximity to groups of problem | 
basins, the objective being to effect economy by applying one network to a 
group of basins. 


Storm Surveys 


After the precipitation networks had been selected, their records of daily 
precipitation were scrutinized to select the major rainstorms of each water 
year. A base value of daily precipitation was arbitrarily selected (usually 
one inch) and whenever one or more stations in the network exceeded the base 
value, precipitation was tabulated for all stations for the entire storm period, 
including days with less than the base value. If fewer than 10 or 12 storms 
were obtained in any one water year, the base was lowered and the record re- 
examined. 

The runoff for each major rainstorm was then roughly estimated by means 
of a rainfall-runoff relation (Fig. 3) as described by Kohler and Linsley. 

The four or five storms producing the highest runoff in any one water year 
were then set aside for a detailed analysis of the time distribution of the rain- 
fall, because the relatively small size of the problem basins required a break- 
down of rainfall into increments no longer than 6 hours for accurate estimates 
of runoff and computations of peak discharge. 


Time Distribution of Rainfall 


The time distribution of rainfall in the several storms selected for each 
year was determined from the records of a recording-gage station in or near 
the network involved. The hourly data were scanned to determine the maxi- 
mum 6-hour rainfall increment for a storm, and the entire storm was divided 
into 6-hour increments synchronizing with the maximum. The average storm 
rainfall for the network was distributed into 6-hour increments in accordance 
with the indicated distribution for the recording-gage station. Six-hour incre- 
ments of runoff were then obtained from a rainfall-runoff relation. For basin: 
smaller than about 15 sq. mi., the peak discharges generally result from 
short-duration, high-intensity rainfall. Consequently, 1-hour increments of 
rainfall were used for such basins, the 1-hour runoff increments being esti- 
mated by means of a special rainfall-runoff olettanec 


Rainfall-Runoff Relations 


Rainfall-runoff relations were available for very few of the problem basins 
Development of these relations is a time-consuming procedure, In order to 
keep the cost of the study to a minimum, available rainfall-runoff relations 
from other basins were utilized. In order to determine which of several 


| 
ASCE PAULHUS - MILLER 1451-5 


relations was best suited to a particular problem basin, the runoff for 25 or 
more storms over the basin was computed from the hydrographs. The run- 
cif for the same storms was also estimated by means of the rainfall-runoff 
relation being tested, and plotted against the corresponding observed runoff. 
The relation yielding values nearest to a 45° line was the most representative 
of the rainfall-runoff characteristics of the problem basin. The mean curve 
for the points yielded by this most representative relation was used to adjust 
all runoff values obtained from that particular rainfall-runoff relation for use 
on the problem basin involved. Fig. 4 presents such a mean curve indicating 
so high a degree of representativity for the rainfall-runoff relation that no 
adjustment was considered necessary. 


Unit Hydrographs 


Unit hydrographs of 1-hour or 6-hour duration were derived from discharge 
data obtained from the original records of the U. S. Geological Survey for 
converting runoff increments into peak discharge. 4) Itis generally recog- 
nized that one unit hydrograph for a basin may net properly reproduce the 
stream hydrograph for all magnitudes of floods. Different unit hydrographs 
are often required for low and high runoff, and sometimes for intermediate 
values. No matter what criteria are established for the choice of the proper 
unit hydrograph, situations will exist where almost equal values of runoff will 
require different unit hydrographs, probably resulting in appreciable differ- 
ences in the peak discharges obtained. In order to avoid such discontinuities, 
Only one unit hydrograph was constructed for each basin, and a curve was 
then developed to adjust the peak discharges yielded by the unit hydrograph. 

The unit-hydrograph adjustment curve for a particular basin was derived 
by first obtaining from observed flood hydrographs the total storm runoff for 
a good range of peak discharges. The next step was the breakdown of storm 
runoff into 1- or 6-hour incremental runoff on the basis of the chronological 
distribution of the rainfall producing the runoff. The peak discharge was 
then computed by applying a unit hydrograph like that of Fig. 5 to the runoff 
increments. The computed peak discharge for each storm was then plotted 
against the observed peak discharge for the storm. The mean curve (Fig. 6) 
was used to adjust the synthesized peak discharges estimated by the unit 
hydrograph. 

_ In addition to eliminating the need for several unit hydrographs for one 
basin, the peak-discharge adjustment curve also allowed for base flow. Base 
flow was not included in synthesizing the peak flow from the unit hydrograph, 
but it was included in the observed values of peak discharge used to derive 
the adjustment curve. Consequently, a part of the adjustment provided by 
application of the curve was for base flow. 


Results 


. 

7 

| 
Comparative Test 
| 


After selection of the most representative rainfall-runoff relation and unit 
oo. and their adjustment curves for a particular basin, peak discharg- 
S corresponding to the rainstorms for the period of discharge measurements 
were estimated. The estimated or synthesized annual peak discharges, i.e., 
the maximum for each water year, were arrayed in decreasing order of 
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magnitude and then compared with the observed, similarly arrayed. If the 
average difference between the two was less than 20 percent, the record of 
synthesized peak discharges was considered to be representative of the 
streamflow record. Annual peak discharges were then synthesized for the re 
mainder of the precipitation record. 


Synthesized Annual Peak Discharges for Codorus Creek, Pa. 


Table 2 shows such a comparison between synthesized and observed annual 
peak discharges (annual floods) for Codorus Creek Basin above Spring Grove, 
Pa., area 74.3 sq. mi. Codorus Creek has a streamflow record starting with 
water year 1933, and annual flood peaks were available through water year 
1954. Peak discharges were therefore synthesized for the corresponding 
period. 

The peak discharges for Codorus Creek were synthesized from major rain- 
storms observed by a 3-station network (Ephrata, Lancaster and Lebanon, 
Pa.) about 45 miles northeast of the basin. The time distribution of the rain- 
fall in the selected storms subsequent to 1940 was determined from either or 
both of the recording gages (Lancaster, Lebanon) and prior to 1940 from the 
recording gage at Harrisburg, about 30 miles west of the network. The run- 
off from these storms was estimated by means of a rainfall-runoff relation 
previously developed for the Anacostia River Basin, about 50 miles south of 
the problem basin. This rainfall-runoff relation was so well adapted to the 
problem basin that no adjustment was considered necessary (see Fig. 4). 

The 6-hour unit hydrograph (Fig. 5) used to estimate the peak discharge 
from the computed runoff was based on five storms of approximately 6 hours’ 
duration in the problem basin. An adjustment curve was required (Fig. 6). 

The first three columns of Table 2 show poor agreement on an annual basis 
between observed and synthesized annual floods. This is not surprising as 
it would indeed be strange if floods of the same approximate frequency oc- 
curred simultaneously over the problem basin and over the network 45 miles 
away. Tests based on two pairs of adjoining basins showed poor agreement 
between observed annual peak discharges for adjoining basins on a year-to- 
year basis. 

The last two columns of Table 2, showing observed and synthesized annual 
floods in decreasing order of magnitude, suggest that the two may well repre- 
sent samples of the same population. The synthesized annual floods are 
generally well within the 20-percent limit of acceptability. A frequency analy 
sis of the observed and synthesized annual floods for the 22-year period by 
the Gumbel method yielded the results shown in Table 3. 

The frequency data of Table 3 might be criticized by some statisticians 
because they were based on all annual floods for the 1933-54 water years. 
Referring to Table 2, it is seen that the 1933 observed and synthesized annual 
floods were both extremely high, being about twice the second highest. Some 
would argue that the 1933 annual floods are not of the same population as the 
lesser floods and should therefore not be included in the frequency analysis. 
Table 4 shows the results of the frequency analysis of observed and synthe- 
sized annual floods for water years 1934-54. This table shows even better 
agreement between observed and synthesized annual floods than when the 
1933 data were included. In any event, it appears that the record of synthe- 
sized annual floods compares reasonably well with that of observed annual 
floods for Codorus Creek Basin. ‘ 
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Table 2 
ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA. 
(74.3 Sq. Mi.) 
CHRONOLOGICAL LISTING STATISTICAL ARRAY 
Water Observed Synthesized Order Observed Synthesized 
Year (cfs) (cfs) No. (cfs) (cfs) 
1954 1910 2890 1 11200 11600 
1953 2180 22ho 2 6070 5370 
1952 1750 2340 3 4880 4900 
1951 3200 4900 k 4180 4800 
1950 1610 2590 5 3200 3910 
i949 1730 1640 6 3180 3670 
1948 1810 1930 7 3070 2890 
1947 1430 1320 8 2710 2810 
1946 3070 4800 : 9 2490 2730 
1945 2490 2810 See LO 2380 2590 
1944 2380 3910 ei 2210 2360 
1943 4180 1750 Sad 2 2180 2340 
1942 2710 5370 -— 13 1910 22k0 
1941 870 1340 : I 1810 1970 
gh 4880 1970 ae 1750 19h0 
1939 1620 1940 3) P16 1730 1930 
1938 3180 2730 os S027. 1620 1780 
S 1937 2210 2360 ie ad 1610 1750 
1936 1380 1780 : 19 1430 1640 
1935 1110 1330 = iS) 1380 1340 
1934 6070 3670 sao: 1110 1330 
| 1933 11200 11600 : 22 870 1320 
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The suitability of the above precipitation network, rainfall-runoff relation 
and unit hydrograph for synthesizing peak discharges for Codorus Creek 
Basin having been established, the next step was to synthesize annual floods 
for the remainder of the precipitation record. The results are given in Tabl 
5. 

The results of a frequency analysis of the 55-year period of synthesized 
annual floods of Table 5 are given in Table 6. Also listed are the results ob- 
tained with the 1933 flood eliminated. Both sets of data suggest that the 1933 
observed flood peak of 11,200 cfs is an unusual event with a return period 
probably much greater than 100 years. 


Synthesized Annual Floods for the Other 35 Project Basins 


Synthesized annual floods for the other 35 project basins (Table 1) provid- 
ed comparable degrees of agreement with the observed for the respective 
basins. Altogether, only 13 precipitation networks and 17 rainfall-runoff 
relations were used in synthesizing annual flood records for all 36 basins, 


CONC LUSIONS 


This study substantiates the conclusion that reasonably reliable flood- 
frequency data may be synthesized from rainfall records. Moreover, the 
precipitation stations need not be in the basin for which the synthesized flood 
peaks are required. The fact that the rainfall record from one network and 
one rainfall-runoff relation can apparently be used to synthesize flood- 
frequency data for a number of basins makes the procedure economically 
practicable, 

One point that must not be overlooked, however, is that there were no ade- 
quate precipitation networks in the project basins. This meant that there 
were no entirely suitable controls for determining whether the synthesizing 
procedure was yielding reliable values of peak discharge storm by storm. It 
is obvious that tests involving even but a few years of streamflow data with 
concurrent rainfall observations from an adequate network in the basin prope 
are to be preferred to the tests employed in these investigations. 

In spite of this shortcoming, the results of the studies made thus far 
strongly suggest that the synthesizing of peak discharges from rainfall is 
practicable and economically feasible. 
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Table 3 


FREQUENCY ANALYSIS OF OBSERVED AND SYNTHESIZED ANNUAL FLOODS 
FOR 1933-1954 


Return period Observed Synthesized 

(years) (cfs) (cfs) 

ce 2530 2730 

5 4850 5040 

10 6380 6570 

25 8320 8500 

30 9750 9940 

100 11180 11360 

Table 4 


FREQUENCY ANALYSIS OF OBSERVED AND SYNTHESIZED ANNUAL FLOODS 


FOR 1934-54 
Return period Observed Synthesized 
(years) (cfs) (cfs) 
2 2280 2470 
5 3610 3730 
10 4500 4.560 
25 5610 5610 
50 6440 6390 
100 7260 7160 
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Table 5 


SYNTHESIZED ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA, 
FOR 1900-1954 


CHRONOLOGICAL LISTING : STATISTICAL ARRAY 
Water Annual Order Annual Water 
Year Flood No, Flood Year 
1954 2890 1 11600 1933 
1953 2240 2 5420 1908 
1952 2340 3 5370 1942 
1951 4900 k 4900 1951 
1950 2590 5 4800 1946 
1949 1640 6 4270 1902 
1948 1930 T 4100 1913 
1947 1320 8 3910 1944 
1946 4800 : 9 3670 1934 
1945 2810 : 10 3420 1916 
1944 3910 : 11 3040 1915 
1943 1750 : 12 2950 1927 
1942 5370 : 13 2890 1954 
1941 1340 : 14 2850 1903 
1940 1970 : 15 2810 1945 
1939 1940 : 16 2730 1938 
1938 2730 : 17 2660 1926 
1937 2360 : 18 2590 1950 
1936 1780 : 19 2360 1937 
1935 1330 : 20 2340 1952 
1934 3670 : 21 2240 1953 
1933 11600 : 22 2210 1906 
1932 2200 : 23 2200 1932 
1931 1390 : 2k 2200 1924 
1930 2030 : 25 2140 1919 
1929 1380 : 26 2120 1907 
1928 1530 : 27 2030 1930 
1927 2950 : 28 1970 1940 
1926 2660 : 29 1960 1900 
1925 1380 : 30 1950 1917 
1924 2200 : 31 1940 1939 
1923 1200 : 32 1930 1948 
1922 1400 : 33 1820 1910 
1921 1220 : 34 1810 1905 
1920 1260 3 35 1780 1936 


(Continued next page) 
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SYNTHESIZED ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA. 
FOR 1900-1954 


‘CHRONOLOGICAL LISTING 


STATISTICAL ARRAY 


Water Annual : Order Annual Water 
Year Flood No. Flood Year 
1919 2140 36 1750 1943 
1918 1520 37 1730 1911 
1917 1950 38 1640 1949 
1916 3420 39 1600 1904 
1915 3040 ko 1590 1914 
1914 1590 KL 1530 1928 
1913 4100 ko 1520 1918 
1912 1250 43 1470 1901 
1911 1730 4h 1400 1922 
1910 1820 45 1390 1931 
1909 1270 46 1380 1929 
1908 5420 47 1380 1925 
1907 2120 48 1340 1941 
1906 2210 4g 1330 1935 
1905 1810 50 1320 194-7 
1904 1600 51 1270 1909 
1903 2850 52 1260 1920 
1902 4270 53 1250 1912 
1901 1470 54 1220 1921 
1900 1960 oy 1200 1923 
Table 6 


FREQUENCY ANALYSES OF SYNTHESIZED ANNUAL FLOODS 


ee 


Return period 1900-54 1900-32 
hoceen) 1934-54 
Pe ge ee ee ee ee 

2 2240 2160 

5 3840 3220 

10 4890 3920 

25 6230 4810 

50 7220 5470 

100 8200 6130 


nea EEE EEEEESSSE SEED 
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(Based on streamflow measurements 
near Riverdale and Hyattsville, Md.) 
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DISCHARGE (1,000 cfs) 


TIME (Hours) 


SIX-HOUR UNIT HYDROGRAPH 
CODORUS CREEK AT SPRING GROVE, PA., (74.3 Sq. Miles) 
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OBSERVED PEAK DISCHARGE (1,000 cfs) 


COMPUTED VS OBSERVED PEAK DISCHARGE 
CODORUS CREEK AT SPRING GROVE, PA., (74.3 Sq. Miles) 
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FLOW CHARACTERISTICS ON THE OGEE SPILLWAY 


Robert B. Jansen,* A.M. ASCE 
(Proc. Paper 1452) 


SYNOPSIS 


Although the ogee spillway has been the subject of more research than per- 
haps any other hydraulic structure, there are still deficiencies in our 
knowledge of its functioning. There is more to be learned about the process 
of energy loss on the spillway face and about the effects of submergence of 
the spillway jet. A general formula is here suggested for calculation of the 
spillway energy loss, and a theoretical derivation is presented defining the 
limit between the true hydraulic jump and the drowned jump. A criterion for 
estimating bottom velocities along a level spillway apron is considered. 


INTRODUCTION 


The determination of energy losses in flow on the face of an ogee spillway 
and the solution of the energy dissipation problem on the spillway apron often 
have been approached by means of hydraulic models. The reliability of such 
tests is universally accepted. The hydraulic laboratory plays a most im- 
s0rtant role in the design of these structures. The wealth of data on spillway 
energy loss produced by the many laboratories offers a good source of ma- 
ferial for use in the development of general mathematical criteria. This 
Japer attempts to bring this potentiality into focus. 

Experience gained thus far in design of the overflow spillway indicates 
hat this hydraulic structure has certain elements which lend themselves to 
standardization. In shaping the ogee crest, for example, it has become wide- 
y accepted practice to conform to the lower nappe of a stream issuing over 
1 sharp-edged weir. Uniformity has been obtainable in this case because of 


ote: Discussion open until May 1, 1958. Paper 1452 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 


_ Vol. 83, No. HY 6, December, 1957. 
‘Supervising Hydr. Engr., Dept. of Water Resources, State of California, 
Sacramento, Calif. 
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the easily predictable nature of the phenomenon of overflow. Probably due to 
the relative difficulties of measurement, the prediction of flow characteristic 
on the spillway face and apron has not been simplified to this extent. How- 
ever, by correlation of the results of experimentation by our hydraulic labor 
tories, such simplification should be possible. 


Energy Loss on Spillway Face 


A simple and exact criterion for determining the energy loss which will 
occur on the face of an overflow spillway would be a valuable design tool. 
Various mathematical procedures are now employed for this purpose. The 
application of the Manning formula has many advocates and is believed by 
them to yield satisfactory results. Solution by this means, though, is not di- 
rect and there are differences of opinion as to how the formula can best be 
applied. For example, in using the Manning equation, a common practice has 
been to determine the friction loss that would occur for flow at the depth with 
out air entrainment, but to allow for an assumed lower actual loss by reducin 
the value of the Manning coefficient of roughness. Yet, in the words of one 
reviewer of this paper, it is not even universally agreed that entrainment of 
air in the water has the effect of reducing the friction loss. 

Flow characteristics on the apron below the dam will depend, among other 
factors, upon the energy that remains in the jet when it enters the tail water, 
The height of the spillway, the head on the crest, the shape of the crest, and 
the position of the tail water will all affect the amount of energy thus involvec 
The theoretical energy content of the jet will obviously vary directly with the 
vertical distance of travel. The amount of water involved, and therefore the 
total energy rate, will depend primarily upon the crest head and the crest 
form, 

It is reasonable to assume that the energy loss on the spillway face is de- 
termined by a simple natural law. The factors involved are few, and their 
influences are relatively uncomplicated. The principal components of the 
formula would be the rate of discharge and the height of fall. The roughness 
and slope of the spillway face also affect the energy loss, and would have to 
be taken into account. The influence of these factors could probably be in- 
corporated into a coefficient in the formula. 

In the Madden Dam spillway prototype tests, (1) Pitot-tube measurements 
of jet velocities were made for a range of crest heads. The loss in energy 
head on the spillway face was taken as the difference in elevation between 
the reservoir water surface and the energy grade line calculated on the basis 
of the Pitot-tube readings. In order to appraise the accuracy of these ve- 
locity tests, the method of construction of the Pitot tubes should be noted. 

As described by Mr. Randolph, each installation consisted of a plate steel 
pedestal embedded in the concrete of the spillway apron, and projecting from 
14 to 20 inches above the floor. The Pitot tube itself was a length of 1-inch 
pipe protruding 4 inches from the nose of the pedestal at a height between 3 
inches and 18 inches from the floor. Although readings on these Pitot tubes 
were thoroughly and carefully taken, it can not be assured that the resulting 
accuracy is comparable with that obtained with laboratory instruments. Als¢ 
since the tubes were located close to the floor, they gave indications of the © 
maximum velocity rather than the mean. Despite these apparent limitations, 
the writer considers the Madden Dam velocity tests to be reliable enough to 
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serve in a preliminary general analysis of basic relationships. 

The head loss between the spillway crest and any point on the spillway 
face, for a given concrete roughness, should be a function of the height of 
fall and the rate of discharge over the crest. It is suggested, therefore, that 
the free-flow loss on the spillway face may be defined by the general ex- 
pression: 


= (P_+ H}@ 
Jeet . (1) 


wherein K is a coefficient, P is the height of crest above the downstream 
apron, H is the total head on the crest, q is the discharge per unit length of 
crest, and “a” and “b” are constant exponents. Fig. 1 shows the relationship 
of the basic dimensions. 

Application of this general formula to the prototype observations on the 
spillway of the Madden Dam shows encouraging possibilities. Using values 
of K = 1/8, a = 3/2, and b = 1/3, an acceptably close correlation can be ob- 
tained. However, the face of the Madden spillway has a roughness somewhat 
greater than average. For a typical concrete spillway of common roughness 
and comparable in other respects to Madden Dan, it is believed that the 
following equation will give a reasonably close measure of the head loss. 


_ 0.10(P + H)9/2 
Ee qi/3 (2) 


Substitution of CH3/2 for q in Equation (2) leads to 


_ 0.10(P + H)9/” 


5 cl/S,1/2 (3) 


in which C is the coefficient of discharge over the crest at the given head. 
Mr. Randolph in his paper attempted to derive an expression for the 
friction loss on the ogee spillway face. The following two formulas, convert- 
ed to the notation used herein, were developed by him for the Madden Dam 
spillway: 


= a 

a ad ead (P + H} (4) 
and 
4 
‘ 
_ 0,00934 P(P + HJ 
: Mae 0.623 (5) 
: H 
; 


Equation (4) is dimensionally correct but Equation (5) is not. The deficien- 
cies of these expressions, either individually or comparatively, are demon- 
strated by the substitution of the Madden Dam spillway height of 172 feet and 
a crest head of one foot. The result in each case is an energy loss greater 
than the actual head available. It is believed that Equation (3) will be found 

to be a more generally applicable conclusion. The implication is not intend- 
ed, however, that Equation (3) itself is a final refinement. It is merely an 
approximation independently drawn from the basic data presented by Mr. 
Randolph in his highly commendable work. 

_ Mr. Randolph, in his paper on the Madden Dam tests, stated that he was 
suggesting the installation of Pitot tubes on the aprons of other dams designed 


cael 
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by the U. S. Bureau of Reclamation. If data from such Pitot-tube measure- 
ments are available for a full range of heights of prototype dams, a formula 
probably of a form similar to Equation (1) can likely be derived which will 
remove the uncertainties from this design problem. The use of data obtained 
from high prototype spillways will enable the extension of the criteria to in- 
clude the effect of air insufflation. 


Type of Flow on Spillway Apron 


In order to provide full assurance against structural failure at the toe of 
an overflow dam, the desirable flow characteristics on the spillway apron 
must be known, along with the means of maintaining these characteristics. 
Formation of a true hydraulic jump may be desired, but not always feasible 
to attain. When the tailwater depth exceeds a certain limit, the jump can be 
expected to become submerged and the efficiency of energy dissipation ac- 
cordingly reduced. Of importance to the designer, therefore, is the determin- 
ation of the limiting condition at which the jump will lose its true character- 
istics. 

The criterion of demarcation between the true hydraulic jump and the 
drowned jump can be established theoretically. This is demonstrated in the 
analysis which follows. 

Using the term H to denote the total energy head on the crest (including 
velocity head) and D, as the depth of flow at the base of the overfall, the ve- 
locity at the base of the overfall is 


Vi = 2e(P +H - Dy - hy ) (6) 


wherein hy, is the total head loss from the crest to the apron, and g is the ac- 
celeration of gravity. 
The depth of flow at the base of the fall can be expressed as 


2 cH/2 


1 /ee(P + H- D, - hy) (7) 


The hydraulic jump will become drowned at a tailwater depth greater than 
the sequent depth Dg which would occur at the downstream end of a true jump. 
The limit between the true jump and the submerged jump can therefore be 
defined by this sequent depth. The theory of conservation of momentum gives 


2 

Do = BL ae sl aa 
2 3 

gD, ‘ 

(8) 


By substitution of terms derived above, the following dimensionless ratio 
is obtained: 


¥ 
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eongensg S| RSE ails 
(9) 


By plotting the ratios H/P versus D9/P, obtained through application of 
Equations (7) and (9), a family of limiting curves for various values of the 
coefficient of discharge can be developed. Fig. 2 serves as an illustration. 
For convenience in determining the depth at the end of the overfall, Fig. 3 
may be used instead of Equation (7). 

In planning for the design of the Boulder Canyon Project, the U. S. Bureau 
of Reclamation conducted an extensive hydraulic model testing program, much 
of which was devoted to a study of performance of overflow spillways. This 
program has provided experimental data relating to the limit between the 
true jump and the drowned jump. 2 Fig. 4 was derived from these obser- 
vations by the Bureau of Reclamation. Comparison of Figs. 2 and 4 indicates 
that the limiting depth observed in the laboratory was somewhat higher than 
that obtained by the theoretical formulas. This difference is noted even when 
the loss on the spillway face is neglected in applying Equations (7) and (9). It 
is suggested that perhaps Fig. 2 defines the limit at which the jump first tends 
to move upstream, and that the limiting curve of Fig. 4 marks the tailwater 
depth at which the drowning process is complete. 


Bottom Velocities on Spillway Apron 


Of primary concern to the designer is the distance traveled by the high- 
velocity stream under the tail water. On this will depend the extent of the 
works necessary for protection of the streambed against erosion. It is im- 
portant that a criterion be established to locate the point at which the velocity 
of the jet has been reduced to a safe magnitude. 

In the testing program for the spillway of Madden Dam, a set of measure- 
ments was taken of bottom velocities along a horizontal apron below a model 
overflow dam. An analysis has been made of the data reported from these 


tests. The relations drawn in Fig. 5 show the most promise in interpreting 


the experimental results. The graph is dependent upon limited data and there- 


fore should not be applied to general use. It is presented here to illustrate 


the possibilities of such analysis. The curves in Fig. 5 apply to a level apron 
without baffle piers or deflector sill. The variables used in the graph are 
defined as follows: The term “x” is the horizontal distance from the inter- 
section of the downstream spillway slope with the level apron (neglecting the 
curve of the bucket); Vg is the velocity of flow at the point where the stream 
on the spillway face passes the elevation of the tail water at the downstream 
end of the apron; d, is the depth of flow, measured perpendicular to the spill- 
way face, at the same point as Vo; hyo is the velocity head corresponding to 
V,; Dis the depth of tail water at the downstream end of the apron; and V is 
the bottom velocity of the jet on the apron. 


CONC LUSION 


This presentation has been exploratory in nature. The primary factors 
affecting the characteristics of spillway flow have been considered, along with 
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suggested means of relating them. Basic relations developed here are de- 
pendent upon a limited experimental foundation and are therefore intended to 
serve only for comparison with other model or prototype observations. It is 
expected that the discussions which are invited will produce further records 
from observed spillway flows that may enable the formulation of more wide- 
ly usable criteria. 

Particular emphasis has been placed herein on the determination of ve- 
locities of flow on the spillway face and apron, and on the phenomena influ- 
enced by these velocities. It is believed that a simple natural law governs 
flow on the face of an ogee spillway, although the exact relationship has not 
yet been written. The formula for energy loss here offered for comparison 
may help in the writing of a dependable rule that can be applied generally to 
overflow spillways of any shape or size. 

It is hoped that this paper will provoke discussions which will contribute 
knowledge not yet available to the profession at large. If the skeleton of cri- 
teria appearing here is given sustenance drawn from the experience of others, 
the purpose of this work will be served. 
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DISCHARGE CHARACTERISTICS OF RECTANGULAR 
THIN-PLATE WEIRS 


Carl E, Kindsvater,! M. ASCE, and Rolland W. Carter,2 A.M. ASCE 
(Proc. Paper 1453) 


SYNOPSIS 


A comprehensive solution is proposed for the discharge characteristics of 
ectangular, thin-plate weirs. The solution is based on a simple equation of 
ischarge and experimentally derived coefficients which account for the influ- 
nce of the fluid properties and the physical characteristics of the weir and 
1e weir channel. The effects of viscosity and surface tension are related to 
n increase in the effective head and an increase or decrease in the effective 
otch width. Thus, the combined effects of the fluid properties are accounted 
r with adjustment coefficients which are applied to measured values of the 
ead and width. Consequently, the coefficient of discharge is defined as a 
inction of the width-contraction ratio and the head-weir-height ratio only. 

Values of the coefficients required to describe the discharge of water over 
ectangular, thin-plate weirs are determined from original tests covering a 
ide range of the significant geometric ratios. Additional data needed to con- 
rm the proposed method of analysis are obtained from the published results 
f experiments by other investigators. The proposed discharge equation is 
Impared with some of the most widely used formulas. 


INTRODUCTION 


The rectangular, thin-plate weir is defined in this report as a basic weir 
rm which includes rectangular-notch weirs as well as full-width or 


te: Discussion open until May 1, 1958. Paper 1453 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 6, December, 1957. 
Regents Prof., School of Civ. Eng., Georgia Inst. of Technology, Atlanta, 
Ga., and Consultant, Geological Survey, U. S. Dept. of Interior, Washington, 
De C. 
Chief, Research Section, Surface Water Branch, Water Resources hee 
Geological Survey, U. S. Dept. of Interior, Washington, D. c 


1453-1 


1453-2 HY 6 December, 194 
! 


“suppressed” weirs. A typical rectangular-notch weir in a rectangular 
channel is shown in Fig. 1. 


Fic. l— THe Basic THIN-PLATE NotcH WEIR 


For many years hydraulicians have sought a comprehensive equation thé 
would describe the discharge characteristics of rectangular weirs for a ful 
range of fluid, flow, and geometric variables. Because weir discharge in- 
volves free, curvilinear flow and the combined influences of several fluid 
properties, it is not subject to complete mathematical description. Most i 
vestigators have resolved their difficulties by restricting their research tc 
the flow of water over full-width weirs of considerable height and width. lL 
fortunately, however, different workers continue to disagree regarding the 
physical limitations as well as the discharge characteristics of even this s 
called standard measuring weir. Thus, there still exists a variety of weir 
formulas, each based on a different set of experimental data, which disagr 
so much that the weir is questioned as an accurate measuring instrument. 

Most published weir formulas are related to a quasi-rational analysis i 
which the weir is described as a limiting example of the two-dimensional 
orifice. The traditional form of the basic equation of discharge is obtaine 
by integrating an approximate velocity equation over an approximate area. 
The result is a complex equation which must be modified to make it practi 
able and adjusted to make it agree with a particular set of experimental d: 
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he fact that a large number of different empirical formulas have been evolved 
this procedure is evidence that the “theoretical” equation is inadequate, 

he obvious reason for the inadequacy is the fact that the basic assumptions 

re inaccurate, 

Because the flow pattern for the weir is not yei subject to mathematical de- 
eription, the most direct solution for the discharge function involves a com- 
ination of dimensional analysis and experiment. Such an approach to the 
roblem is used in this report. 


Analysis of the Discharge Function 


lescription of the Weir 

Fig. 1 shows a typical example of the basic weir form which is the subject 
f this investigation—a symmetrical, sharp-edge, horizontal-crest, rectangular 
otch in a smooth, vertical, thin plate located in a smooth, horizontal, rec- 
angular channel. The nappe is fully ventilated and unsubmerged. The channel 
pstream from the weir, described hereafter as the approach channel, is as- 
umed to be of sufficient length to develop the normal velocity distribution for 
li discharges. This ideal condition is seldom realized in practice. In fact, 
1e use of short channels with various arrangements of baffles to control the 
elocity distribution is a major cause of disagreements in the results of ex- 
eriments made by different investigators. 


imensional Analysis 

The geometry of the flow pattern shown in Fig. 1 is described by the width 
f the notch (b), the width of the approach channel (B), the height of the weir 
rest above the bottom of the channel (P), and the piezometric head (h) re- 
srred to the level of the crest and measured in the uniform flow section up- 
fream from the weir. The fluid properties involved are the specific weight 
F), the density (?), the viscosity (4), and the surface tension (¢). Only one 
idependent flow characteristic is involved, and it can be represented by either 
le discharge (Q) or the head (h). However, h is already involved as an inde- 
sndent geometric variable, and Q is conveniently selected as the dependent 
iriable. Thus, a complete statement of the discharge function will include 
oth h and Q, as in the expression 


| On £,(b, B, P, h,T, Ps Mo) (1) 
From Eq. 1, dimensionless ratios which describe the discharge function 
in be formed as follows: 


Deans 


BEOie ae 
bh Vgh 


i which g = 7//, the acceleration due to gravity. The dependent ratio in Eq. 
is a coefficient of discharge. The first three independent ratios describe 

e geometry of the flow boundaries, and the last two symbols represent the 
2ynolds and Weber numbers. 

‘In American engineering practice the acceleration due to gravity is includ- 
|in the definition of the coefficient of discharge. Thus, for practical pur- 
es, a convenient definition of the coefficient (C) is 
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2 by beh ( 
= = £3 (2) 2? BPR, Ww) 
pho/2 By hie? 
from which 
Q =c bho/2 ( 


The dimensions of C in Eqs. 3 and 4 are the dimensions of Vg. Because 
its obvious simplicity, and despite its lack of dimensional purity, Eq. 4is us 
as the basic discharge equation in this report. 

The discharge function represented by Eq. 3 cannot be evaluated by snail 
cal procedures. Accordingly, the relative influence of each of the independe: 
ratios must be evaluated by experiment. 


Significance of the Geometric Ratios in Eq. 3. 


One of the most significant ratios in Eq. 3 is b/B. It is a measure of the 
channel-width-contraction characteristic of the weir. In combination with h, 
it is also an area-contraction ratio. The influence of b/B is similar to that 
of the corresponding width or diameter ratios which are used to describe the 
geometry of orifices. A better analogy is the width ratio which is used to de 
scribe open-channel constrictions. 1) Strangely, in view of the great volume 
of recorded research on weirs, the relative influence of b/B has received lit 
attention. The weir which has been investigated most extensively in the labc 
ratory is the full-width weir, for which b/B = 1.0. 

The b/h ratio can be described as a shape parameter. The independent 
influence of this ratio is believed to be negligible over the full practical ran; 
of the other variables. An earlier investigation 2) at the Georgia Institute 0: 
Technology supports this conclusion. The fact that most of the published re: 
sults of research on full-width weirs ignore the b/h ratio indicates either th 
it was not recognized as an independent ratio or that its influence was not e\ 
dent from the experimental data. A few recorded efforts to incorporate the 
b/h ratio in discharge formulas are believed to be based on misinterpretatic 
of influences related separately to the magnitudes of b and h. 

The h/P ratio is a measure of the depth-contraction characteristic of the 
weir. It is complementary to b/B as an area-contraction ratio. Its influenc 
is said to be represented by a velocity-of-approach term in several publish 
formulas. It has also been described as a Froude number for the flow in th 
approach channel. Actually, for the boundary conditions specified, the flow 
the approach channel is always tranquil, and flow with minimum specific en 
gy always occurs in the vicinity of the weir crest. In other words, the Frou 
number is not a significant independent ratio. (The Froude number would b 
significant if, for example, the approach channel were sloped or if the weir 
were located at the foot of a spillway or sluice.) Thus, to describe the h/P 
ratio as anything but a geometric parameter is to disguise its fundamental 
significance in the discharge function. 


Significance of R and W in Eq. 3 : 4 
Very little is known about the separate influences of viscosity and surfac 
tension which are represented by the Reynolds and Weber numbers in Eq. 3 
It is known, however, that the relative influence of the two fluid properties ; 
creases as the head on the weir and the size of the weir decrease. 
The Reynolds number (R) is a measure of the relative influence of visco: 
ty. It is usually expressed as 
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VL 

R 2 es = ae (5) 
n which V is a typical velocity, L is a conveniently evaluated and physically 
significant length, and v is the kinematic viscosity of'the fluid. For large, 
ull-width weirs the most significant length is the head. For small, narrow 
veirs, however, the width (and possibly the height) of the weir, as well as the 
ead, are independently significant. In general, therefore, because the Rey- 
olds number must be related to both the head and the width of the weir, there 
ire at least two forms of R. Furthermore, because velocities in the vicinity 
uf the weir crest are proportional to the square root of the head, Vh can be 
substituted for V in Eq. 5, whence Vh(h)/v and Vh(b)/v are independent ratios 
vhich have the significance of Reynolds numbers. Thus, for a given fluid (V 
onstant), the total influence of viscosity ona given weir can be expressed in 
erms of the absolute magnitudes of h and b, or h alone. 

The Weber number (W) is a measure of the relative influence of surface 
ension. It is usually expressed as 


(6) 


n which the surface tension (¢) and the density () are essentially constant 
or a given liquid. Because V is proportional to Vh, and because either h or b, 
r both, can be significant lengths (L), the ratios Vh(Vh)/Vo/p and Vh(Vb)/ 
/o/ are independent ratios which have the significance of Weber numbers. 
‘hus, there are at least two independent forms of the Weber number as well 
S the Reynolds number. Furthermore, the effects of both viscosity and 
urface tension for a given weir and liquid are related to the absolute magni- 
udes of hand b. Therefore, it is impossible to distinguish the separate ef- 
ects of the two fluid properties from experiments with a single liquid. 
Because the Reynolds number for any flow pattern is inversely proportional 
9 a typical viscous-shear force, the relative influence of viscosity decreases 
s Rincreases. Similarly, because the Weber number is inversely proportion- 
1 to a typical surface-tension force, the relative influence of surface tension 
ecreases as W increases. It follows from the foregoing definitions of R and 
V in terms of h and b that the relative influence of the combined fluid proper- 
les diminishes as either h or b becomes larger. Thus, for large heads on 
irge weirs, the influence of viscosity and surface tension is negligible. This 
onclusion is substantiated, for example, by the observation that the coefficient 
f discharge for large weirs is essentially constant for all heads above a 
ertain minimum. In other words, C is independent of the fluid-property para- 
1eters for large values of h and b. Here the distinction between “small” and 
large” must be based on systematic experimental investigations of the inde- 
endent relationships between C and h and between C and b. Many investigators 
ave attempted to determine the relationship between the absolute head and 
e discharge coefficient for full-width weirs. The writers know of no system- 
tic effort to evaluate the independent effects related to the absolute width of 
le weir. 
If attention is henceforth restricted to water at ordinary temperatures, it 
lows from the preceding discussion that the quantities h and b can be used 
ffectively to replace both R and W in Eq. 3. Furthermore, if the b/h ratio 
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is discarded on the basis of experimental evidence that it is insignificant, the 
discharge function represented by Eq. 3 can be changed to 


c= £,@ 2 h, b) (7 


which must be evaluated by experiment. 


Evaluation of Eq. 7 

The function represented by Eq. 7 involves four independent variables, a 
condition which precludes the presentation of experimental data in simple 
graphical form. However, if the relatively minor effects represented by h 
and b are described by an approximate, independent, empirical equation, the 
coefficient of discharge (C) can be expressed as a function of b/B and h/P 
alone. In other words, if the parameters representing the influence of viscos: 
ty and surface tension are effectively “removed” from Eq. 7, the relationship 
between the coefficient of discharge and the principal geometric ratios can be 
shown by a family of plane curves. 

The effects of viscosity and surface tension cannot be described by exact 
physical equations. Nevertheless, any procedure whereby these effects are 
eliminated from the discharge function must be based on a general understanc 
ing of the manner in which the fluid properties influence the flow pattern. 


The Influence of Viscosity 

A flow pattern which is determined by boundary conditions alone is de- 
scribed as a potential-motion pattern. The influence of viscosity in a real 
fluid motion can be illustrated by comparison with its potential-motion counte 
part. 

One of the effects of viscosity which can be ignored in many accelerated 
fluid motions is the energy loss which results from internal shear. The total 
loss of energy between the section of head measurement and the crest of a 
weir, for example, is generally insignificant. Therefore, the principal effect: 
of viscosity on weir flow are those which are associated with flow-pattern 
modifications due to boundary resistance and separation. 

A separation zone occupied by a large eddy occurs in the corner between 
the weir plate and the bottom of the approach channel, as shown in Fig. 2(a). 
The effect of this occurrence on a low weir, in comparison with the potential 
pattern for the same weir, is similar to the effect of sloping the weir plate 
downstream. That is, compared with its potential-motion counterpart, the 
vertical trajectory of the nappe is lower and the coefficient of discharge, bei 
virtually a coefficient of contraction, is larger. The comparable effect on thi 
discharge, therefore, is the same as an increase in head. 

A separation zone similar to that shown in Fig. 2(a) also occurs in the 
corner between the sides of a notch-weir plate and the walls of the approach 
channel. The effect of this occurrence, in comparison with potential motion, 
is to reduce the width contraction of the weir nappe. Its influence on the dis- 
charge, therefore, is the same as an increase in the width of the notch. 

Viscous shear causes the flow to be retarded in the vicinity of the bounda- 
ries. This condition can be described with respect to potential motion by 
means of the discharge-displacement, boundary-layer thickness, 6 * , shown 
schematically in Fig. 2. The cross-hatched areas in the figure represent the 
virtual displacement of the potential motion due to boundary resistance. On 
the surface of the weir plate, near the vertical edges of the notch as well as 
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Potential saree 
Real veiocity = 


(b) 


Fic.2—INFLUENCE OF VISCOSITY AND SURFACE TENSION (SCHEMATIC) 


he horizontal crest, the effect of the boundary layer is similar to the effect of 
Ounding the upstream corners of the notch. Compared with potential flow 
trough a sharp-edged notch, therefore, the coefficient of contraction for the 
eal fluid motion is larger. The effect on the discharge, on the other hand, is 
imilar to the effect of increasing both h and b. 

_ A boundary layer is also formed on the bottom and side walls of the ap- 
roach channel. The effect of this occurrence is particularly significant for 
irge values of h/P and b/B. For example, the flow through full-width weirs 
»/B = 1.0) is characterized by displacement boundary layers on both sides of 
1e channel at the weir section. The effect of this occurrence on the discharge 
; similar to the effect of a reduction in the width of the weir (b). For values 

f b/B less than 1.0 the effect of the side-wall boundary layers is to increase 
le effective value of b/B by decreasing the effective channel width (B). Simi- 
ily, the influence of the boundary layer on the bottom of the approach channel 
2sembles the influence of a decrease in P, The effect of this occurrence on 
ie discharge function, therefore, is similar to the effect of an increase in 

LP, 
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The influence of viscosity has been associated with the occurrence of 
boundary layers and separation zones. It should be emphasized, however, th 
the thickness of the boundary layers and the size of the separation zones are 
not directly proportional to the magnitude of the head or the size of the weir. 
(For this reason, the effects of viscosity are often described as “scale” ef- 
fects.) Thus, for example, considering the small absolute magnitude of 6 *, 
the relative effect of each of the several boundary layers depends on the abs 
lute magnitudes of h, b, and P. This observation is consistent with the previ 
ous conclusion that the relative effects of viscosity diminish as the size of th 
weir and the head on the weir increase. 


The Influence of Surface Tension 

Lindquist,(3) discussing the influence of surface tension on thin-plate, full 
width weirs, described two significant phenomena. He observed, first, that 
the nappe clings to the top edge of the crest as shown in Fig. 2. Thus, com- 
pared with its potential-motion counterpart, a real fluid flow is characterizec 
by a lower trajectory and a larger coefficient of contraction because of this 
surface-tension phenomenon. The effect of the occurrence is similar to the 
effect of rounding the upstream edge of the plate. Therefore, the effect on th 
discharge is the same as an increase in head. 

The second phenomenon described by Lindquist is the result of surface 
tension in both the upper and lower nappe surfaces. The effect of this oc- 
currence is a resultant force, F,, shown schematically in Fig. 2(b), which 
acts in the direction of the center of curvature of the nappe. Because the rac 
of curvature of the nappe surfaces decrease with decreasing head, and becau: 
the surface-tension force varies inversely with the radius of curvature, the 
resultant force Fy increases as the head decreases. Furthermore, because 
the resultant surface-tension force has a dominant downward component, it 
has the same effect on the discharge as an increase in head. 

The surface-tension phenomena described for full-width weirs are equally 
effective on notch weirs. In addition, similar phenomena occur at the vertic: 
sides of the notch as well as the crest. Thus, the total influence of surface 
tension on the discharge through rectangular-notch, thin-plate weirs can be 
compared with an increase in the width of the notch as well as an increase in 
the head on the crest. It should be emphasized, however, that the relative 
influence of surface tension depends on the absolute magnitudes of h and b. 
Therefore, like the influence of viscosity, the influence of surface tension di: 
minishes as the size of the weir and the head increase 


A Comprehensive Discharge Equation : 

By comparing the flow of a real fluid with the potential flow for the same 
boundary conditions, the influence of the fluid properties has been related to 
an increase in the effective head, an increase or decrease in the effective 
width of the weir, and a decrease in the effective height of the weir. It has 
been suggested also that the coefficient of discharge can be expressed in ter 
of the principal geometric variables alone if the effects of the fluid properti« 
are described by approximate, empirical equations. It is proposed, therefor 
that the effective values of h and b be defined as 


he 


Ml 


h + ky 


ya ae Bey ( 
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in Which the quantities kj, and k, represent the combined effects of the several 
phenomena attributed to viscosity and surface tension. It is also proposed 

hat a coefficient of discharge which is independent of the size of the weir or 
he magnitude of the head be defined as 


Ps QoS igi fey. 
ce be hg3/2 5G > (10) 


tis assumed that ky and kp as well as C., can be defined by experiment. Thus, 
rom Eq. 10, the equation 


ere (11) 


S proposed as a comprehensive equation of discharge for all rectangular, thin- 
late weirs. The adequacy of this equation, and, in fact, the adequacy of the 
ffective-head and effective-width concepts, must be examined in the light of 
xperimental data. 


Sources of Verification Data 


ixperiments of Bazin 

The classic experiments of Henri Bazin(4; 5) on full-width, thin-plate weirs 
vere made in 1886 and 1887 in a flume located on the bank of the Canal de 
sourgogne, near Dijon, France. Bazin’s setup included a sluice-controlled 
orebay and a concrete flume 6.56 ft wide and approximately 700 ft long. His 
feirs consisted of a crest plate mounted on a timber bulkhead. The nappes 
vere fully ventilated on the downstream side. The crest pieces were described 
S being iron, carefully straightened, 0.276 in. thick, sharp-edged, and with 
op surface perpendicular to the plane of the weir. Several different weirs 
vere actually used in making the tests, but they were described as being 
perfectly comparable with each other.” All of the weirs were of the full- 
fidth type (b/B = 1.0). In order to cover a full range of values of h/P, Bazin 
aade ten series of tests. The scope of his work is shown in Table 1. 
_ Series 1, 2, and 3 of the Bazin experiments were described as reference 
ests, because they alone involved volumetric measurements of discharge. 
‘or all other tests the weirs were calibrated by comparison with the reference 
reirs. Head measurements were made with hook gages in stilling wells con- 
ected to large orifices in the walls of the flume 16.4 ft upstream from the 
eirs. 
In general, Bazin’s instrumentation was good and his technique was meticu- 
Jus. His results differ considerably from those of most other investigators, 
Dwever, and this is usually attributed to the conclusion that, for different 
easons, the crest pieces were not effectively thin and sharp-edged. 6 


eriments of Schoder and Turner 

One of the most extensive investigations of full-width weirs reported in the 
terature is that performed under the direction of E. W. Schoder and K. B. 
urner at Cornell University from 1911 to 1920. 7) Nearly 2500 tests were 
ade in this investigation, but many of these involved velocity distributions 
hich were not uniform and crest plates which were not sharp-edged. The 

vo series of tests which have been selected for use in this report are believed 


| 
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TABLE 1. 


- SCOPE OF THE BAZIN ZXPERINENTS 


Series P in feet b in feet h in feet 


2 Bele 3.28 0.19 to 1.34 
3} 3.30 1.64 0.19 >to 1.76 
4 2.47 6.54 0.48 to 1.28 
5 2.47 6.55 0.49 to 1.44 
6 1.64 6.54 0.31 to 1.22 
7 1.64 6.55 0.34 to 1.40 
8 1.16 6.53 0.29 to 0.94 
9 1.14 Gs55 0.30 to 1.34 
10 0.79 6.55 0.30 to 1.34 


to have been made with “normal” velocity distributions. The scope of the 
tests used (Series E and F) is shown in Table 2, 


TABLE 2. - SCOPE OF THE SCHODER-TURNER SERIES E AND F TESTS (b/B = 1.0). 


Series P in feet b in feet h in feet 
E 0.50 4.22 0.029) to 1.41 
F 1.00 4.22 0.029 to 2,21 


The Cornell tests described in Table 2 were made in a wooden flume 4.2% 
ft.wide and approximately 30 ft.long. Because of the short length of the flur 
baffles were provided at the inlet to control the velocity distribution. A moy 
able bottom was used to obtain different weir heights. Discharges were 
measured volumetrically. The crest piece for this group of tests was made 
rolled brass, 1/4 in. thick. It was described as having a straight, square, a 
sharp upstream top corner, with a machined 3/32-in. wide top surface and a 
beveled downstream face. The crest piece was mounted on a timber bulkhe: 
Head measurements were made with a float gage connected to piezometers 


11.74 Ne nba from the weir. Complete details are given in the paper 
cited. 
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‘xperiments of the U. S. Bureau of Reclamation 

Experiments on full-width weirs were performed by the U. S. Bureau of 
teclamation in connection with hydraulic investigations for the Boulder Canyon 
‘roject. 8) A distinctive feature of the equipment used for the tests was a 
nife-edge, stainless-steel crest piece. The purpose of this innovation was 
) prevent the clinging-nappe phenomenon described previously in this report. 
‘he crest piece was mounted on a steel plate located in a sheet-metal-lined 
ooden flume 2 ft.wide and effectively 18 ft-long. A movable floor was used 
) obtain different values of weir height. Baffles located at the upstream end 
ere used to produce a uniform velocity distribution in the flume. Discharges 
ere measured with venturi meters in the supply line. Heads were measured 
fith a hook gage in a stilling well connected to a piezometer in the movable 
loor 9.9 ft.upstream from the weir. The scope of the USBR tests is shown in 
‘able 3. 


TABLE 3. - SCOPE OF THE USBR EXPERIMENTS (b/B = TsO). 


Series P in feet b in feet h in feet 

er eee ee 
1 5.00 2.007 0.14 to 1.26 
2 4.60 2.007 0.46 to 1.25 
3 3.50 2.007 0.14 to 1.22 
4 2.50 2.007 0.14 to 1.25 
3 2.00 2.007 0.24 to 1.19 
6 1.50 2.007 0.28 to 1.20 
7 1.01 2.007 0.25 to 1.18 
8 0,57 to 0.65 2.007 0.44 to 1.11 
9 0.32 to 0.41 2.007 0.40 to 1.04 

» 10 0.13 to 0.14 2.007 0.35 to 0.89 

Bil 0.08 to 0.09 2.007 0.34 to 0.71 


ew Data on Rectangular Weirs 

Original experiments on rectangular weirs were performed in 1955 in the 

draulics Laboratory of the Georgia Institute of Technology at Atlanta. For 
Georgia Tech tests the weirs were located at the downstream end of a 

eel flume 30 in. deep, 10 ft. wide, and 25 ft-long. The floor of the flume was 

vel and smooth. Its width was varied by means of movable walls faced with 

4-in, aluminum plates. For tests on full-width weirs the wall plates were 

ade to project through the notch and past the crest a distance of 6 in. Baffles 

the forebay and curved guide walls at the upstream ends of the movable 

were adjusted to produce a uniform velocity distribution. Surface floats 
re used at the entrance to damp waves in the flume. 
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The basic weir bulkhead was made of 3/8-in. aluminum plate. The width 
and height of the weir notch was varied by means of lap-joined plate sectio 
made of the same material. The notch edges, including the crest piece, wer 
made of 1/8-in. stainless-steel plates which were fastened to the aluminum 
plates with the upstream surfaces flush. The edge pieces were machined on 
a flat-bed planer to be smooth and flat on the upstream side and beveled on 
the downstream side. The top surface of the notch edges was 1/16-in. wide 
and perpendicular to the upstream face of the weir. The upstream corners 
the edge pieces were sharp. 

Water was supplied to the flume from a constant-head system. Discharge 
were measured in an electrically operated weighing tank at the upstream enc 
of the flume. Heads were measured with a hook gage in a stilling well which 
was connected to piezometers located in both movable walls 2 in. above the 
floor and 5 ft.upstream from the weir. 

The purpose of the Georgia Tech tests was to obtain the additional experi: 
mental data necessary to describe the discharge function over the full, pract 
cal range of the variables shown in Eq. 7. The scope of the 249 tests which 
were made at Georgia Tech is shown in Table 4. 


TABLE 4. - SCOPE OF THE GEORGIA TECH TESTS. 


P in feet b in feet b/B h in feet 
1.44 2.68 1.00 0.16 to 0.5 
1.44 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.14 to 0.5: 
0.56 2.68 1.00 0.10 to 0.6 
0.56 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.10 to 0.7 
0.30 2.68 1.00 0.17 to 0.6 
0.30 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.10 to 0.7 
0.30 1.20 0.13 0.11 to 0.6 
0.30 0.577 0.06 0.28 to 0.7 
0.30 0.281 0.03 0.27 to 0.6 
0.30 0.118 0.01 0.08 to 0.6 
0.30 0.400 1.00, 0.80, 0.50 0.15 to 0.7 
0.30 0.200 1.00, 0.80, 0.50 0.10 to 0.6 
0.30 0.100 1.00, 0.80, 0.50 0.10 to 0.6 
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Verification of the Proposed Solution 


ecapitulation 

The coefficient Cg in the proposed discharge equation (Eq. 11) is defined in 
aq. 10 in terms of the effective head, he, and the effective notch width, be. It 
s also defined as a function of the geometric ratios h/P and b/B. Unlike C 
Eqs. 4 and 7), therefore, Ce. is expected to be independent of the fluid-property 
‘fiects represented by the magnitudes of h and b in Eq. 7. 

It is assumed that C, can be determined from laboratory data if the 
neasured values of h and b are adjusted as indicated in Eqs. 8 and 9. In Eq. 
} the effective head is defined as the measured head plus a quantity kph. Simi- 
arly, in Eq. 9 the effective width is defined as the measured width plus a 
juantity kp). It is reasonable to conclude that the proposed discharge equation 
vill be practicable and convenient to use only if ky and kp are constants or 
unctions of not more than one independent variable. Thus, verification of the 
iroposed solution is critically related to the experimental evaluation of the k 
luantities. 


ivaluation of ky, 

The quantities kj, and kj, are evaluated from experimental data by a trial 
Omputation procedure. For a series of tests in which h is the principal vari- 
ble, kj, is determined by successive approximations as the quantity which will 
remove” the correlation between the coefficient of discharge and the head. 
‘he procedure is illustrated by its application to the results of the Schoder- 
‘urner tests. 

Fig. 3(a) shows C as a function of h/P for two series of tests made by 
choder and Turner on full-width weirs (Table 2). The scatter in values of C, 
hich is most apparent at lower values of h/P, is systematically related to 
le magnitude of h. If a small quantity (k,) is added to h, the value of the coef- 
icient is reduced by an amount which is inversely related to the value of h. 
‘alues of the adjusted coefficient (Ce) from a trial computation, using Ky = 
.003 ft, are shown in Fig. 3(b). A second trial, with k, = 0.004 ft, is shown 
1 Fig. 3(c). From Figs. 3(a) and 3(c) it is apparent that a constant value of 
h = 0.004 ft.effectively eliminates the correlation between the coefficient of 
ischarge and the head for Series E and F in the Schoder and Turner tests. 
rom Fig. 3(c), furthermore, Ce is conveniently related to h/P by a straight- 
ne equation. 

The Schoder and Turner Series E and F were selected for two reasons, 
irst, they included some tests with very low values of h. For these tests, 
lerefore, the effects of viscosity and surface tension which are related to the 
lagnitude of h are relatively large. Second, the width of the weir used for 
le tests was large enough (4.22 ft) that the influence of kp and Ce could be 
ssumed to be negligible. Proof of this conclusion is based on results of 
eorgia Tech tests which are described subsequently (See Fig. 12). 

Fig. 4(a) shows C as a function of h/P for the Bazin tests on full-width 
eirs (Table 1). Fig. 4(b) shows the same data in a plot of Ce as a function of 
/P, The value of kp used for computing C, for these tests was 0.012 ft. This 
alue was selected on the basis of trial computations as the value required to 
iminate the correlation (evident in Fig. 4(a)) between the coefficient of dis- 
large and the head. The effect of the width-adjustment factor, ky), was as- 
imed to be negligible for all the Bazin tests. This assumption is believed 
‘have had a small effect on the computed values of Ce for only one group of 
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tests, Series 3, for which b was 1.64 ft. From Fig. 4(b) it is apparent, again, 
that C, is a simple, linear function of h/P. 
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Fig. 5 shows C, as a function of h/P for the USBR tests on full-width weir 
(Table 3). For these tests, kj, was determined to be 0.003 ft. The scatter in 
the points for larger values of h/P is believed to be indicative of experiment: 
errors resulting from observational difficulties related to low values of P. 
Symbols are used in Fig. 5 to indicate different values of P. The width of the 
flume was approximately 2 ft. The effect of k, was assumed to be negligible. 

Only full-width weirs (b/B = 1.0) were involved in the tests made by Schoc 
and Turner, Bazin, and the USBR. Therefore, one of the principal purposes 
of the Georgia Tech tests (Table 4) was to determine the validity of the pro- 
posed solution as applied to notch weirs. Tests made specifically to evaluate 
kj, included a complete series for each of the following values of b/B: 0.20, 
0.40, 0.60, 0.80, 0.90, and 1.00. For the tests on notch weirs, the width (b) 
was 1.80 ft. For the tests on full-width weirs, b was 2.68 ft. Thus, for allo 
the kj tests the relative influence of effects related to b were negligible. Bo 
h and P were varied in each series of tests in order to cover a full range of 
values of h/P. Values of h varied from 0.10 to 0.76 ft, and values of P varie 
from 0.30 to 1.44 ft. Sufficient tests were made to define the coefficient of 
discharge from h/P = 0 to h/P = 2.4. Fig. 6 shows the results of a typical 
series (b/B = 0.80). 

The successive approximations procedure was used to evaluate ky for eac 
of the series of tests described above. Remarkably, the value of k, requirec 
to satisfy the effective-head criteria was the same for all values of b/B and. 
all values of h/P included in the Georgia Tech tests. In other words, for a 
full practical range of the significant variables a constant value of kh (0.003: 
is adequate to compensate for the fluid-property effects related to the head. 
This conclusion was subsequently verified by successfully applying the same 
value of kj, to additional tests covering values of b/B from 0.01 to 1.00 and 
values of b as small as 0.10 ft. Thus, the effective-head concept is substant 
ated as an important part of the proposed solution. 
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Differences in values of kj and in the Ce curves determined from Schoder 
and Turner, Bazin, USBR, and Georgia Tech tests on full-width weirs are be- 
lieved to have been caused by differences in test equipment. The larger value 
Mf ky, which characterizes the Bazin tests, for example, is believed to be an- 
ther indication that Bazin’s weir plate was not effectively thin and sharp. 
[his conclusion is supported by the reasoning that the effect of crest-rounding 
S similar to some of the effects of surface tension and viscosity. The results 
yf the USBR tests are doubtless influenced by the fact that the crest piece of 
he weir was knife-edged. The Schoder and Turner tests may have been influ- 
nced by the fact that their flume was very short and entrance conditions were 
oor for large values of P. Most of Bazin’s tests, on the other hand, were 
nade in flumes long enough to develop a normal, open-channel velocity distri- 
ution. Thus, the results of each investigation are characterized by the equip- 
nent used in making the experiments. 

_ A constant value of ky, was found to be adequate for all weirs tested in the 

reorgia Tech investigation. Different values of ky were determined for full- 

yidth weirs tested by other investigators. Therefore, although it is indicated 
hat the same value of kj, can be used for all values of b/B, the value of k;, to 

used is also dependent on the physical characteristics of the weir and weir 
hannel. 
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Evaluation of kp 
The weirs investigated by Bazin, Schoder and Turner, and the Bureau of 


Reclamation were too wide to reveal the fluid-property effects related to 
small values of b. Therefore, the results of these tests could not be used to 
evaluated k,. Some tests on weirs of relatively small width were included in 
the Georgia Tech experiments. It was the purpose of these tests to evaluate 
kp and to investigate the validity of the effective-width concept. 

In Fig. 6 the total displacement of values of Ce with respect to correspond 
ing values of C is a measure of the combined effects of kj, and k,. However, 
only a few points at low values of h/P represent heads so small that the adjus 
ment caused by k, was appreciable. The nearly uniform displacement at high 
er values of h/P, therefore, is the result of using a value of kp = 0.014 ft in 
the computation of C,. This value was determined from trial computations as 
the value of kh which would eliminate the correlation between b and C, in the 
results from a special series of tests with b/B = 0.80 (not shown in Fig. 6). 


Legend 
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Fig. 7 shows the results of special tests made to evaluate kp for another ~ 
value of b/B. Actually, several values of b/B are represented on the figure, 
but the results of the tests indicated that kp is essentially constant for values 
of b/B less than 0.2. In Fig. 7, as in Fig. 6, the effect of the kp adjustment i 
shown by the displacement of values of Ce with respect to corresponding | 
ue of C. The value of kp used in computing Ce for this example was 0.00 

Fig. 8 shows the results of all tests made to evaluate kp. The figure 
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Fic. 8.— VALUE OF kp DERIVED FROM GeorGIA TECH TESTS 


ndicates that k, is a function of b/B, that it increases with increasing values 
af b/B until it reaches a maximum positive value of 0.014 ft at b/B = 0.80, and 
hat it decreases rapidly thereafter, reaching a minimum value of -0.003 ft at 
/B = 1.0. The authors acknowledge some uncertainty regarding the generality 
Mf the curve in Fig. 8. It is believed that the value of kp obtained from the 
seorgia Tech tests may have been influenced, to an unknown extent, by the 
hysical characteristics of the equipment used. Tests made especially to 
valuate kp, involved very small values of b and correspondingly small values 
f Band Q. Thus, unavoidable errors of measurement may account for com- 
aratively large relative errors in the computed results. Nevertheless, the 
idequacy of the effective-width concept was demonstrated by the successful 
limination of b as an independent variable in the experimentally determined 
elationship between Cg, h/P, and b/B. 

_ A remarkable example of the experimental evidence which substantiates 
oth the effective-head and the effective-width concepts is shown in Fig. 9. In 
his figure a single, straight line satisfactorily represents the relationships 
etween Ce and h/P for full-width weirs involving values of b from 0.10 to 

-68 ft and values of h from 0.10 to 0.72 ft. 


valuation of Ce 

The results of the tests from all sources indicate that Ce, the “effective” 
efficient of discharge, is a function of b/B and h/P alone. Thus, in confir- 
aation of Eq. 10, Cg is independent of the influences of viscosity and surface 
ension. Furthermore, the results show that the relationship between C, and 
/P for all values of b/B is conveniently defined by a family of straight lines. 

Fig. 10 shows a summary of the results of the Georgia Tech tests. The 
iformation on this figure, in combination with values of kp and kp (Fig. 8) de- 
srmined from the same tests, is the information needed to evaluate the 


1453-20 HY 6 December, 195 


Legend 
Symbol P in feet b in feet 
0.30 2.68 
0.56 2.68 
1.44 2.68 
0.30 0.40 
0.30 0.20 
0.30 0.10 


kp = +0.003 ft 
kp = — 0.003 ft 


Value of Ce 


Value of h/P 


Fig. 9—Georcia TECH Tests (b/B = 1.0) 


ASCE 


value OTF Le 


4.2 


40 


| 3.0 


KINDSVATER - CARTER 1453-21 


Legend 


——- From Experimental Data 
—-— Interpolated 


N 
PX 


K 


A 


2 
Pe 


a 
mee 


Value of h/P 


Fic.10.—Ce For ALL VaLuEs oF b/B, FROM GeorRGIA TECH TESTS 


1453-22 HYG December, 195 


proposed discharge equation. Thus, Eq. 11, complemented by the effective- 
head and effective-width concepts, is confirmed as a comprehensive equation 
of discharge for all rectangular, thin-plate weirs. 


Summary 


Experimental Results ; 

The proposed discharge equation (Eq. 11) and the concepts of effective he 
and effective width (Eqs. 8 and 9) have been substantiated by experiments fro! 
various sources covering a wide range of test conditions. It has been demon- 
strated that C, is a linear function of h/P for each value of b/B, that ky is a 
constant for all values of b/B and h/P, and that kp is a function of b/B. 

From the results of tests made at the Georgia Institute of Technology, 
values of Cg, k, and kp have been determined for water at ordinary tempera- 
tures. Fig. 10 shows Ce for a full range of values of b/B, with h/P from 0 to 
2.0. From the information contained in this figure, simple equations for Cg 
can be written if desired. The constant value of ky determined from the 
Georgia Tech tests is 0.003 ft. Values of kp are shown in Fig. 8. 

Comparisons of values of Cg and kp obtained from tests made by different 
investigators of full-width weirs confirm the general conclusion that weir dis 
charge is critically influenced by the physical characteristics of the experi- 
mental setup. Especially critical are the weir-notch edges and the channel 
properties which control the velocity distribution and turbulence upstream 
from the weir. 

Neither superior equipment nor superior technique is claimed for the 
Georgia Tech tests. However, it is believed that the over-all quality of these 
data is as good as that of any other data available. Furthermore, this is the 
only investigation known to cover a full, practical range of the significant var 
ables. Therefore, as a reasonably accurate solution for the comprehensive 
discharge characteristics of rectangular weirs with effectively sharp, thin 
crests, the authors recommend the use of the proposed discharge equation an 
the coefficients derived from the Georgia Tech experiments. 


Effect of Neglecting kp and kp 

It has been established that the relative influence of viscosity and surface 
tension diminishes as both the head and the notch width becomes larger. In 
the proposed discharge equation the fluid-property effects are represented bj 
the head- and width-adjustment factors, ky and kp. Thus, the relative influ- 
ence of ky, decreases as the head increases in magnitude, and the relative 
influence of k, decreases as the width increases in magnitude. It follows tha 
for “large” heads on “large” weirs the influence of the adjustment factors is 
negligible. The experimental results can be used to evaluate the influence of 
kp and kp in relation to the magnitudes of h and b. 

A convenient approximation for the discharge can be obtained from Eq. 11 
using Cg from Fig. 10, with the unadjusted (measured) values of hand b. Th 
relative accuracy of this approximation is a significant measure of the influ- 
ence of ky, and kp. It is also a measure of the influence of the combined effec 
of viscosity and surface tension. Thus, if the absolute error (AQ) due to 
neglecting kp and kp is 


AO = C.b,h,°/2 = C,bh3/2 (1 
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jthe relative error is 


SL ch 3/2 
Q Bike (13) 


The influence of kj, alone can be demonstrated by assuming that the notch width 
ig very large, whence b = be in Eq. 13. Then, from Eqs. 8 and 13, 


ce Rai ya ade 
Q ae (14a) 


which is shown in Fig. 11 for kh = 0.003 ft. Similarly, the separate influence 
of k, can be demonstrated by assuming that h= he. Thus, from Eqs.-9-and 
13, 


Re eS 


Q b+ k, 


) 
(14) 


which is shown in Fig. 12 for values of ky from Fig. 8. Figs. 11 and 12 are 
recommended for use in connection with Eq. 11 as a means of determining 
when (depending on the accuracy desired) the influence of the fluid properties 
may be neglected. 


Upper Limit of h/P 

Boss has shown that critical depth will occur in the uniform flow up- 
stream from a weir with h/P greater than about 5. For values of h/P greater 
than 5, therefore, the weir is not a control, and the ordinary weir discharge 
equations are not applicable. A lower, practical limit on h/P is imposed by 
observational difficulties and measurement errors. For precise measure - 
ments the writers recommend that h/P be limited to values less than 2.0. 


Influence of Weir Height 

The weir height, P, was acknowledged as a possible third significant length 
9arameter in the discharge function. There is disagreement in the literature, 
lowever, regarding experimental evidence of a systematic correlation between 
P and the coefficient of discharge. Bazin 9) warned against the use of his dis- 
sharge Faia for a “very low weir. . . which should always be avoided.” 
Rehbock‘10) acknowledged “some disturbance, which is not yet clearly under- 
stood (and because of which) the coefficient C is not exactly the same for any 
wo sharp-crested weirs of different absolute heights, P, even if the ratios, 
1/P, are identical.” Weirs as low as 0.08 ft were used in the USBR tests in 
yrder to achieve the higher values of h/P. The results of these tests show a 
lefinite correlation between Ce and P for values of h/P greater than 2 (not in- 
luded on Fig. 5). The Georgia Tech tests, however, failed to show a similar 
correlation for values of P as low as 0.30 ft and values of h/P as high as 2.4. 

It is generally acknowledged that observational difficulties and errors in 
measurement may account for appreciable errors in the results of experiments 
nvolving small values of P and large values of h/P. For this reason alone, 
veirs lower than 0.3 ft are not recommended for precise measurement of dis- 
harge. Nevertheless, although questions regarding the influence of weir 
eights less than 0.3 ft may be academic, additional experiments are needed 
9 resolve the contradictions to be found in the literature. 
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Comparison with Other Experiments and Formulas 


: Some Classic Formulas 
: One equation is frequently described as the basic equation of discharge for 
all weirs. It is founded on an assumed analogy between the weir and the two- 
dimensional orifice. In its derivation an approximate velocity equation is inte- 
grated over the approximate area limits of the nappe in the plane of the weir. 
The effects of viscosity and surface tension are ignored. The result, in the 
form usually ascribed to Weisbach, is 
Q = C, 2V2¢ b [(h+ ae aR 
3 28 28 (15) 


in which C, is a coefficient of contraction and V1 is the average velocity in the 
approach channel. 

Some writers have erroneously implied that Eq. 15 is a comprehensive and 
theoretically correct equation of discharge. It is claimed, for example, that 
the velocity-head terms account for the influence of approach-channel geome- 
try. From this viewpoint h/P is eliminated as an independent variable, and 
C, is a function of b/B alone. These conclusions are not substantiated by ex- 
periment. 

The inadequacy of Eq. 15 can be attributed partly to the fact that the effects 
of viscosity and surface tension are ignored in its derivation. However, ex- 
periments show conclusively that these effects are appreciable only when h or 
b, or both, are relatively small. Eq. 15 is deficient largely because its deri- 
vation was based on several erroneous assumptions. 

Unjustifiable alterations of Eq. 15 have been made for the purpose of im- 
proving agreement between the equation and experimental data. Thus, it is 
quite common to omit the second velocity-head term in the brackets, and, 
subsequently, to attribute great significance to the “total head” term which 
remains. The result is an equation such as that proposed by Hamilton Smith, 
Jr., in 1884,(12) 

3/2 


2 
aN ers aS 
Q=C¢, 3 28 b [h + ky Ge )) (16) 


in which both Cg, and ky are coefficients determined from experiment. Smith 
described the coefficient kj as a number which represents “the portion of the 
kinetic energy [which] produces a useful effect at the opening of discharge.” 
He recognized that ky is a variable which depends on the velocity distribution 
in the approach channel. Nevertheless, he recommended a constant value of 
1.33 for all full-width weirs and 1.40 for all “fully contracted” weirs. He 
described C, as a discharge coefficient, and he concluded that it was a function 
of the head, the notch width, and the width-contraction ratio (b/B). Actually, 
both k, and Cg are functions also of the ratio h/P, and they are critically de- 
pendent on the physical characteristics of the weir and the approach channel, 
Therefore, Eq. 16 is little more than a basis for the experimental definition 
of the discharge function. It is obviously not the most convenient formula for 
this purpose, nor is it the comprehensive solution which is needed. 

Hamilton Smith’s exhaustive study and detailed exposition of the influence 
of the velocity of approach was subsequently ignored or misinterpreted by 
others. Some writers have inferred that his kj coefficient can be associated 
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with the kinetic-energy coefficient (#1). Accordingly, they have perpetuated 
the misconception that the coefficient of discharge is a constant when the 
quantity within the brackets (Eq. 16) is evaluated as the true total energy head 
in the approach channel [h + 1 (V42/2g)]. 

The coefficient @ is defined as the ratio of the true average velocity head 
in a cross section to the velocity head computed on the basis of the average 
velocity in the section. It is a dimensionless measure of the degree of velocity 
nonuniformity, but it does not account specifically for the shape of the velocity 
distribution curve. Smith, Schoder and Turner, ( ) and others who have studied 
the influence of the velocity distribution in the approach channel have shown 
conclusively that the pattern of nonuniformity is more critical than the degree 
of nonuniformity. Thus, experiments performed in a large variety of weir 
channels have failed to reveal a satisfactory correlation between a, and kj. 

Most of the defects of Eq. 16 are common also to a large number of 
equations in the form suggested by Weisbach in 1844, 


Q=¢, 3 V2elk, + 3 (2 =) Yon?/? (17) 
in which Cy, Ky and kg are experimentally determined coefficients which ac- 
count for the influence of the fluid properties, the velocity distribution, and 
the physical characteristics of the weir and approach channel. The principal 
advantage of Eq. 17 is that the velocity-head term in Eq. 16 is replaced by a 
term which involves the h/P ratio. The effect of this modification is to avoid 
a successive-approximations solution for discharge. Equations of this gener- 
al variety include some of the most widely used formulas for full-width weirs 
(b/B = 1.0). Among these are the Bazin formula, (5) 


3/2: 


Q = (3.25 + O07) U1 + 0.55(5 25) Ibh (18) 
the Rehbock formula of 1912,(13) 
Q = (3.24 + 0.43 Be 2.018) ph3/2 (19) 
and the S.I.A. (Swiss Society of Engineers and Architects) formula, (14) 
Q = (3.29 + FS Orggs)[L + 0.565257) vn? (20) 


The complex expressions which are the coefficients of the quantity bh3/2 in 
each of these equations are equivalent to C in the simple equation of discharge 
(Eq. 4). The numerical quantities are experimentally derived coefficients 
which account for the influence of the fluid properties as well as the weir and 
channel characteristics. Specifically, a term involving and having the di- 
mensions of h is now recognized as a means of compensating for the combined 
effects of viscosity and surface tension. Because it is a dimensional quantity, 
the h term restricts the application of Eqs. 18, 19, and 20 to the fluid used 
for the experiments (water at ordinary temperatures). However, it is apparent 
that As influence of the head term is negligible except when h is relatively | 
small, 


ee formula for full-width weirs was proposed by Rehbock in 
1929 5 
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Q = (3.22 + 0.44 3) bagi? (21) 


| This was the first formula known to involve the effective head (he = h + 0.004 
ft). The quantity 0.004 ft, corresponding to kj, in Eq. 8, was apparently derived 
by Rehbock from an analysis of his own experiments. 

Rehbock did not associate the h term in Eq. 19 or the effective-head term 
in Eq. 21 with the influence of the fluid properties. He described both terms 
aS a means of accounting for a disturbance related to the magnitude of P. He 
hesitated to endorse Prandtl’s suggestion that they were related to the effect 
‘of capillarity. (10) 

Various restrictions were placed on Eqs. 18, 19, 20, and 21 by their auth- 
ors. Bazin\9 specified that “the arrangement of our standard weir {must be] 
exactly reproduced.” Rehbock(15) required that the head be large enough to 
ensure a free nappe, but he claimed a high degree of accuracy for both of his 
formulas when applied to full-width weirs of all sizes. The S.I.A., in its 
Standards for Water Measurement, 14) restricts the application of Eq. 20 to 
values of P not less than 1 ft, values of h between 0.08 and 2.6 ft, and values 
of h/P not greater than 1.0. 

One of the new formulas developed for rectangular-notch weirs (b/B < 1.0) 
is that proposed by the S.I.A.: 14 


0.063 - 0.053(2)? 


Q = [3.09 + 0.20(2)2 4 
B h + 0.005 


b. 4 h 2. aya (22) 
[1 + 0.5C) C—) ] bh 


The application of this formula is restricted to values of P not less than 1 ft, 
values of h between 0.08(B/b) and 2.6 ft, values of h/P not greater than 1.0, 
and values of b/B not less than 0.3. Eq. 22 is widely quoted in European 
technical literature, but correspondence and library research failed to reveal 
the source of the data on which it was based. The formula is virtually unknown 
in this country. 

In American practice, approximate formulas for notch weirs are adapted 
from formulas for full-width weirs by means of the contracted-width equation 
ascribed to James B. Francis:(16, 17 

b. = b+ 0.2h (23) 


| 
4 


In this equation be is the width of the contracted stream, b is the gross width 
of the notch, and the quantity 0.2 is a coefficient recommended by Francis for 
“fully contracted” rectangular-notch weirs. Thus, from Eqs. 4 and 23, the 
Francis discharge formula for notch weirs is 


: . 
| Q = C, (b - 0.2 h) h?/2 (24) 


in which C is the coefficient of discharge in any of the formulas for full-width 


weirs (b/B = 1.0). 
Hamilton Smith, Jr., after reviewing the experiments of Francis, Fteley 
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and Stearns, and others, described the requirements for “full contraction” as 
follows:(17) The weir height (P) should be not less than twice the head (h) or 
not less than 1 ft; the width from the side of the approach channel to the end 
of the notch should be not less than twice the head or twice the width of the 
notch (b). Thus, Smith’s recommendations included the requirement that h/P 
be not greater than 0.5 and that b/B be not greater than 0.2. Francis 17) add- 
ed the requirement that b/h be not less than 3.0. It is recognized now that the 
concept of “full contraction” is inconsistent with the true significance of b/B 
as an independent variable. 


Comparison of Experiments and Formulas for Full-Width Weirs 

The experiments of Schoder and Turner, Bazin, and the Bureau of Recla- 
mation, as well as the Georgia Tech experiments, were used to confirm the 
method of analysis and the discharge equation proposed in this report. Conse- 
quently, values of Cg and ky, were determined, and these values can be used 
to compare the results of the tests made by the different investigators. Com- 
parable values of k, could not be determined because only the Georgia Tech 
experiments included tests on weirs narrow enough to reveal the fluid-property 
effects associated with small values of b. 

From the results of the Georgia Tech tests on full-width weirs, kj, = 0.003 
ft, k, = -0.003 ft, and C, is given by the simple equation, 


r h 
Ce = 3.22 + 0.40 B (25) 


From the Schoder and Turner tests (Series E and F) in Fig. 3(c), ky = 0.004 
ft, and 


Gs h 
Ce = 3.2) +0245 P (26) 


The results of the Bazin tests on full-width weirs are shown in Fig. 4. The 
value of k, required for the Bazin tests is considerably larger than that de- 
termined from the other experiments. From Fig. 4(a), kj, = 0.012 ft, and 


i h 
Ce = 3.25 + 0.44 5 (27) 


The results of the USBR tests (Fig. 5), in comparison with the others, are 
characterized by smaller values of C, at the larger values of h/P. For the 
USBR tests, kp = 0.003 ft, and 


Ce = 3.25 + 0.35 B (28) 


Rehbock’s 1929 formula is also based on the effective-head concept. 
Furthermore, his equation for C, is in the linear form of Eqs. 25 to 28, in- 


clusive. Thus, for comparison with those equations, Rehbock’s 1929 formula — 
gives 


C, = 3.22 + 0.44 B (al 


with kj, = 0.004 ft. Rehbock did not recognize the existence of kj, although 
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ome writers have indicated that his experiments were made ina very narrow 
lume. Therefore, it is possible that Eq. 29 would be modified if kp were taken 
nto account. This conclusion could not be verified because Rehbock’s experi- 
mental data were not available. 


| 
| 4.2 


| 
| Rehbock (1929) formula (Eq.29) 
| | | 
Schoder and Turner tests (Eq.26) 


Bazin tests (Eq.27 


~— 


Georgia Tech tests (Eq.25) 


= a 


tests (Eq.28) 
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Fic. 13- COMPARISON OF EXPERIMENTAL Data (b/B=I.0) 


Eqs. 25 to 29, inclusive, are compared in Fig. 13. The Schoder and Turner, 
3azin, and Rehbock curves agree reasonably well. The curve determined from 
he Georgia Tech tests lies below these curves and above the USBR curve. It 
vas suggested previously that the differences in the results of the experiments 
shown in Fig. 13 were caused by differences in the test equipment. Thus, the 
listinctive character of the C, curve from the USBR tests can be associated 
vith the fact that the crest piece for these tests was knife-edged. Similarly, 
he larger value of kj, determined from Bazin’s tests can be related to the 
Onclusion that his crest was not truly thin and sharp. Such details as these, 
requently unknown or uncontrollable, are believed to be the cause of the vari- 
nce shown in Fig. 13. Nevertheless, it is noteworthy that the equations for 
le derived from various sources are all simple, linear functions of the h/P 
atio. 

Fig. 14 shows a comparison of the results of the Georgia Tech tests with 
ome of the most widely used formulas for full-width weirs. The comparison 
; made for two values of h and a full range of values of h/P. The Bazin 
xrmula (Eq. 18), the Rehbock 1912 formula (Eq. 19), and the S.I.A. formula 
Eq. 20) are in the non-linear form suggested by Weisbach. For these formu- 
as C was computed as Q/bh3/2 from Eqs. 18, 19, and 20. For the Georgia 
‘ech tests and the Rehbock 1929 formula, C was computed from the relation- 
hip 


Qa. = Saboth + tey)??  b tik 3/? 


Tg AS eee h : (30) 
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Fic. 14.- COMPARISON OF FORMULAS (b/B=1.0) 
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which b is assumed to be so large that the influence of kp can be neglected, 
hence b>be. Values of Ce in Eq. 30 were computed from Eq. 25 for the 
reorgia Tech tests and from Eq. 29 for the Rehbock 1929 formula. 
| Fig. 14(a) shows a comparison of the different formulas for h = 0.2 ft. For 
nis condition the influence of the fluid properties and effects of differences in 
e experimental equipment are relatively large. Fig. 14(b) shows a compari- 
on of the formulas for h equal to 1.5 ft. For this condition the effects related 

small values of h are negligible (See Fig. 11). The curves in Fig. 14(b) 
how reasonably good agreement. 

In general, the comparisons of experiments and formulas shown in Figs. 13 
ind 14 cast doubt on the adequacy of some of the widely used formulas and on 
he reproducibility of the test conditions from which they were derived. In 
articular, the relatively complex Weisbach-type formulas are not justified 

comprehensive discharge equations. This conclusion is emphatically demon- 
trated in Fig. 4. From the family of curves shown in Fig. 4(a), it is apparent 
hat the Bazin formula compares unfavorably with the results of the Bazin 
ests for all but the intermediate values of h included in the experiments. In 
‘ig. 4(b), however, a single, straight-line curve fits the data reasonably well 
or all values of hand h/P. Thus, the simple, linear equation for Ce, comple- 
nented by the effective-head concept, is substantiated as an improved method 
f defining the discharge function. 


‘omparison of Formulas for Notch-Weirs 

Fig. 10, based on the Georgia Tech tests, shows the relationship between 
e, h/P, and b/B for notch weirs. A comparable family of curves can be 
lotted from the S.I.A. formula for notch weirs (Eq. 22) if the term involving 
he absolute head is omitted. For relatively large values of h (and b), the 
nfluence of the h term is negligible, and values of C from Eq. 22 can be com- 
ared with values of C, from Fig. 10. The comparison is shown in Fig. 15. 
the maximum difference between corresponding values of Cg shown in Fig. 15 
$s about 4 percent. BOWES, within the limits of applicability specified by 
he S.I.A. (b/B 20.3, h/P © 1.0), the maximum difference between the two fami- 
ies of curves is about 2 percent. It is emphasized that Fig. 15 does not in- 
lude a comparison of the different methods of compensating for the fluid- 
roperty effects. However, it is believed that the advantages of the effective- 
ead concept were clearly demonstrated in the previous comparisons of formu- 
as and experiments for full-width weirs. 
The results of the Georgia Tech tests on notch weirs (Fig. 10) indicate that 
‘e is a function of b/B and h/P. In the Francis contracted-width formula 
Eq. 23), however, a correction for “end contractions” is assumed to be a 
inction of h alone. Because of this fundamental difference, the Francis dis- 
harge formula (Eq. 24) cannot be compared with the Georgia Tech tests or 
1e S.I.A. formula on a graph such as Fig. 15, Nevertheless, it is possible to 
how that the Francis formula is not substantiated by the curves shown in Fig. 
{ Eq. 23 can be written in the form 


y 
: b, = b - Kh (31) 
which, according to Francis, K is a constant equal to 0.2 for “fully contract- 
d” rectangular-notch weirs. From Eqs. 11, 24, and 31, 


| 
| 


| 


| 
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| Q = Cebghe?!? = Coi(b - Khyn?/? ey 
a which C, is the coefficient of discharge for notch weirs and C,, is the corre- 
Onding coefficient for full-width weirs (b/B = 1.0). If the fluid-property ef- 
cts related to small heads and widths are neglected, h = he and b ~ be. 

hus, from Eq. 32, 


c 
RD CE (33) 
fig. 16 shows K as a function of h/P and b/h. Values of K shown in this figure 
vere computed from Eq. 33, using values of C, and Cg, from Fig. 10. Itis 
pparent that the constant value of K recommended by Francis is not substanti- 


ted, even within the limits specified by Smith and Francis for “fully contract- 
d” weirs. 


Value of K 


(Francis) 


eae K=0.2 


Curves based on 
Georgia Tech Tests 
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Fic. |6.- CoMPARISON WITH FRANCIS FormULA FOR NoTcH WEIRS 


The comparison shown in Fig. 15 indicates that the S.LA. formula for notch 
eirs agrees reasonably well with the results of the Georgia Tech tests. How- 
yer, the complexity of the formula is a deterrent to its use, and the limits of 
plicability specified by the S.I.A. prevent its being a comprehensive solution 
ra full, practical variety of notch weirs. Furthermore, the effective-head 
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and effective-width concepts which are complementary to the simple C, 
equations derived from the Georgia Tech tests are believed to be superior t 
any other method of accounting for the influence of viscosity and surface 
tension. 


CONC LUSIONS : 

The flow pattern for rectangular, thin-plate weirs is not subject to com- 
plete mathematical analysis. Consequently, an analytical solution for the dis 
charge characteristics has not been developed. A comprehensive solution 
based on dimensional analysis and experiment is presented in this report. I 
provides a simple, direct solution for the discharge and a convenient method 
of compensating for the influence of viscosity and surface tension. It is sub- 
stantiated by experimental data from various sources. 

Widely used formulas adapted from the classic Weisbach equation are not 
justified as comprehensive solutions. In fact, they do not satisfy the need fo: 
a simple equation which can be used to correlate experimental data covering 
a full range of the significant variables. Unjustified emphasis has been plac 
on the significance of the orifice analogy in the derivation of weir equations. 
Furthermore, there is no logical basis for the assumption that formulas invc 
ing the true total-energy head are independent of the h/P ratio. The Francis 
contracted-width formula is fundamentally inadequate as a means of adapting 
formulas for full-width weirs to notch weirs. 

Weir discharge is critically influenced by the physical characteristics of 
the weir and weir channel. It is especially dependent on features which cont 
the velocity distribution in the approach channel. For this reason the result 
of experiments made by different, capable investigators do not agree, and 
formulas based on a particular set of data reflect the individual characterist 
of those data, For this reason, too, a truly reproducible, standard measurinr 
weir and a precise, universal discharge formula are impractical. 
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SYNOPSIS 


Air Binding is a potentially serious cause of excess head loss in water 
pipelines flowing under less than atmospheric pressure. This paper dis- 
cusses the phenomenon of air binding and methods of control, with emphasis 
on the accumulation of air in sloping pipes. The specific installations 
described are circulating water systems of steam-electric stations, but the 
principles discussed are applicable to any water pumping station with similar 
profile. 


INTRODUCTION 


Air binding in a pipeline is the trapping of air within the water passage in 
1 manner which prevents the pipe from flowing full. In poorly laid out or 
J0orly vented pipelines air binding may be the source of head losses which 
“an Seriously reduce the pipeline capacity. 

_ The accumulation of air is particularly severe in a pipeline flowing as a 
siphon under less than atmospheric pressure. Such a condition may occur, 
or example, where the pipe passes over a hill. 

Important water pumping installations which normally utilize a siphon are 
he circulating water systems of steam-electric generating stations. These 
umping plants may -move as much as 150,000 gallons per minute through the 
ondenser of each large generating unit. Large-diameter circulating water 
ipelines, several hundred feet long, are often designed to flow under a 
acuum between minus 20 and minus 30 feet of water. The phenomenon of air 
inding is illustrated in this paper by two examples of circulating water sys- 
ems which were seriously obstructed by air. Excess head losses due to 
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entrapped air were 15 feet and 19 feet for systems designed originally to 
pump against normal pressures of 43 and 36 feet. 

It will become apparent by study of these examples that the air binding 
problem may be applicable to any other pipeline with similar profile. 


Basic Elements of the Air Binding Problem 


Air accumulation in a pumping system may originate from several 
sources: air exists in the empty pipeline before it is filled, air is entrained 
in the water carried into the system and, at sections under negative pressure, 
it can leak in at joints and fittings. This air will be present in varying con- 
centrations from very small, widely dispersed bubbles to large pockets. 

The air is carried along the pipeline by the moving water. However, wher 
the pipe slopes downward from a high point, the bubbles tend to collect along 
the top of the pipe. When the buoyancy of the air accumulation becomes 
greater than the water force pushing it down the pipe the pipeline becomes 
air bound. 

The most severe head losses encountered in circulating water systems 
flowing under vacuum are the result of large air pockets lying along down- 
ward slopes and creating a condition similar to open channel flow. The ener- 
gy gradient thus roughly parallels the actual slope of the pipeline instead of 
representing only the friction losses in a pipe flowing full. Figure lisa 
diagrammatic view of the phenomenon of airbound slopes. The air pocket 
terminates in a hydraulic jump which absorbs further energy. 

Generally speaking, for a given water velocity and pipe slope air binding 
difficulties increase with the pipe size. For example, an 8-foot-per-second 
velocity may be sufficient to carry all accumulating air out of a 12-inch pipe 
on a given slope, but is not sufficient to move such air through an 84-inch 
pipe on the same slope. For large-size pipe there is very little published in- 
formation relating the factors of air entrainment, air release, water velocity, 
pipe size, pipe slope and natural air removal. 

The author has been concerned particularly with conduits of such a size 
that water velocities up to approximately 10 feet per second have not been 
sufficient to move air down a slope. Pipes of five feet or more in diameter 
fall into this class. The discussion is, of course, applicable to smaller con- 
duit where water velocities are low. 

Experience indicates that where a large diameter pipe normally flows un- 
der vacuum there will be air binding in any non vented section which slopes 
downstream from the high point in the system. Accumulating air must be 
continuously removed or the efficiency of the system will be impaired. Ex- 
ceptions to this general rule have been at slopes immediately following sharp 
horizontal bends; the turbulence induced by the horizontal bend may break up 
large air pockets into more easily transported smaller bubbles. 


Self-Priming Vertical Drops 


When a change in grade is made with a vertical or nearly vertical drop the 
possibilities of air binding are much reduced. Such a vertical drop is gener- 
ally “self-priming.” The turbulence created by the abrupt change in water . 
direction breaks up and disperses large pockets of air. The smaller bubbles 
are then carried away by the moving water. See Figure 2. 
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Vertical Versus Sloping Pipelines 


| it is apparent from the comments above that air evacuation equipment 
must be used unless the change in grade is made with a vertical drop. This 
equipment, discussed briefly below, is expensive and requires constant main- 
enance. However, vertical drops have disadvantages also, and in many cases 
ay be more costly to install and operate than a continuously evacuated slop- 
ing pipeline. 
_ First, the turbulence which entrains and carries away the air at a vertical 
ie will also create a permanent loss of head which raises the cost of pump- 
ng water. This loss for one vertical drop with two reverse 90° bends may 
umount to 1 foot of head and represent an annual excess pump operating cost 
of several hundred dollars and a capitalized cost of several thousand dollars. 
One company solved a potential air binding problem in one of their older 
lants by traversing a downgrade with a series of three vertical drops, a total 
of six 90° bends. At a large plant the cost of head lost in such an arrange- 
ment would be considerably more than the cost for the installation and main- 
enance of an air removal system on a sloping pipeline. 

A second serious objection to the vertical drop is the thrust produced at 
109 bends. Massive thrust blocks must be provided to anchor such a bend 
securely in the event of water hammer surges. At one large power station 
he thrust block required at a vertical drop in twin 84-inch conduits contained 
wer 400 cubic yards of concrete and required complicated formwork. A new 
‘irculating water system at this same station has a sloping conduit equipped 
vith an air removal system at a fraction of the cost of the thrust block. 


Methods of Air Removal 


_ Air taps at high points are standard practice on pipelines, but it is appar- 
nt from the previous discussion and a knowledge of laboratory observation 
hat the air carried by moving water will be found along the downstream 
lope. Air removal systems on the conduits discussed below were not effec- 
ive when they were connected to the high point taps alone. Accordingly, a 
ystem of taps along the downstream slope was provided. These taps were 
Onnected to a common header which in turn was connected to the air removal 
ystem. 

Where a vertical drop is used a single air removal connection should be 
laced at the high point at the outside edge of the bend. Generally, this tap 
ill not be required during operation unless water velocities are very low. 

Where an airbound pipe flows under vacuum some type of air evacuator 
‘ill be required to remove the air accumulations. Such an installation may 
onsist of one or more vacuum pumps, a control system, power supply and, 
or water-sealed pumps, a continuous supply of clean water (or a closed- 
ircuit sealing system using some other fluid). The author’s firm has recent- 
, installed a water jet evacuation system which is considerably less complex 
4 costly than vacuum pumps. Operating water for the jets is supplied from 
le power station elevated water tank, but could be supplied by booster pump 
king suction from the pipeline itself. 

_ The primary consideration in sizing air removal equipment is to provide 
hough capacity to maintain the line free of air during operation. A second- 
ry consideration is the rate of initial priming of the system. 


| 
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The determination of capacity requirements is a difficult matter. Uncer- 
tainty is caused by several factors. A first unknown is the quantity of air 
which may leak in through joints and fittings in long conduits under vacuum. 
Also, experimental data for these large-size pipes is not sufficient to permit 
determination of the amount of entering air which will be carried out natural- 
ly for given velocities, slopes and pipe Sizes. 

It may be desirable to install one fairly large vacuum pump for initial 
priming, then a smaller pump in parallel to maintain the prime. The smaller 
one can be selected and installed after actual operating experience with the 
large vacuum pump has indicated the air binding characteristics of the pipe- 
line. 

The details of two typical circulating water systems are described below. 
Both of these are older stations at which air binding problems had to be cor- 
rected. For all new stations the original design should provide for air re- 
moval if the profile indicates a potential problem. 


Station A - 132,000-KW Power Plant Installation on a Large River 


Station “A” has a 1,750-foot condenser discharge conduit, all of which 
should normally operate under a partial vacuum. The majority of the pipe is 
78-inch diameter, the design velocity is 7.3 feet per second. The profile of 
this conduit is shown on Figure 3. The head recovery in this discharge 
siphon line was very poor and caused considerable concern because of the 
loss in water pumping capacity. 

The vacuum measured at the condenser discharge should be minus 22.9 
feet of water for river elevation 718. The actual measurement at the north 
condenser was minus 7.9 feet, an unexplained loss of 15 feet. It was deter- 
mined that with river water levels at 719 or above air accumulated along the 
57:1 sloping section shown on Figure 3, giving a head loss at water level 719 
of about 7.5 feet. This represented the approximate difference in elevation 
between the top and bottom of the slope, and indicated that the pressure 
gradient was parallel with the pipe slope as could be expected for open chan- 
nel flow down a similar incline. At water levels under 719 the air pocket als 
extended along the entire top of the horizontal section and into the 459 sloping 
condenser discharge pipe, thereby reducing the siphon an additional 7.5 feet. 


When the river level was at 709, flow through the pipeline was recorded as 
follows: 


Actual Flow 


Design Flow in Pipeline Total Loss 

With Proper Venting No Venting of Capacity 
Two Pumps 70,000 gpm 58,000 gpm 12,000 gpm 
Three Pumps 89,000 gpm 74,000 gpm 15,000 gpm 


The points of air evacuation had been installed in accordance with origina 
design, but proved to be inefficient. A 95-cfm, 15-inch Hg vacuum pump was 
located at Point A where the discharge line passed the intake structure. It 
worked continuously and was of insufficient capacity. Two-inch air-powered 
ejectors at the condensers were also ineffective. 

As a result of investigation, a 200-cfm vacuum pump was installed at 
Point B on Figure 3. Additional taps were also installed on either side of 
Point B, and on the upper end of the 45° condenser piping. These all take 
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uction from the same pump. The 95-cfm pump at the intake structure was 
Iso replaced with a 200-cfm pump. 


| This arrangement was effective in eliminating the air binding problem. 
Station “B” - 166,000-KW Power Plant Installation on a Large River 


_ The circulating water system of this station has three pronounced high 

oints, the usual one at the condensers, and two where the intake and dis- 
charge conduits pass over a flood protection levee. The principal airbound 
section is 66-inch diameter, design velocity is 10.7 feet per second. An ap- 
roximate profile is shown on Figure 4. 

The pressure at the crest of the discharge levee crossing was about minus 

-) feet without venting. The design pressure was minus 23 feet, a difference 
bf 19.5 feet representing excess lost head. A vacuum pump taking suction 
irom a Single tap at the crest had been able to maintain a maximum vacuum 
of slightly less than 15 feet, even though the gage glass at the air receiving 
ank showed full prime at this point. Taps were therefore drilled at three 
udditional points on the downstream slope of the pipe and connected to the 
yacuum pump. Maximum design siphon of minus 23 feet was obtained within 
ess than one hour and could be maintained indefinitely with a 6-cfm pump. 
[he air had been trapped on the slope beyond the levee crest as could be ex- 
yected from previous discussion in this paper. 

With the increased vacuum provided by the air evacuation system four- 
jump flow in the system was increased from 115,000 gpm to 134,000 gpm. 
[This increased flow was so great that the pump impellers had to be rede- 
signed to prevent excessive impeller pitting due to cavitation. 

The amount of air to be removed from this system was very small since 
ittle leakage of air could take place. The only section under vacuum was 
hat section of welded steel pipe over the levee. 

When only one or two of the four circulating water pumps were operating 
it low water conditions a slight vacuum also occurred at the intake levee 
rest and a similar air-binding condition arose. Vacuum pump connections 
yere run to this point also. A completely separate parallel circulating water 
ystem has recently been installed for Units 3 and 4 of this power plant. A 
ank of vacuum pumps with interconnected controls and suction lines draws 
ir from both circulating water systems at the intake and discharge levees. 
‘he maximum length of air suction line is about 600 feet. 


CONCLUSION 


These two examples clearly indicate the extent of capacity loss which may 
e experienced in pipelines flowing downhill under a vacuum. This phenome- 
on probably occurs much more frequently than records show. In many 
team plant circulating water systems overdesigned pumps have made up for 
ick of hydraulic efficiency. Air binding may well have gone unnoticed at 
tation B if the extent of head loss had been a little less severe. At Station A, 
n the other hand, air binding threatened to reduce the efficiency of the eoure. 
ower station as summer weather brought a greater demand for cooling 
ater. 
There is a particular need for experimental air-binding data applying to 
irge diameter pipelines. Observations at Station “B” for several different 


—— 


pumping conditions have indicated that the pipe may be self-evacuating at 
velocities slightly above the design. Empirical formulas based on enough 
carefully evaluated field data may well permit the confident design of a self- 
evacuating sloping pipeline flowing under vacuum. 
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FLOW THROUGH CIRCULAR WEIRS 


J. C. Stevens,! M. ASCE 
(Proc. Paper 1455) 


SYNOPSIS 


There are certain advantages inherent in the use of the circular weir. How- 
ver it has not come into general usage in this country. It is more extensively 
sed in Europe where there has been developed a rigid, accurate formula by 
hich the flow can be determined. 

There is presented herein a formula for the theoretical flow through sharp 
rested circular weirs which involves the use of elliptic functions. Most of 
1e European treatises on such weirs are not generally available now. Tables 
ppended give (1) a summary of available experimental data with discharge 
efficients, (2) the flow through weirs of various diameters based on these 
defficients and (3) values of the elliptic functions involved by which the flow 
lrough a circular weir of any diameter may be found. 


INTRODUCTION 


Weirs with circular crests have the advantage that the crest can be turned 
id beveled with precision in a lathe. Moreover the weir crest does not have 
) be leveled, hence the point of zero flow is readily determined. 

When wells are tested by pumping the pump discharge is usually measured 
J an orifice fixed to the end of the discharge pipe. The flow is determined 
7 the orifice formula in which the down-stream head is constant—that of one 
mosphere. When the pipe is not running full the orifice becomes a circular 
eir and it is often important to know the flow under such conditions. 

The problem of circular weir discharge first thrust itself upon the author’s 
tention in connection with a rather novel fishway where the discharge from 
was to be carried in a pipe to a power canal a short distance away. This 
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flow had to be measured. The circular orifice in the end of the pipe would be 
running partly full most of the time, and therefore became a circular weir. 


Experimental Data 


A few laboratory experiments are available on the flow through circular 
and semi-circular weirs. A summary of all such experiments the author coul 
find appears in Appendix A, Doubtless others were made in European countri 
but are not available to the author. There is also given in Appendix B a table 
of flows through such weirs from 0.25 to 5.0 feet in diameter. The theoretica 
flow was first determined by the rigid formula presented herein which was 
then modified by an average discharge coefficient disclosed by the above labo- 
ratory experiments. Appendix C presents a table of the elliptic function in- 
volved in the formula by which the theoretical flow for a weir of any diameter 
can easily be tabulated. 

The following laboratory experiments to determine the flow through circul: 
and semi-circular weirs have been found in available technical papers. 


Experiments by H. V. Cone(1) 

Experiments on the flow through two circular and two semi-circular weirs 
are reported. The results are given graphically by 4 curves. No attempt was 
made to develop a flow formula for them. 


Experiments by F. V. Greve(2) 

Extensive experiments on the flow through circular, parabolic, and triangu 
lar weirs were conducted by the late F. V. Greve at the Engineering Experi- 
ment Station of Purdue University, Lafayette, Indiana. 

The differential formula he developed for the flow through such weirs, he 
states, could not be solved but the flow through semi-circular weirs was foun 
by mechanical integration of his differential equations for several diameters, 

The flows in Greve’s bulletin are all given in pounds per second but they 
have been converted to cubic feet per second herein. 


Experiments by Eldon R. Dodge(3) 

These were made to determine the flow through circular orifices in steel 
plates fixed to the outlet end of sewer pipes intended for storm drain culverts 
The orifice plates were attached to the pipe with various eccentricities, that 
is with the centers of some orifices concentric and others not concentric with 
the centers of the pipes. Only the flows for concentric orifices are tabulated 
herein. Flows were determined under heads with orifices flowing partly full 
and also under surcharged heads wherein the discharge changes from weir 
flow to orifice flow. 

Mr. Dodge presents a trigonometrical series formula, aided by curves, by 
which the theoretical flows through his weirs were approximately found. He 
also presents curves of his discharge coefficients. 

It is doubtful whether any of the above authors and experimenters were 
aware of the work done in Europe on flow through circular weirs. It is prac- 
tically certain that the first one in this country to make recorded tests of suc 
weir flow was H. V. Cone who made his experiments in 1916 at Fort Collins, 
Colorado. Dr. Eng A. Staus’ first publication of which the writer has knowlec 
appeared in 1926, while G. H. Gulliver made flow tests in England on small _ 
circular weirs as early as 1909. (See References) 


CE STEVENS 1455-3 


xperiments by G. H. Gulliver(4) 
These were made on a small circular weir 2-1/2 inches diameter. The 


ifferential formula developed is 
: 


h 
: q = 12.05 [ (a(R-h) [D-(H-h)] dh 


where D = diameter of the circular weir-inches; H = head above the lowest 
ov of the weir-inches, q = discharge - gpm; h = distance to elementary area 
ve the circular weir invert. The formula was mechanically integrated and 
e results are given by curves. Flow coefficients found were between .625 
d .630. The curves were on too small a scale to be read with accuracy and 


is results are not included in this paper. 


Experiments byob, Thijsse(9) 

Circular weirs have been in use for many years in his laboratory to 
measure the flow through hydraulic models of structures and other devices. 
They were adopted because of their ease of construction but more particularly 
because they did not have to be leveled and the point of zero is so readily de- 
termined. 

The formula given is 


Q=m VgR = kH2 


Where Q, flow in liters per second 
g, the metric gravity constant 
R, radius of Circular Weir - centimeters 
H, head on the Weir invert - centimeters 
m, discharge coefficient 
The value of m was found by calibration in place. Curves were prepared for 
values of k vs H2 from the calibration for each size of weir. These curves 
however have no general validity if the conditions of use depart materially 
from those under which they are obtained. 
These data are given in metric units but they have been converted to English 
units in the tables herein. 


The Elliptic Formula 


- The articles by Dr. Eng Staus do not give the transformation to the Elliptic 
orm which is rather involved. This however was developed entirely inde- 
endently by the author with substantial assistance from Dr. W. E. Milne, re- 
ired, formerly head of the mathematics department of Oregon State College. 
[his transformation is given in Appendix D. 

Referring to Fig. 1 it is obvious that 


2 Yi eae = 2N x 
: 1=2 (2) : - y) fy) = y) (1) 
j 
Phe head on the elementary area 1 dy is H - y and the elementary theoretical 
low becomes 


} dQ,= 2 \2e \(O-y) y (H-y) dy 


his is true whether H is greater or less than the radius of the circle. Then 


} 
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i 
ye = Z e Aig - )y¥ (i - y) dy (2) 
Q, 2 \fee {Ni y) ¥ SS 


This can be developed into a series and integrated(6) term by term obtain- 
ing as was done by E. R. Dodge but not in a trigometric form 


bt hee Mia et Re peed > 
Q, =2 veg 2°/4n(a2*-42?-gegz*----) pierce 214, 21m 


Where z = H/D 
This, however, is an infinite series and laborous to deal with. 
The elliptic formula for the theoretical flow through circular weirs is 


4, 1.2) ~2) ze 
a, = Ayaan’? [20 - 4 k JE = (2-k*)} (1-k5) x | (3) 


Application of the Elliptic Formula 


In equation (3) the modulus k depends solely on the head and diameter of 
the weir. In brief the formula may be written 


Q =BA 
where 
B = 4/15 V2gD9/2 = 2.14p9/2 


and Ais the quantity in[ ] 

In evaluating (3) for various diameters a table of k2 vs A was prepared (se 
Appendix C). This table gives A to 4 decimal places for values of D/H to 3 
decimal places. If we multiply the values of A by 2.14 we obtain the theoretic 
flow for a 1-foot diameter weir since 


B = 4/15 Vag D9/2 = 2.14 


The value of the factor B for any other diameter is obtained by 
B = 2,14D5/2 


Discharge Coefficients 


The tabulation of experimental data in Appendix A gives the discharge co- 
efficients obtained by 


c = Q/Q 


These values are plotted in Fig. 2. While there is some evidence of a non 
linear trend it is not definable in the evidence at hand. No coefficients have 
been discarded. The average of each experimenter’s coefficients for each we 
diameter is also indicated in Fig. 2. From this it appears that an average co 
efficient of 59 per cent is justified for the experimental data presented. 

The coefficients from the Thijsse experiments indicate that the flows wert 
influenced by velocity of approach. However this would not affect the validity 


a STEVENS 1455-5 


oi his experiments since his weirs were all calibrated in place and the flows 

| observed include this factor. The Dodge experiments also appear to show 
|} this same effect since the orifices were in the end of drain pipes. 
The experiments by H. V. Cone are all near the 59 per cent line except 
| those for the 3-foot weir for which there is no evident explanation unless his 
zero point was in error. 

The Greve experiments are the most consistent of all and are believed to 
_be entirely trustworthy. 

It is hoped that some laboratories in this country will make more experi- 
ments on the flow through circular sharp-crested weirs and thereby obtain 
more discharge coefficients. 


Notation 


In this paper the letter symbols for Hydraulics, approved by the American 
Standards Association January 1942 has been used as far as practicable. 
The foot-seconds system has been used throughout. See also Fig. 1. 


iD, diameter of sharp crested weir opening 


H, head on weir above the invert measured upstream of the surface curve 


of approach 
Q:, theoretical flow through weir 


Q, observed flow 

(ee coefficient of discharge = Q/Q 

g, acceleration of gravity 

1; length of chord of elementary area 


y; distance above circle invert to an elementary area 
k, modulus of the elliptic integrals = VH/D 

K, the complete elliptic integral of the first kind 

E, the complete elliptic integral of the second kind 
A, Value of the elliptic function involved such that Qt = BA and Q = BCA 
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.211 


.25 


131 


133 
137 
.146 
176 
«L79 
.240 
244 


.336 


.128 


142 
-161 
.168 
.219 
229 
235 
-243 
-262 
.273 
.294 
.299 
302 
321 


E. R. Dodge 
H/D Q, 
.398 .018 
474 = 025 
-ol2: 029 
-626 .042 
.677  =.048 
«111 =. 052 
120 053 
.802 .063 
.810 .064 
-839 .068 
696 075 
915 .078 
-958 .082 
1.000 .088 
F, W. Greve 
029 86.046 
.933 = .048 
.048 .050 
576 .055 
-704 = .078 
311602, 080 
960 .127 
.976 .129 
E. R. Dodge 
-380 §=.054 
-423 = .069 
-478 =.085 
-500 .089 
652 .144 
.682 155 
-700 =—.162 
123 174 
COO MmEsL OD 
-812 208 
875 233 
888 .239 
-900 .244 
-955 .265 


Q 


.012 
.016 
.019 
.026 
-031 
.031 
-032 
.036 
.039 
.043 
-046 
-050 
.051 
-055 


.027 
-028 
.029 
-032 
.045 
.047 
.074 
.076 


-034 
-040 
.050 
.055 
-084 
.097 
.096 
101 
114 
.122 
-136 
.139 
151 
.156 


C 


666 
.640 
.655 
.619 
.646 
.096 
.604 
ae: 
.619 
.633 
.614 
.640 
.623 
625 
624 


.087 
.584 
.580 
.082 
DTT 
.588 
.583 
-590 
585 


-628 
.580 
.588 
618 
.583 
625 
.593 
.580 
584 
.087 
.585 
.983 
-619 
589 
-602 
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D 


-50 


-90 


50 


-20 
25 
30 
35 
-40 
45 
-90 


145 
.219 
.259 
317 
-376 
-407 
-461 
-467 


F, W. Greve 
H/D Qt 
-202 5) 000 
.314 .101 
s300 105 
412 .168 
-424 177 
.046 .282 
-610 .344 
.688 .424 
S1US ona 
880 .633 
H. V. Cone 
-40 e459) 
.50 -241 
-60 .334 
.70 -437 
.80 .045 
.90 -655 

1.00 spe 

E. R. Dodge 
-289 §=.085 
-437 = 189 
sol6 2256: 
<033° 2.369 
5150). 3493; 
-814 .563 
-920 .680 
-933 .694 
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Q 


.043 
.060 
.061 
-098 
-104 
.164 
-200 
.244 
.304 
.368 


-10 
13 
19 
-25 
.30 
.38 
44 


-054 
113 
154 
-218 
-297 
-340 
421 
421 


6 


-606 
.594 
.580 
.583 
.588 
.080 
.580 
A ds, 
.580 
580 
.585 


628 
540 
-568 
212 
.550 
.980 
580 
585 


.635 
598 
-602 
592 
.603 
-605 
619 
.606 
-607 


15 


333 
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1.00 


1.167 
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— 


.293 
-316 
-390 
-492 
.603 
- 108 
-825 
947 
1.166 


15 
-30 
52 
-80 
1.10 
1.41 
1.76 
2.15 
2.50 


247 
377 
.953 
.870 
1.300 
1.684 
2.168 
2.700 
3.788 
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J. Th. Thijsse F, W. Greve 
D H H/D Qt Q C D H H/D Qt Q Cc 
138k .236 2200 .363.° 2210, “Jo78 1.50 214 143° 5345 3207 _. 260g 
.20 212 .406 .246 8.605 226. .151 .384. .230) | 2508 
.30 204 .0%6_ 350 607 248 1.165 .456 272 596 
.30 2296, Solio Gaia, sOLZ sold. «182. 001, Gos0L aoe 
.40 096 990. O13. 616 .309. 206 .700 .417 .596 
45 .o01 1.246-" N67 ~ 2616 -349 .232 .879 .526 .598 
.50 424 1.521 .944 .620 -388 .259 1.086 .643 .592 
200 -466 1.813 1.155 638 424. .283 1.289 755 3586 
-60 .007 2.116 1.365 .645 470. .314 1.574 .916 582 
.65 -550 2.452 1.590 .649 503... 2335 15779. 1.035 Sooee 
-70 -5992 2.796 1.860 .664 048 =.365 2.091 1.205 .576 
ths .635 3.163 2.110 .666 7990 .393 2.401 1.383 576 
.80 -676 3.526 2.410 .680 642 .428 2.815 1.620 S576 
-630 674 .449 3.077 1.780 .579 


-720 .480 3.481 2.003 .575 
920 .614 5.422 3.070 .567 


F, W. Greve 1.050 .700 6.804 3.800 .559 

1.196 .798 8.456 4.780 .565 

1.333 .280 .210 .540 .322 .597 1.406 .938 10.830 6.250 .576 
322 .242 .709 .417 .588 -982 


415.311 1.149" ~/675" 2588 
:495"° 372) 1.611 "935, 2580 
O17 =.433 2.139 1.244 .580 1-748: 3159.5 .091 2204557 128 Weeeoe 


582 .437 2.176 1.260 .580 207 118 .345 .2145 (62) 
.703  .527 3.069 1.797 .586 -204 145 .520)  .3167 “608 
-814 .610 3.986 2.296 .575 -321 .184 .825 .482 .584 
967 .726 5.380 3.055 .568 -388 .222 1.183 .710 .600 
1.124 .844 6.873 3.900 .567 435 1249 1.476 .875 .593 
1.146 .860 7.077 3.995  .565 487 .279 1.837 1.080 .589 
1.332 1.000 8.780 5.070 .578 -528 .302 2.141 1.250 .584 
-580 -588 .345 2.756 1.540 .560 


647 .370 3.143. 1.830 582 
-703 .402 3.670 2.140 .583 


H. V. Cone .750 .429 4.142 2.390 .577 
.801 .458 4.677 2.710 .579 

1,50. "220. 22185) SO s0gy in” dap .853 .488 5.257 3.050 .580 
25 .167 .467 25 535 .907 .519 5.883 3.420 .582 

80" *.200.” 660 9.38 875 591 


35 .233 ©.886 = .51 .975 
40 -267 1.152 .69 .998 
-50 -333 1.758 1.03 -086 
.60 -400 2.480 1.45 584 
-70 -467 3.308 1.90 .075 
-80 -533 4.211 2.44 .580 


J. Th. Thijsse 

H H/D Qa 
223 .126 ©.406 
25 st4aP 507 
30 69a ri22 
40 -226 1.263 
50 pap 2 eto) 
60 -339 2.749 
70 -395 3.663 
80 -452 4.707 
90 908 5.835 
00 .965 7.068 

F. W. Greve 

paso) 092 "293 
-231  .118 .483 
sai3 .136 +640 
.3845 .172 1.014 
-396 .198 1.328 
-470 .235 1.846 
POL Oe eee Odea. LOD 
sd0%  _.288 2.333 
sOl3—  .306, 3.071 
-661 .330 3.544 
-730 .365 4.286 
-7178 .389 4.828 
-873  .436 5.971 
-960 .480 7.134 
-988 _.494 7,520 
e2c0) 124 715 
foot, 156 25128 
413 .183 1.534 
70290 .238 25441 
627 .279 3.451 


Q 


-246 
. 306 
-437 
. 768 
1.185 
1.700 
2.330 
3.060 
3.910 
4.920 


-186 
- 300 
.392 
-615 
-192 
1.097 
1.305 
1.560 
1.808 
2.070 
2.504 
2.810 
3.485 
4.223 
4.430 


-450 

.690 

.840 
1.50 
2.02 


STEVENS 


Cc 


-606 
.604 
.605 
-608 
.613 
.618 
.636 
.650 
-670 
695 
-630 


.635 
622 
612 
.607 
096 
-595 
095 
Oe 
-589 
Ais} }5) 
085 
.583 
.583 
-092 
589 
.095 


.630 
-611 
-548 
614 
585 
-592 


D 


2.36 


2.495 


3.00 


J. Th. Thijsse 
H H/D Qt 


-223 .095 473 
20 -106 .590 
30 -127 = 846 
40 -169 1.484 
-90 -212 2,294 
60 -254 3.248 
-70 -296 4,362 
-80 -339 5.646 
-90 -382 7.066 
1.00 -424 8.583 
1.10 .466 10.23 
1.20 .908 11.98 
1.30 -950 13.85 
1.40 -994 15.87 
1.50 -635 17.85 
1.64 -695 20.86 


F. W. Greve 


311 .125 .942 
-409 .164 1.607 
453.181 1.943 
-490) 196 2.265 
-026 .211 2.608 
-563  .226 2.978 
-985 .235 3.209 
-620 .248 3.560 
-665 .267 4.105 
699 .280 4,498 
-723  .290 4.814 


H. V. Cone 


-30 -100 = .953 
-40 .133 1.686 
50 -167 2.640 
.60 -200 3.730 
-70 .233 5.006 
80 .267 6.506 
-90 -300 8.142 
1.00 .333 9.932 


Q 


-283 
354 
.907 
.894 
1.380 
1.975 
2.680 
3.470 
4.360 
5.430 
6,580 
7.900 
9.330 
10.90 
12.60 
15.20 


.073 

. 966 
1.175 
1.359 
1.550 
1.765 
1.895 
2.118 
2.424 
2.648 
2.823 


.48 

-81 
1.22 
1.72 
2.30 
2.95 
3.67 
4,44 


1455-13 


GS 


-600 
.600 
.600 
.603 
-603 
-608 
FOS) 
-615 
-617 
-633 
.643 
.659 
2675 
685 
. 706 
728 
636 


.608 
.601 
604 
-600 
.094 
093 
.090 
.094 
590 
.089 
566 
.594 


.504 
.481 
.463 
-462 
-459 
454 
452 
447 
-467 
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APPENDIX B 


Flow Through Circular Weirs as Calculated from the Following Formula 


Q = 1.262 D9/2 2(1 - k2 - k4) E - (2 - k2)(1 - k2) K | 


where 

Q is the flow in ft3/sec 

D, the diameter of the Circular Weir - feet 

H, the head on the Weir invert - feet 

xd. H/D the head in terms of diameter 

E and K taken from tables of complete elliptic integrals 

b=4 C Vag D5/2 where C = .59 

15 
based on the experiments in Appendix A 
Values of b 
D p?/2 b D p°/2 b 
25 .0312 .039 2.00 5.656 7.14 
50 1768 223 2.50 9.882 12.47 
15 .4871 615 3.00 15.588 19.67 
1.00 1.000 1.262 4.00 32.000 40.38 


io) 
ie3| 


27] 
TT 


Cy 6 IT 


ee 


OFNOOODHLNOOD & wy 


Na A A. Nf ef te 


crn 


STEVENS 1455-15 


D= 50 9D =.75 D= 1.00). D = 1.50 D= 2:00) D=2.50 Di = 3:00 


H Q Q Q Q H Q Q Q 
-05 .006 .008 .009 O11 el .051 .054 .063 
10 = .025 .031 -036 .045 2 .204 227 .202 
15.054 -069 .081 .099 3 458 .014 .563 
20 .094 121 141 .176 4 WARE .908 .995 
29 86.142 .184 217 2210 -2 1.227 901.395 15558 
SOc Od .259 308 .389 [Onl 44eee 1.061 220! 
35 257 .346 -413 -022 of 2.336 2.667 2,954 
40 .321 -441 -930 .679 -8 3.000 3.445 3.840 
45 8 ©.385 .044 .661 -850 9 3.739 4.303-4.305 
00 446 .655 -802 1.036 1.0 4.535 5.241 5.862 
.05 _.769 953 1.244 1.1 5.390 6.259 17.038 
-60 -888 1.113 1.462 1.2 6.294 17.352 8.267 
-65 1.008 1.281 1.694 1.3) “7.245 (8.511 9.583 
-70 1.124 1.454 1.950 1.4 8.222 9.725 11.03 
15 1.235 1.633 2.210 1.5) 95237 921500) 12:50 
-80 1.815 2.482 1.6 10.26 12.32 14.04 
.85 2.000 2.774 1.7 10.98 13.68 15.69 
-90 2.181 3.067 1.8 12.33 15.06 17,35 
95 2.357 3.369 129) 13.33 16:50 19:06 

1.00 2.524 3.690 2.0 14.28 17.93 20.87 
1.05 4.007 2.1 19.38 22.66 
1.10 4.330 2.2 20.82 24.49 
1.15 4.671 2.3 22.25 26.42 
1.20 5.001 2.4 23.63 28.29 
1.25 9.335 2.5 24.94 30.17 
1.30 5.676 2.6 32.10 
1.35 6.010 2.7 33.99 
1.40 6.333 2.8 35.82 
1.45 6.657 2.9 37.65 
1.50 6.956 3.0 39.34 
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D=4.00 D=5.0 


a 


Q 


.081 
.295 
.646 
1,155 
1.809 
2.592 
3.497 
4.518 
5.669 
6.941 
8.342 
9.865 
11.49 
13.22 
15.04 
16,97 
19.00 
21.25 
23.36 
25,66 
28.04 
30.49 
33.02 
35.61 
38.26 
40.99 
43.73 
46.52 
49.38 
92.25 


OARDARWNKPOHODARTUAWNHODODARTAWNE 


BNO (38) BS SD GS) Do) ae esi Gin oe ee 


CoS 
oOo wo 


Q 


.092 
.325 
.720 
1.284 
2.018 
2.907 
3.951 
5.136 
6.448 
7.895 
9.482 
11.21 
13.08 
15.09 
17.24 
19.49 
21.86 
24.35 
26.94 
29.65 
32.47 
35.41 
38.46 
41.60 
44.83 
48,15 
91.54 
95.22 
58.58 
62.21 


HY 6 


DeMaria hmastet So ne OO: P88. OO: LO G9 
SCODARMNARWNHOHOAARIPwWNYH 


D = 4.00 


December, 195' 


D = 5.00 


Note: Although Q is given for the complete circle of the weir, the flow begin: 
to increase at about 70% of the diameter beyond that given by the 
formula, thus, partaking of a transition between weir flow and orifice 


flow as the weir becomes surcharged. 
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APPENDIX C 


Values of A = 2(1 - k2 - k4) E - (2 - k2)(1- k2) K 
where k2 = H/D; E and K to be taken from tables 
of complete elliptic sa tag py 


. 


WAT 85-¢ 
ee: for ¥2 values the 8a decimal is a te top. a 
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» ae ’ : 
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io 
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k2 0 1 2 3 4 5 6 7 8 9 Diff 
for 2 
= .001 

OP «0 0000 .0001 .0001 .0002 .0002 .0002 .0003 .0003 .0004 0.4 


01 .0004 .0005 .0006 .0007 .0008 .0009 0010 .0011 .0011 .0012 0.9 
02 .0013 .0015 .0016 .0018 .0019 .0020 0022 .0023 .0025 .0026 1.4 
03 .0027 .0028 .0030 .0032 .0034 .0035 0037 .0039 .0041 .0043 1.9 
04 .0046 .0048 .0050 .0053 .0055 .0058 0060 .0063 .0065 .0068 2.5 
05 .0071 .0074 .0077 .0080 .0083 .0087 0090 .0093 .0096 .0099 3.1 
06 .0102 .0106 .0109 .0113 .0117 .0120 0124 .0128 .0132 .0135 3.7 
07 .0139 .0144 .0148 .0152 .0156 .0160 0165 .0169 .0173 .0178 4.3 
08 .0182 .0187 .0192 .0197 .0202 .0206 0211 .0216 .0221 .0226 4.9 
09 .0231 .0236 .0242 .0248 .0253 .0258 0264 .0270 .0275 .0281 5.5 
10 .0286 .0292 .0298 .0304 .0310 .0316 0322 .0328 .0334 .0340 6.0 
11 .0346 .0353 .0359 .0366 .0372 .0379 0386 .0392 .0399 .0405 6.6 
12 .0412 .0419 .0462 .0433 .0440 .0448 0455 .0462 .0469 .0476 7.1 
13.0483 .0491 .0498 .0506 .0514 .0522 0529 .0537 .0545 .0552 7.7 
14 .0560 .0568 .0576 .0585 .0593 .0601 0609 .0617 .0626 .0634 8.2 
15 .0642 .0651 .0659 .0668 .0676 .0685 0694 .0702 .0711 .0719 38.6 
16 .0728 .0737 .0746 .0755 .0764 .0773 0783 .0792 .0801 .0810 9.1 
17 .0819 .0828 .0838 .0847 .0857 .0866 0876 .0886 .0895 .0904 9.5 
18 .0914 .0924 .0934 .0944 .0954 .0964 0974 .0984 .0994 .1004 10.0 
19 .1014 .1024 .1035 .1046 .1056 .1066 1077 .1088 .1098 .1108 10.5 
20) ELIT 130 14 eLLO2 16S L174 1185 .1196 .1207 .1218 11.6 
21 .1229 .1240 .1252 .1264 .1275 .1286 1298 .1310 .1321 .1332 1408 
22 .1344 .1356 .1368 .1380 .1392 .1404 1416 .1428 .1440 .1452 12.0 
23.1464 .1476 .1489 .1502 .1514 .1526 1539 .1552 .1564 .1576 12.5 
24 .1589 .1602 .1615 .1628 .1641 .1654 1667 .1680 .1693 .1706 13.0 


25) .1719 .1732 .1746 .1760 1773 1786 =.1800 .1814 .1827 .1840 13.5 


.26 .1854 .1868 .1882 .1869 .1910 .1924 -1938 .1952 .1966 .1980 14.0 
-27  .1994 .2008 .2023 .2038 .2052 .2066 -2081 .2095 .2110 .2124 14.5 
-28 2139 .2154 .2169 .2184 .2199 .2214 .2229 .2244 .2259 .2274 15.0 
29.2289 .2304 .2320 .2335 .2351 .2366 .2381 .2397 .2412 .2428 15.4. 
-80 .2443 .2459 .2475 .2490 .2506 .2522 -2538 2554 .2569 .2585 15.8. 


-381 .2601 .2617 .2633 .2650 .2666 .2682  .2698 .2714 .2731 .2747 16.2 
382.2763 .2780 .2796 .2813 .2829 .2846  .2863 .2879 .2896 .2912 16.6 
-383 .2929 .2946 .2963 .2980 .2997 .3014 .3031 .3048 .3065 .3082 17.0 
-34 3099 .3116 .3134 .3151 .3169 .3186  .3203 .3221 .3238 .3256 17.4 
.85  .3273 .3291 .3309 .3326 .3344 .3362  .3380 .3398 .3415 .3432 17.8 


386.3451 .3469 .3487 .3506 .3524 .3542  .3560 .3578 .3597 .3615 18.2 
37 = .3633 .3652 .3670 .3689 .3707 .3725 .3745 .3763 .3782 .3800 18.6 
-38 .3819 .3838 .3857 .3876 .3895 .3914 .3933 .3952 .3971 .3990 19.0 
-389 .4009 .4028 .4048 .4067 .4087 .4106 .4125 .4145 .4164 .4184 19.4 
-40 .4203 .4223 .4243 .4262 .4282 .4302  .4322 .4342 .4361 .4381 19.8 


-41 .4401 .4421 14441 14462 .4482 .4502 .4522 .4542 .4563 .4583 20.2 
42 .4603 .4624 .4644 .4665 .4685 .4706 .4727 .4747 .4768 .4788 20.6 
-43° .4809 .4830 .4851 .4872 .4893 .4914 .4935 .4956 .4977 .4998 21.0 
44 .5019 .5040 .5062 .5083 .5105 .5126 .5147 .5169 .5190 .5212 21.4 
45° .5233  .5255 .5277 .5298 .5320 .5342 .5364 .5386 .5407 .5429 21.8 


46.5451 5473 .5495 .5517 .5539 .5561 .5584 .5606 .5628 .5650 22.1 
47 5672 .5694 .5717 .5739 .5762 .5784 .5806 .5829 .5851 .5874 22.4 
-6100 22.7 
-6262 .6285 .6308 .6331 23.1 
-6448 .6471 .6494 .6518 .6541 .6565 23.4 


.6588 
.6825 
- 7064 
-7306 
21551 


AUTRE 
.8050 
-8304 
-8560 
-8818 


.9079 
.9342 
-9608 
- 9876 
1.0147 


1.0420 
1.0694 
1.0969 
1.1246 
1.1524 


1.1804 
1.2085 
1.2368 
1.2653 
1.2939 


1.3226 
1.3514 
1.3802 
1.4091 
1.4380 


1.4670 
1.4960 
1.5250 
1.5540 
1.5830 


1.6120 
1.6410 
1.6699 
1.6988 
1.7276 


1.7561 
1.7844 
1.8125 
1.8403 
1.8678 


1.8950 
1.9219 
1.9484 
1.9744 
2.0000 


1 2 
-6612 .6635 
-6849 .6873 
-7088 .7112 
.1330 7355 
1576 7601 
-7824 .7849 
-8075 .8101 
-8330 .8355 
.8586 .8612 
-8844 .8870 
-9105 .9132 
-9369 .9395 
-9635 .9662 
-9903 .9930 

1.0174 1.0202 
1.0447 1.0475 
1.0721 1.0749 
1.0997 1.1024 
1.1274 1.1302 
1.1552 1.1580 
1.1832 1.1860 


1.2113 1.2142 
1.2396 1.2425 
1.2682 1.2710 
1.2968 1.2996 


1.3255 1.3284 
1.3543 1.3572 
1.3831 1.3860 
1.4120 1.4149 
1.4409 1.4438 


1.4699 1.4728 
1.4989 1.5018 
1.5279 1.5308 
1.5569 1.5598 
1.5858 1.5880 


1.6149 1.6178 
1.6439 1.6468 
1.6728 1.6757 
1.7017 1.7046 
1.7304 1.7333 


1.7589 1.7618 
1.7872 1.7900 
1.8153 1.8181 
1.8430 1.8458 
1.8705 1.8732 


1.8977 1.9004 
1.9246 1.9272 
1.9510 1.9536 
1.9770 1.9795 


-6659 
.6897 
1137 
-7380 
~7625 


-7874 
.8126 
.8381 
.8637 
.8896 


.9158 
-9422 
- 9688 
-9957 
1.0229 


1.0502 
1.0777 
1.1052 
1.1329 
1.1608 


1.1888 
1.2170 
1.2453 
1.2739 
1.3025 


1.3312 
1.3600 
1.3889 
1.4178 
1.4467 


1.4757 
1.5047 
1.5337 
1.5627 
1.5917 


1.6207 
1.6497 
1.6786 
1.7074 
1.7362 


1.7650 
1.7928 
1.8208 
1.8485 
1.8760 


1.9031 
1.9298 
1.9562 
1.9821 


STEVENS 


-6683 
.6921 
Sey 
- 7404 
-7650 


-7899 
-8152 
.8406 
.8663 
.8922 


.9184 
-9448 
-9715 
-9984 

0256 


Ue 
1.0530 
1.0804 
1.1080 
1.1357 
1.1636 
ile 
i 


1916 
2198 
1.2482 
1.2767 
1.3054 


1.3341 
1.3629 
1.3918 
1.4207 
1.4496 


1.4786 
1.5076 
1.5366 
1.5656 
1.5946 


1.6236 
1.6526 
1.6815 
1.7103 
1.7390 


1.7674 
1.7956 
1.8236 
1.8513 
1.8787 


1.9058 
1.9325 
1.9588 
1.9846 


-6706 
-6944 
-7185 
-7428 
«7675 


.71924 
8177 
.8432 
.8689 
.8948 


.9210 
.9475 
-9742 
1.0012 
1.0284 


1.0557 
1.0832 
1.1108 
1.1385 
1.1664 


1.1944 
1.2226 
1.2510 
1.2796 
1.3082 


1.3370 
1.3658 
1.3946 
1.4236 
1.4525 


1.4815 
1.5105 
1.5395 
1.5685 
1.5975 


1.6265 
1.6554 
1.6843 
1.7132 
1.7418 


1.7702 
1.7984 
1.8264 
1.8540 
1.8814 


1.9084 
1.9352 
1.9614 
1.9872 


-6754 
.6992 
1233 
1478 
1725 


~1975 
.8228 
.8483 
.8741 
.9001 


.9263 
-9528 
.9796 
1.0066 
1.0338 


1.0612 
1.0886 
1.1163 
1.1441 
1.1720 


1.2001 
1.2283 
1.2568 
1.2853 
1.3140 


1.3428 
1.3716 
1.4004 
1.4293 
1.4583 


.6778 
-7016 
.7258 
.7502 
. T7749 


.8000 
.8253 
.8509 
-8766 
.9027 


.9289 
-9555 
-9822 
1.0093 
1.0365 


1.0639 
1.0914 
1119 
1.1468 
1.1748 


1,2029 
1.2311 
1.2596 
1.2882 
1.3169 


1.3456 
1.3744 
1.4033 
1.4322 
1.4612 


1.4873 1.4902 
1.5163 1.5192 
1.5453 1.5482 
1.5743 1.5772 
1.6033 1.6062 


1.6323 1.6352 
1.6612 1.6641 
1.6901 1.6930 
1.7190 1.7218 
1.7475 1.7504 


TAYE SES UTRSHTS 
1.8041 1.8069 
1.8320 1.8347 
1.8596 1.8623 
1.8868 1.8896 


1.9138 1.9165 
1.9404 1.9431 
1.9667 1.9692 
1.9923 1.9949 


-6801 
. 7040 
- 71282 
.7526 
~TTT4 


-8025 
.8279 
.8534 
.8792 
-9053 


-9316 
-9581 
.9849 
1.0120 
1.0393 


1.0667 
1.0941 
1.1218 
1.1496 
1.1776 


1.2057 
1.2340 
1.2624 
1.2910 
1.3197 


1.3485 
1.3773 
1.4062 
1.4351 
1.4641 


1.4931 
1.5221 
1.5511 
1.5801 
1.6091 


1.6381 
1.6670 
1.6959 
1.7247 
1.7532 


1.7816 
1.8097 
1.8375 
1.8650 
1.8923 


1.9192 
1.9458 
1.9718 
1.9974 
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APPENDIX D Seamed ara 


Transformation of the basic hydraulic equation (2) to the « liptic 
(3) in the text; (20) in this appendix. 


Es e 
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| Equation (2) may be transformed into a form 


1455-21 


Transformation to Elliptical Integrals 


of complete elliptic integral 


ad a solution found to almost any degree of precision desired. Tables are 
yailable of these elliptic functions to 8 or more decimal places. 

Following are a few fundamental relationships that are essential to this 
sansformation. The elliptic integral of the first kind is 


F(G,k)=u=f£_o8 Lk <1] , [x= sing] 
4A - 2 sin 6 
ae ie dx 
ee Dwight 750 
SUE = pote Sl) 
is the amplitude and k the modulus 751.1 
samu 751.2 
ing = snu=x 751.3 
9s ¢ = cnu= Vi -x2 751.4 
1u= V1 - k2x2 751.5 
m= u=- amu 153.1 am 0=0 754.1 
ie U = -snu 753.2 sn 0=0 754.2 
- U = ch u 753.3 Crim ts 754.3 
wm u=dnu 753.4 dn 0=1 754.4 
2u+cn2u=1 755.1 
2 u +k sn2u = 1 755.2 
u - k2cn2u = 1 - k2 751.7 and 755.3 
sn x = cnx dnx 768.1 
< 
he elliptic integral of the second kind is 
E(¢,k) » [eit a3 
[4 - x2 x2 dx [x = sin ¢ | 
Sag lin Ke" Dwight 771 


If the integration is taken between the limits 0 and 7 /2 we obtain complete 


liptic integrals thus 


™% d¢g 
K = | V1 = k2sin2g 


Dwight 773.1 
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and 
E : | Vi-k*sin?¢ ag Dwight 774. 
Let 
y = Hsn2x =D k? sn” x (3A 
H = Dk2 (4 
The first derivative of (3) is 
dy = 2Hsnxd snx (5 
dx 
and from (4) and Dwight 768.1 
dy = 2D k* snx cnx dn x dx (€ 
from(3) and 755.1 
H - y = H(1 - sn2x) = D k*cn2x (7 
from (4) and 755.2 
D -y = D(1 - ksn2x) = D dn2x ( 


From (3), (4), (6), (7) and (8) 


|(D -y) -y) 2 dy 
‘ 


i 
wate x ere x D ke cnx | (2Dk sn x cn x dn x) dx 


= [4D5x8sn4 x eee x dn+ Zid 


( 
(3) may be written y = Hsn2(sin-1 ¢ ). Hence wheny = 0, sin-! 6=x=0 
also late H, from Dwight 751.3 since y/H = 1 


=sn -K the complete elliptic integral of the first kind. 
Adding the limits of integration and substituting (9), (2) becomes 


K 
Q,=4V2g p3/2 4 f er Mes (1 
: ° 
But by Dwight 755.1 and 755.2 


en2x = 1 - sn2x and dn2x - 1 - k2sn2x 
substituting in (10) 


a: 4 Vag v®/2 x4 1 (1 


where 
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: 
a K 
L= ee eee 
45 hie: “0 0 (12) 


rk 


2 | RAP bs 


: 

: ; (6) 

| Reverting to Pierce‘*’ we have 

| (m +1) | sn™ x dx = (m+ 2)(1+k*) {sn ™t2x ax 
) 


m+1 


x dx + sn 


-(m+ 3)k* | x 4 


xcnxdnx, Pierce 567 


which for m = 2, and adding the limits of integration becomes 

Pies - rh 7K 
) 3, sn™x dx = 4(1+k*) jsntax - 5k2 jsnOx dx + sn3x cn x dn x 
fe) 


“pn Ne 


v 


since by Dwight 751.3 and 754.2 sn 0 - 0 and cn K =x = 0 theintegrated terms 
(the last one) of Pierce 567 disappears and there remains 


UK /k -k 
5k*/ sn®x dx = 4(1+k?) | sin*x dx -3 J on?x dx (13) 
+ fe) () 
also for m = 0 in Pierce 567 

K K 
| fidxp= c(h ics) lomexids —93k [sntx ay 
| : (14) 


/ 
° 


Since the left hand term = K we obtain by Pierce 564 
K 
Ve Oe | 
J en?x = <p (K- B) (15) 


° 


since 
E(o, k) = 0 and E(am k, H) = E 


substitute (15) in (14) 
K 


3x2 [ subs dx = 2( 14 12 ) (X22 )i-k (16) 
2 Ke 
or 
K 
J sntax = 2(1 + k2) (K-E)-K (17) 
oO 3k4 3k 


substitute (15) in (17) and (13) 


K 
512 | sno dx = 4(14k*) -[ 2(¥4k*) -(K - E)- kK 
0 3k? Ike 


= 8(1+k2)4 (K - E) - 4(1+k?) K - 3 (K- E) 
3k4 3ke k2 


K 


2 
: 12 |enbx ax = Balin ke )oas <3 K-E| -4 1+k?K 
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: 


1455-24 HY 6 December, 1957 


substitute (15), (17) and (8) in (12) and obtain : 


2 
171 (K-£) -(1+k?) wee (K-E) -,K, 
K 


collecting E - terms 


ra 


B.(,- 156? st 10+ 20k” eh 16ke 8k oF 19kA) 


Sg! 


ae 2(1-k?+k4)E 


collecting K - terms 


K  (15k* - 10 - 20k* - 10k + 5k2 + 5k44 84 16k2 
15 k4 


+ 8k* - 9k2 444 4k2) = 1 x 


(2 - 3k“ + k*) 


15 k# 
combining E - and K - terms 


Ys Ste ara 2) wie ee K | 


19 
15k4 (19) 
substitute (19) in (11) and obtain finally 
= i%2 yep DS/4-- Piisakei+ hE? eS) 
Siig. g [ ( ) K] (20) 


a rigid formula for the theoretical flow through a sharp edged circular weir 
with pee ais contractions and no appreciable velocity of approach. 
In (20) k2 = H/D the head on the invert of the circle in terms of the diamet 
E and K may be taken from tables (see references). 
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PROPORTIONAL WEIRS FOR SEDIMENTATION TANKS? 


. Closure by J. C. Stevens 


| 
| J.C. Stevens,! M. ASCE.—The author is grateful indeed to Prof. Guido 
Di Ricco of the University of Rome, for his contribution to the subject of 
proportional weirs. 
| He is quite correct in assuming that the author was entirely unaware of 
ihe very extensive work done by him in this field. This must be charged to 
ihe handicap of the language barrier that separates us. To illustrate: when 
Prof. Ricco’s discussion appeared any attempt to obtain copies of his papers 
nly met with frustration. From his footnotes, they could not be found in the 
ngineering Societies Library in New York nor any other library where in- 
quiries were made. They were finally obtained through ASCE for whom Prof. 
Ricco had photostat copies made of some of the papers cited, which were 
sassed on to the writer. Two translations of the principal paper were ob- 
ained from two local Italian sources, but both had to be retranslated by the 
wuthor into the King’s English no doubt because neither translator know hy- 
lraulic terms — more frustrations. 

The credit for inventing the proportional weir must undoubtedly be given 
o the late Dr. Oscar Van Pelt Stout, M. ASCE, while he was Professor of 
sivil Engineering at the University of Nebraska. He designed, tested and put 
0 practical use such a weir in 1896 and presented a paper on it at the first 
neeting of the Nebraska Engineering Society held in Lincoln in 1897. 

These facts were submitted to the author by Prof. Steponas Kolupaila of 
Votre Dame University, who also submitted a bibliography on the subject 
vhich, with his permission, is being appended hereto. It is always best to 
‘ive full credit where credit is due. 

The writer was a student under Prof. Stout. He faintly remembers being 
old about the Stout proportional weir while in college, but did not know or 
ad forgotten that its mathematical details had been made public some 8 
ears prior. During the intervening 50 years it was easy to forget until Prof. 
<clupaila’s most excellent system of references brought the facts cogently to 
nind, for which grateful acknowledgment is hereby made. 

There is a typographical error in the footnote to Prof. Ricco’s discussion 
S printed in No. HY 1, Feb. 1957 p 1177-41. The date of his paper in 
/Energia Elettrica should be October 1936 instead of 1956. 

_ Beginning with the author’s basic equation (7) Prof. Ricco presents his 
eneralized formula(2) in which the lateral width of the weir (x) for any depth 
y) is evolved in his equation( 4) in terms of the gamma function. 


. Proc. Paper 1015, June, 1956, by J. C. Stevens. 
. Cons. Hydro. Engr., Stevens & Thompson Engrs., Portland, Ore. 
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Had the author known of this he might have saved himself much of the 
laborious series developments with which his paper is adorned. 

The use of a series of chords instead of a purely curvilinear profile for 
the weir is also an innovation that may well be adopted in this country. 
Compare Figs 4 and 5 with the author’s Fig 2. It appears that by this expe- 
dient the proportional weir can be designed to produce almost any desired 
pattern of outflow from a tank or reservoir and might, with profit, be applied 
to some cases of river regulation. 

It is hoped that this paper and the discussion will stimulate a more exten- 
sive use of the proportional weir than exists at present in this country. It has 
many advantages and could well be used more extensively. 
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WHE PROBLEM OF RESERVOIR CAPACITY FOR LONG-TERM STORAGE2 


Closure by Aly Fathy and Aly Salem Shukry 


| ALY FATHY! and ALY SALEM SHUKRY2.—The authors find that the dis- 
jussions of their paper have brought forth no consequential information on the 
jubject under study with the exception of Mr. Hanna’s exposition on the ques- 
ion of sub-annual storage. That exposition fills a gap in the inquiry that was 
assed over rather lightly. 
| In Mr. Mostafa’s discussion an attempt was made to demonstrate the in- 
\dequacy of a normal capacity of 45 milld. m3, as arrived at through applica- 
ion of the authors’ method to the projected long-term reservoir on the Nile 
it Aswan. In that attempt, however, sub-annual effects were introduced. In 
he light of Mr. Hanna’s observations it will be seen that when these effects 
ire taken into consideration the extra capacity provided for flood relief should 
included in the working capacity. According to Mr. Mostafa’s own calcula- 
ions the total capacity of 75 milliards (45 + 30) would be more than sufficient. 
The discussion by Messrs. Simaika and Boulos is devoted mainly to 
emonstrating the unreliability of conditions derived from past observations 
§ a guide in forecasting the future behaviour of a natural event. In a field of 
Nquiry where probability considerations play the leading part such an attempt 
s really superfluous. As far as hydrological phenomena are concerned, the 
rincipal criteria to which reference may be made are: 


1. The Mean. 

2. The Standard Deviation. 

3. The Maximum Deviation of a single observation from the Mean. 

4. The Maximum Deviations of the arithmetic mean for groups of observa- 
tions of varying number. 


Of these criteria, the Standard Deviation is the most erratic and the least 
menable to provisional adjustment. Thus a method of treatment which de- 
ends primarily on that parameter would be the least reliable. 

In Dr. Hurst’s discussion a slip has been made in the attempt to recon- 
trruct the expression of the ideal range (Eq. b). As the length of the critical 
sriod is supposed = N/2, the derived expression should be multiplied by V2 : 
hat would raise the numerical factor to the same value as given by the 
ithors in Eq. 15. 

Dr. Hurst considers it faulty to adopt the mean value N/2 as a standard for 
e length of the critical period. It would certainly be more legitimate to do 
) with regard to a parameter that fluctuates about a definite mean than to 


. Proc. Paper 1082, October, 1956, by Aly Fathy and Aly Salem Shukry. 
. Part-time Prof., Faculty of Eng., Alexandria Univ., Alexandria, Egypt. 
_ Asst. Prof., Faculty of Eng., Alexandria Univ., Alexandria, Egypt. 
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one that varies progressively such as the range R. According to Dr. Hurst’s 
“Century Storage” theory a standard value is assigned to R corresponding to 
a working period of 100 years. It is here that disregard of the variability of 
R may be questioned. 

Incidentally, it may be noted that by assigning a standard value of 1/2 to 
the relative length of the critical period the resultant expression of R is 
brought to the same form (though not the same value) as arrived at by Dr. 
Hurst himself through a very elaborate mathematical analysis. Mr. Hanna’s 
observations on that point will be found most illuminating. 

The authors wish to stress the fact, however, that in the practical applica- 
tion of their method of treatment no reference at all is made to the ideal 
range which is an indeterminate quantity. The capacity S needed to guarantee 
a draft less than the mean (as would always be the case) is derived straight- 
forward from the deviation-curve. The tying up of S with R, as proposed by 
Dr. Hurst, greatly complicates the problem. Both of his expressions of S 
(Eq. 8 in the original paper and Eq. 27 in Mr. Hanna’s discussion) will be 
found to lead to absurd results when the difference between the guaranteed 
draft and the mean is large. According to (8), S becomes zero when (M - B) = 
1.04 o . That quantity is usually much in excess of the maximum deviation of 
one observation or of the mean for a small group. For the Nile at Aswan, 

M = 93 and o = 19.8. Thus S becomes zero when B = 72.4. That figure is 
30.4 milliards greater than the observed minimum supply of one year, which 
is 42. Evidently, the storage provided should cover the deficit in one year at 
least. 

The logarithmic relation (27) yields a value of S = 4.9 for a guaranteed 
draft of 60 millds./ann. Actually, the storage needed to cover the deficit in 
the worst year would be 18 milliards. 

The above examples bear sufficient evidence of the need for a totally dif- _ 
ferent line of approach. It is the authors’ belief that the line they have sug- 
gested is the only alternative. 


Corrections: Fig. 7 of the original paper did not contain the construction 
referred to in the text. The correct figure is reproduced herewith. 
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TECHNICAL PROBLEMS OF FLOOD INSURANCE2 


Closure by H. Alden Foster 


: 
| 


: 

) H. ALDEN FOSTER,* M. ASCE.—The discussion by Mr. Alexander em- 
phasizes some of the difficulties encountered in preparing probability esti- 
mates for floods. From a practical viewpoint, it is impossible to prepare an 
absolutely reliable probability curve from a limited flood record. The record 
itself may be represented by a “standard” flood-frequency curve; but the cor- 
responding curve for any future “sample” of equal length may be appreciably 
different within so-called “tolerance limits” or “confidence limits.” 

| It has been pointed out by numerous writers that the probability method is 
not suitable for the determination of the flood magnitude for design of the 
spillway of an important dam. As pointed out by the writer, however, in the 
present paper and on other occasions, (1) probability methods have a proper 
and useful application in certain fields, particularly in estimating average an- 
nual flood benefits for a flood-control project or in the design of cofferdams 
for sites subject to floods, where precision of the estimate is not of primary 
importance. For satisfactory determination of flood-insurance rates, some 
kind of probability estimate is essential. In this case, however, the possible 
range of the probability estimate between the confidence limits may add ap- 
preciably to the risk assumed by the insuring agency. 

Mr. Alexander points out the differences of opinion among mathematicians 
as to the proper method of estimating the confidence limits for a particular 
record which is essentially a statistical problem. A fellow-countryman of 
Mr. Alexander(2) has discussed methods for determining such confidence 
imits for the log-normal and Pearson Type III distributions. 

The writer agrees with Mr. Alexander that the method of estimating floods 
Ny use of the “maximum possible precipitation” requires clarification, par- 
icularly as to its relation to probability estimates. It is understood that a 
laper by Mr. Alexander discussing this question, together with the interesting 
subject of the relation of probability estimates to “factors of safety” for dams, 
vill soon be published by the Institution of Civil Engineers of Australia. 

Professor Kolupaila’s remarks about the writer’s 1924 paper are greatly 
ippreciated. The problem of “sampling errors” and confidence limits was 
jot well understood by most engineers and hydrologists 30 years ago. An ap- 
reciation of the importance of these problems would doubtless have pre- 
ented unsuitable application of the theoretical probability methods in hydro- 


Ogical studies. 


. Proc. Paper 1165, February, 1957, by H. Alden Foster. 
*Prin. Associate, Parsons, Brinckerhoff, Hall & Macdonald, New York, 
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A HIGH-HEAD CAVITATION TEST STAND FOR 
HYDRAULIC TURBINES#® 


Discussion by Robert T. Knapp 


| ROBERT T. KNAPP, 1 M. ASCE.—The authors are to be congratulated upon 

e design and construction of a fine cavitation laboratory for the testing of 
tydraulic turbines. Their careful planning should pay dividends for a long 
lime. The design of the dynamometer and the method of mounting the test 

achine are particularly noteworthy. There are many other nice things that 
san be said about this new laboratory. However, in a discussion it sometimes 
seems more fitting to raise differences in viewpoint. 
| One of the possibilities sacrificed in the unique turbine mounting is that of 
risual observation of cavitation in the wicket gates or in any part of.the instal- 
| ion except the draft tube and the discharge edges of the runner vane. This 
may be serious in the use of the laboratory as a development tool. Another 
infortunate limitation, in the writer’s opinion, is that the test machine cannot 
”@ Operated as a pump and that transient conditions cannot be investigated. 
This seems unfortunate in light of the growing interest in pump-turbines. 

_ The writer is surprised at the weight put on the simulation of field test 
echniques. Such techniques are often dictated by necessity rather than by 
lesirability. It should not be difficult, with properly selected instruments, to 
ret at least equal accuracy and reliability with one operator, with or without 
in assistant, and possibly require less time to make the same run. In this 
‘onnection, it would be interesting to know whether the laboratory has a 
llanned program for the regular recalibration of the venturi meters by the 
jame agencies using the same test installations that were employed in the 
riginal calibration. Such a program would appear necessary if the quantity 
neasurements are to be maintained at the same high standard of accuracy as 
hat obtainable with the other test measurements. 

The writer does not feel that it is permissible to adjust the air content 
olely on the basis of preventing the accumulation of air in the manometer 
onnections. Unfortunately, there is as yet not enough known about the rela- 
ionship between air content and cavitation performance to establish a better 
riterion. It is hoped that the laboratory operators will bear this possible 
elationship in mind so as not to become too enamoured of their current 
>chnique. 

The samples of test results presented in Figs. 4 and 5 are very interesting. 
‘he authors are to be complimented on not using a “step-up” factor on their 
1easured efficiencies. At the Reynolds numbers at which these tests were 
1ade a step-up would be questionable. The paper seems to imply that 


‘ Proc. Paper 1201, April, 1957, by W. G. Whippen and G. D. Johnson. 
. Prof., Hydr. Eng., California Inst. of Technology, Pasadena, Calif. 
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continuous operation of the turbine at “critical o ” would not produce signifi- 
cant cavitation damage. Further comments on this point would be appreciated 
It would be interesting to see “sigma break” curves for a series of turbines 
of differing designs. The writer is under the impression that it is not uncom- 
mon to obtain appreciably steeper breaks. In this connection it would seem 
desirable to emphasize the authors’ warning that if emergency operation at 
very low sigma values is anticipated, the permissible limits of operation must 
be based on a complete and reliable set of model test curves covering the 
sigma values likely to be encountered. 

The writer agrees with the authors that it is most desirable to test the ef- 
ficiency and cavitation performance of hydraulic machines at prototype heads 
and velocities. He made the same choice in 1933-34 during the design of the 
Hydraulic Machinery Laboratory of the California Institute of Technology and 
has yet to regret it. However, it must be remembered that this approach 
ignores several facets of the cavitation process. For example, the vertical 
dimensions of large machines are great enough to affect the extent of the 
areas subjected to cavitation. To simulate these effects model tests at proto- 
type Froude numbers are necessary. But experiments show that such tests 
also fail to agree with the cavitation performance of the prototype, primarily 
because they fail to alter the properties of the test liquid to correspond to the 
changed scale. It is feared that this dilemma will continue to confront not 
only the authors, but also all other workers in cavitation laboratories for 
some time to come. 

In their concluding paragraphs the authors propose that accurate model 
tests, carried out in adequate laboratories such as the one they describe, can 
be substituted for the more expensive and the consuming field tests. In 
general, the writer agrees with this proposal. He would like to emphasize, 
however, that such model tests to be acceptable must include all of the hy- 
draulically significant features of the field installation. Furthermore, the 
model turbine must be a true model of the machine as installed in the field. 
In many cases this would require the construction of a new model turbine, 
since the model used in the design study may not be geometrically similar 
to the turbine as installed. Two common deviations from similarity are: 

(1) design changes introduced between the model test and the prototype con- 
struction, and (2) errors introduced in the manufacture and installation of the 
full-sized machine that cause it to differ from the final design drawings. In 
the opinion of the writer, if model acceptance tests are to replace field tests, 
the laboratory responsible for making these tests should also have the 
responsibility of ascertaining and certifying to the geometric similarity of 
the model and field installations, as well as to the laboratory accuracy of the 
test measurements. 
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THE EFFICACY OF FLOOR SILLS UNDER DROWNED 
HYDRAULIC JUMPS? 


Discussion by Alvin J. Peterka, 
Alfred C. Ingersoll, and S. Leliavsky 


) 
| 
| 


ALVIN J. PETERKA,! M. ASCE.—The conclusions drawn by Mr. Shukry 
jare contrary to general practice procedures which have been developed over 
ja period of years from hydraulic laboratory tests, prototype tests, and other 
jprototype experiences. From a limited number of velocity measurements on 
ja small hydraulic model, Mr. Shukry has concluded that a sill placed near the 
imiddle of the apron is in the optimum position. His conclusion is based on an 
junproved assumption that “The efficiency of any sill against bed erosion is 
indicated by the rate of adjustment of the flow to the normal distribution in 
jthe downstream channel.” In addition, the tests were conducted on a fixed- 
lbed model where erosion could not be measured, and the conclusions are con- 
trary to test results obtained from a movable bed model of the same struc- 
ture, Figure 2. 

It is not clearly stated in the paper whether the difficulty in reconciling 
erosion test results was between several models of the same structure or 
between the models and the prototypes. The “erratic results” of the proto- 
type scour tests are not described. It is not clear, therefore, why Mr. 
Shukry abandoned the conventional movable bed tests and used velocity 
traverses to indicate scour tendencies. It is also not clear why Mr. Shukry 
has mathematically defined a scour hole since he used a fixed-bed model in 
his tests. However, certain impressions are obtained in reading this paper, 
and it is on these impressions the writer must base his discussion. 

The fact that agreement between model and prototype erosion tests was not 
found may be due to any one of several factors. It should be remembered that 
models are mechanical devices and that the interpretation of the model test 
results is often as important as the model itself. A 3-bay model will repre- 
sent a 46-bay prototype spillway only if all factors concerning operating pro- 
cedures are considered. By proper selection of model scale and bed material 
size, it is possible to predict qualitatively where the most severe erosion will 
occur in a prototype structure. It is also possible to construct several or 
more models of the same structure and obtain identical erosion patterns. 
Contrary to the statement in the paper, it is possible, therefore, to obtain 
qualitative similarity between model and prototype and between different 
models of the same structure. On the other hand, it is not generally believed 
that an ordinary movable bed model can be used to determine quantitative 


2. Proc. Paper 1260, June, 1957, by Ahmed Shukry. 
l. Hydr. Engr., Bureau of Reclamation, U. S. Dept. of the Interior, Denver, 
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erosion test results that are directly transferable to the prototype structure. 
Mr. Shukry infers in his paper that this was to be expected. 

If prototype erosion was greater than that predicted by the original sec- 
tional model, unsymmetrical operation of the 46 bays could have produced al- 
most any conceivable scour pattern. Degradation of the prototype riverbed 
could also account for lack of agreement between model and prototype erosion 
patterns. Clear water released from a reservoir tends to pick up a silt load 
as soon as possible. The riverbed below a dam is sometimes attacked in this 
manner when the bed material is very fine. This action has nothing in com- 
mon with the erosive action encountered downstream from stilling basins. 
Prototype riverbeds have been lowered several feet in a few years as a result 
of degradation. 

Mr. Shukry states that for the Edfina barrage and others that the best loca. 
tion for the sill is in the middle of the floor. This statement is meaningless 
unless the length of the solid floor is defined in terms of the flow depth or in 
some dimensionless number. Sufficient data are not presented to allow the 
reader to compute the probable length of apron required, but it appears from 
Figure 1 that the aprons are longer than those usually required for discharge: 
and heads of this magnitude. For example, it appears that a floor extending t 
Section IV or V, Figure la, would provide a sufficiently long paved floor. A 
sill at the end of the apron and riprap downstream from the sill would ordi- 
narily provide ample protection to the riverbed, since the total energy in the 
flow is relatively small. 

A sill is most effective in preventing scour when it is located so that mild 
bottom velocities strike the sill and are directed away from the riverbed in a 
positive and gentle manner. Near-zero velocity then prevails at the bed and 
water interface and optimum conditions for the prevention of scour are 
present. With the sill located in the middle of the apron, a redistribution of 
velocity occurs before the flow contacts the riverbed and possibly before a 
boundary layer is established, resulting in greater scour. Broad sloping sill; 
usually produce weaker ground rollers than steep sills, reducing the circula- 
tion and loss of bed material at the end of the apron. A properly shaped 
dentated sill also helps to produce zero velocity at the interface and a negli- 
gible ground roller. Scores of stilling basins with dentated sills have been 
tested in various laboratories and have been found to be superior to solid 
sills. 

It should also be pointed out that the action described on the Edfina apron 
should not be classed as hydraulic jump action. It appears (sufficient data 
are not given) that the Froude number of the incoming flow is very low, per- 
haps no greater than 2. In a hydraulic jump, the energy loss is only 6 percer 
of the total when the Froude number is 2. If the jump is drowned, the loss ~ 
would be even less. The action on the apron in the Edfina tests would probab 
be better described, therefore, as a roller action where energy loss is relate 
to the turbulence developed on the apron. Statements concerning the shape ai 
placement of sills should therefore be applied, at best, to the structures 4 
tested and should not be considered to be generally true for hydraulic jump 
aprons. 


ALFRED C. INGERSOLL,! M. ASCE.—In this paper Mr. Shukry has con- 
tributed greatly to the hudraulic designer’s understanding of the effectivenes 
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tof barrage floor sills at various locations in minimizing downstream scour. 
The need for a clarification of the previously existing multifarious and con- 
tflicting design criteria is admirably demonstrated by the tests on a movable 
jbed model of the old Esna Barrage. The model tests predicted scour be- 
thavior that has been largely refuted by experience on full-scale structures. 
| As sound as the author’s principal conclusion - that anti-scour sills can 
tbe investigated by means of fixed-bed model studies - appears to be, it is 
montheless unfortunate that no case is cited wherein either a movable-bed 
‘model or a prototype structure behaves clearly according to the fixed-bed 
ipredictions. It is true that the fixed-bed studies on the Edfina Barrage seem 
ito bear out the earlier of the two sets of movable-bed studies, but the more 
irecent set of tests, shown in Fig. 2, gives results diametrically opposed to 
ithe conclusions of the fixed-bed studies, with the mid-floor (or edge) position 
jof the sill being the best by one test method and the worst by the other. It 
)would seem that this all-important question cannot be truly resolved until the 
behavior of this and similar prototype structures can be observed. 
_ In Eq. (6) the author appears to have evolved an expression for the shear 
istress that is independent of the viscosity or Reynolds number, and he there- 
fore proceeds to apply the Froude law to the fixed-bed model studies, com- 
aring one model test qualitatively with another, to predict which prototype 
istructure will experience the least scour. The writer feels that there is a 
danger in this, and that it lies in the effect of the Reynolds number upon the 
velocity distribution (especially that downstream) and consequently upon # and 
therefore upon T,. 

Although the analogy is admittedly rough, one may seek a relation between 
£& and TR from information on smooth pipes. Inasmuch as the scale of the 
Edfina Barrage model is 1:66.7, the Reynolds number scale, when tested by 
the Froude law, is 1:(66.7)3/2 or 1:543. As shown by Rouse,(1) the variation 
to be expected in § for this change in Reynolds number is on the order of 
from 1.01 at the prototype scale to 1.03 at the model scale. 

While this magnitude of difference in itself may be negligible, it may be 
observed that it probably applies more to the downstream condition, e.g. 
Station VI, than to the upstream condition at Station I. That is, the velocity 
distribution immediately downstream from the gate is affected most by the 
gate opening and position, the viscous forces having little opportunity to come 
into play in the short distance before Station I. As the flow continues over 
the long floor and the bed-protection works, however, the viscous forces tend 
to establish a more normal velocity distribution, as is shown in Fig. 6. Thus 
when the slightly differing values of B 9 are subtracted (as in Eq. 6) from a 
value of 8 4 which is the same for both model and prototype, the resulting dif- 
ference in T, could be significant. The net effect of this would be for the 
prototype to develop a higher relative bed shear stress than the model and it 
could conceivably occur that fixed-bed model studies of competing prototype 
designs would lead to an erroneous evaluation on this account. 

The curves of Fig. 7 are especially interesting to the student of open- 
channel hydraulics. The rising surface profiles for tranquil flow over a 
horizontal bed are, of course, impossible for normal channel velocity dis- 
tributions. For the marked maldistributions observed, however, the validity 
of these profiles may be confirmed by an energy analysis applied to the . 
reach between Stations I and VII. 
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Considering condition (a), with discharge passed under the gate, the 
velocity distribution curve at Station VII may be analyzed to yield a value 
for the energy correction factor, 


SV JP v Faz ay 


Vv. ? bY (1) 
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of 1.32. With a depth of 4.50 meters and average velocity of 2.13 ms. per sec 
the specific energy, 


. (2) 


is 4.81 ms. 

Allowing for a reasonable loss of, say, 5 cm between Stations I and VI, 
the specific energy at Station I must be about 4.86 ms. From Fig. 7(a) the 
depth at I may be scaled at 0.68 x 4.50 or 3.06 ms. The kinetic energy must 
then be 1.80 ms. Inasmuch as the average velocity at this section is 3.14 ms. 
per sec., the velocity head is only 0.501 ms. and clearly this will require a 
value of a of 1.80/0.501 or 3.59. 

Turning to the velocity distribution diagram for Station I of Fig. 6(a), it is 
possible to compute a value for a by scaling the diagram, provided it be 
assumed that this diagram is representative of the velocity distribution 
across the entire 10.5-meter width of one vent and pier. The value of a so 
obtained is approximately 2.92, not sufficient to account for the high kinetic 
energy of 1.80 meters. 

The explanation for this discrepancy seems to lie in the lateral distribu- 
tion of velocity as it is badly distorted at the downstream edge of the 2.5- 
meter wide pier, just at Station I. If, then one assumes that the distribution 
diagram shown applies only to the 8-meter width of the channel upstream, 
while the average velocity, V, applies to the full 10.5-meter width down- 
stream, the value of a becomes 4.85. This is indeed more than enough to 
account for the kinetic energy. The difference between this and the minimum 
requirement for @ of 3.59 isreadily absorbed in considering lateral 
velocity distributions between the two extreme cases postulated. 

This paper is extremely valuable in opening new vistas of research, where 
in the elusive scour phenomenon can be studied with the fixed-bed model and 
the pitot tube. It is to be hoped that the author will continue these significant 
studies and that he will, whenever possible, relate the predictions so ob- 
tained to observations of scour on prototype structures. 


REFERENCES 


1. “Engineering Hydraulics,” edited by Hunter Rouse, John Wiley & Sons, Inc 
New York, N. Y. 1950 Fig. 9, p. 401. 


S. LELIAVSKY*.—For over twenty years the writer held the post of 
Director, Designing Service, in the Department in charge of the construction 


*Civ. and Hydr. Eng., formerly in the employ of the Egyptian Government. 


Asce DISCUSSION 1456-19 


ind remodelling of the barrages shown in the authors Fig. 1, and in this 
¢apacity, he ventures to make some constructive remarks about Professor 
thukry’s interesting and timely paper. 

) The authors’ argument is roughly this: since it is impossible to reproduce 
m the model the true ratio of the structural dimensions of the prototype to 

he size of the particles of the granular soil, let us drop this material al- 
sogether, and carry out the model tests within rigid boundaries, deriving our 
fonclusions from the observed behaviour of the flowing water. 

This is not a new approach, but it was usually believed to be particularly 
suitable for simpler structures, such for instance, as the Gully Control works, 
which were the subject of the experiments of B. T. Morris and D. C. 
fohnson(1); but whether it is also applicable to works of the size and im- 
vortance of the Nile Barrages—on the safety of each of which may easily 
lepend the welfare of as much as a million rural population—is a matter re- 
juiring much thought, and consideration of the various aspects of the relevant 
problems. 

_ Had the author been personally responsible for such designs, he would 
wave certainly realised the complexity of these problems and might have 

en, very probably, abstained from veiled complaints against the alleged 

onotony of the adopted design solutions. (Which remark should not be inter- 
yreted as adverse criticism, but only as an attempt to clarify the writers 
.rgument.) 
_ Let us begin with the Assiut Barrage (d in Fig. 1) the remodelling of which 
vas the first among the large Nile Works, to be designed according to model 
ests completed before the drawings were prepared and actual work began. 
[The author has envisaged, in his experiments, three alternatives—flow above, 
n between and below the gates—but he has not considered the fourth alterna- 
ive, i.e., flow above the undergate weir; and yet, such a weir formed (accord- 
ng to the model tests referred to) a significant part of the anti-scour-devices 
rrangement, adopted for the remodelling of this work. 

Originally, the object of this weir was a matter of gate design. In fact, the 
ew regulation order, for which the Barrage was then being remodelled, 
alled for a greater total height of obstruction than the old one, but the per- 
nissible height for the new gates was limited by considerations pertinent to 
he housing space available behind the parapet. Hence the weir. 

But when the consulting engineers, Messrs. Coode, Wilson, Mitchell and 
‘faughan Lee, suggested metal trestles above the parapet level and in spite 
f strong opposition the suggestion was adopted (see Fig. 2),(2) the housing 
pace problem was ipso facto solved, and the weir was therefore capable of 
eing dispensed with. 

This should have been, apparently, the proper solution from the hydraulic 
tandpoint, since it yielded a more regular pattern of flow lines, but a 
eneral analysis of all the then available evidence, derived from numerous 
10del tests with granular material, tended to show the merits of the under- 
ate weir in fighting tail-erosion; and the results in practice, confirmed 
ntirely this conclusion, notwithstanding the fact that it had been derived 
*om experiments with sand, which the author appears to condemn. 

In the pre-modelling period, the Inspector in charge lived, during the 
ood months, in constant terror of what that blessed Barrage would do next, 
r apart from a permanent scour-hole (“bayarah”), local scour in closer 
icinity to the apron was not unfrequent, and had to be immediately levelled 
it with rubble-stone dipped from barges, of which more than 130,000 cub.met. 
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had been thrown into the water, since the barrage had been built. Little 
wonder, therefore, that when we saw that after the remodelling, the surface 
of the earthen channel behind the barrage remained almost scourless, we be- 
ean believing in miracles. 

Of course, model tests with sand are not to be interpreted too literally, but 
with a certain experience in the subject, you may very well distinguish in 
*xamining the information they supply, the basic and substantial, from the 
ccidental and inaccurate. In fact, inconsistencies and inaccuracies in such 
»xperiments when they occur are not so much the result of differences in 
rompactness and cohesion and other systematic causes, as the author seems 
jo believe, but are chiefly due to the fact that it is difficult (if not impossible) 
j0 reproduce accurately in a very small size experiment, the true case his- 
jory of the phenomenon. 
| The only way of dealing with such essentially accidental inaccuracies is to 
increase the scale of the tests; which means that for important works, such 
K s the Nile Barrages, the scale of an experiment with sand, must not be less 

an 1/10. This explains that, immediately after the completion of the Assiut 
Sarrage Remodelling, we included in the project of Mohammad Aly Barrages, 
4 generously proportioned test-laboratory, in which such large scale experi- 
ments could easily be carried out. It may be unfortunate that, for administra- 
live reasons, this laboratory was not completed, until after the second world 
war. 
Reverting to the author’s tests, it might be of interest to find out whether 
i. method is sensitive enough to perceive the effect of the undergate weir 
m downstream erosion? Additional experiments to clarify this point appear 
10 be, certainly, desirable. 

The writer wishes to emphasize that he has no intention whatever of ob- 
ecting against the author’s solution, but he believes that the more usual type 
of model-test with sand, must take precedence on it, and that the two systems 
ean then be used in parallel, according to a consistent, self-contained pro- 
rram, aiming at mutual verification of the results. 

The next important work on the Nile, after the Remodelling of Assiut Bar- 
“age was completed, was the construction of the Mohammad Aly Barrages. 
[he hydraulic lay-out of this project, for which the writer was particularly 
-esponsible, incorporated two main requirements: - 


A) To take advantage of the lesson of Assiut Barrage. 

B) To take, also, into account that the backwater curve of the new barrage 
extended beyond Cairo, a large city with over a million population and 
many valuable old buildings resting on very dubious, if any, foundations. 


Thus, as the yearly flood began to subside and the water table started to 
‘eceed, frequent instances occurred in which such buildings fell, causing loss 
f life and property. It was, therefore, essential not to widen the range of 
vater level variations, i.e. reduce the additional heading-up, consequent on 
he construction of the new barrage, to a minimum. 

The success of the solution of this problem deserves, perhaps, being 
laced on record: during the highest flood after construction, the afflux on 
his barrage did not exceed 2 cm (3/4"'). To achieve this result the writer 
‘emoved (or smoothed out) all solid surfaces exposed to the direct action of 
he shock of water, or causing separation of flow from solid boundary. 

Thus, though the Remodelling Project of Assiut Barrage was taken as 
‘a the sill of the hypothetical undergate weir was joined with the 
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downstream and upstream edges of the floor, by straight lines; thus yielding 
two sloping surfaces on each side of the highest point under the gates. 

Better to explain this semi-intuitive solution, attention is called to Fig. p(é 
in which the profiles of the three main Nile Barrages are superposed upon 
one another. In this drawing, the method of the sloping flour of the Mohamma 
Aly Barrages, as used in loco of the weir of Assiut Barrage, is quite obvious, 
and the solution seems to be evident but. . . the writer was, at the time, 
severely criticized for using it, because it was then alleged that sloping floor, 
intensified downstream scour, and the designs for Gebel Aulia Dam, on the 
White Nile, which was then in course of construction, were modified accord- 
ingly. 
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The heresy, in fact, arose from a wrong interpretation of another set of 
model experiments with sand, and since subsequent practice fully vindicated 
the writers views on the matter, this problem is another subject on which 
Prof. Shukry’s experiments could throw more welcome, additional light. In 
fact sloping floors seem not to have been considered by him at all. 

It remains to analyse the conditions prevailing at the barrage of 
Nag-Hammadi (see Fig. 1-e). This is a sui generis case, because the floor — 
of this barrage is built to three different levels (58.50, 59.50 and 60,50 metre 
above the conventional irrigation datum, which at some time was believed tg 
be the mean sea level at Alexandria). 

As regards downstream scour, this floor arrangement is by no means a _ 
negligible contingency; in fact, as shown by both, the writer’s calculations of 
the probable depth of scour (see Fig. c),(4) and also by actual soundings take 
in the channel of the river, downstream of the apron, the depth of scour de-- 
creases gradually with the discharge in the part with the highest floor a BS 
but in the deepest portion of the barrage, it begins to increase again as the — 
river drops to its summer levels. 

This is another aspect of the scour problem for the Nile Barrages, which 
might be covered by Professor Shukry’s valuable tests. 
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But, the main point on which further experiments of the same type appear 
be necessary, is that the scour-hole itself will basically affect the flow 
ittern, and therefore, to be able to draw significant conclusions from Prof. 
ukry’s tests, we must possess, in addition to the results obtained for a 
vel channel-bed, full information obtained in the same manner for scour- 
les of various depths; viz., using available empirical information on the 
ofiles of the most probable scour-holes, the solid bed, which in these ex- 
riments replaces the earthen channel, must be shaped accordingly, and 
en, tests carried out and diagrams prepared as for a level channel for vari- 
S assumed depths of the scour-hole. This will then supply exhaustive in- 
rmation on the matter and may, possibly, yield a convincing result. 

In fact, it must be realised that scouring action may be started at any time 
d for any accidental reasons, and we must be sure that if this occurs, the 
le will be filled in naturally, or will not attain a dangerous depth, before we 
ve had time enough to fill it in with rubble. 
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Broadly speaking, the Nile Barrages must not be visualised as granite 
castles built on a rocky foundation, but rather, as islands made of friable 
material, which are held in precarious equilbrium, afload on an immense 
sea of almost oily, semi-liquid misture of water and very fine sediment. 
Granted that this description is an obvious exaggeration, from the designers 
standpoint the exaggeration is in the right direction. 

Every precaution must, consequently, be taken to ensure stability of 
granular material, and departures from established practice must be care- 
fully considered. 

The writer, therefore, concludes that it is too early to dispense with the 
standard type of model-tests with sand altogether, replacing them by 
Professor Shukry’s interesting method; but this last method must be used to 
interpret, and sometimes explain, the older type of such experiments. 
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STILLING BASIN EXPERIENCES OF THE CORPS OF ENGINEERS? 


Discussion by Alvin J. Peterka 


_ ALVIN J. PETERKA,! M. ASCE.—The paper by Mr. Berryhill is an excel- 
int summary of “how things have worked out.” With a few exceptions, it ap- 
vars that design procedures based on hydraulic model tests have provided 
latively trouble-free operation. 
_ The writer wishes to raise a question, however, concerning the proper 
ppth of tail water, Dg, for a stilling basin containing baffle piers. It has 
sen common practice in some offices to reduce the conjugate depth, Dg, by 
’ to 20 percent because of the restraining effect of the baffle piers. Certain 
ronomies result; the basin floor may be set at a higher elevation reducing 
<eavation costs, and the training walls may be correspondingly reduced in 
sight. In some cases these savings result in a lowered basin cost. In other 
uses, however, the savings, in terms of other costs, are minor. 
In either case, the writer believes that full conjugate depth or more (about 
percent more) should be used when setting the floor elevation, for the fol- 
wing reasons: 


1. The safety of the entire structure may hinge on whether the baffle piers 
smain fully effective after they have been subjected to the wear of abrasion 
id/or the erosion caused by cavitation. Rounded baffle piers are considera- 
y less effective than square-edged piers and only nominal damage may 
nder the baffle piers ineffective. With the baffle piers ineffective, the apron 

too high and the tail water too low to produce the proper jump action. 
angerous riverbed erosion may result. Difficulties typical of a high apron 
‘e illustrated by Lock and Dam No. 1, Missouri River and Bull Shoals Dam. 


2. Due to degradation of the river channel after the dam is constructed, 
e tail water elevation becomes progressively lower over a period of years, 
ducing the factor of safety against jump sweepout and lowering the effi- 
ency of the jump. On the other hand, an apron designed for conjugate depth 
ll operate properly even if the baffle piers are not replaced. 


3. Experience has shown that actual tail water elevations are often lower 
an the computed values, used for design or model testing of the basin. This 
ay have been true for Kanopolis Dam, which was found by the author to be 
erating at 0.7 Dg. Unless the tail water curve is known to be correct, 
tra tail water depth should be provided in the design. 


4. In the temporal sense, the actual tail water depth lags the depth shown 
the tail water curve during periods of increasing discharge. On some 
Proc. Paper 1264, June, 1957, by R. H. Berryhill. 

Hydr. Engr., Bureau of Reclamation, U. S. Dept of the Interior, Denver, 
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rivers, it may take several hours for the tail water to stabilize after the dis- 
charge is set. When it is necessary to open spillway gates rapidly, the basin 
must be capable of performing satisfactorily with tail water elevations below 
the stable values. 


5. A slightly drowned jump, in the writer’s opinion, produces a better all- 
around energy dissipator than one at conjugate depth. The slight loss in effi- 
ciency (energy dissipation) is more than offset by the smoother downstream 
water surface conditions and the reduction in splash and surging. A slight 
excess of tail water also moves the toe of the jump against the sloping face 0: 
the spillway (if one is present), stabilizing the jump and producing smoother 
overall operation. 


On the Center Hill, Conchos, Norfork and Fort Gibson Dams, Mr. Berryhi 
describes damage caused by bed meterial or other debris circulating on the 
apron. The writer has also witnessed damage of this type on other structure: 
and feels that abrasion damage is more wide spread than generally believed. 
Every effort should, therefore, be made, particularly by those engaged in 
hydraulic model testing to eliminate “dead” areas in basins where debris ma 


‘ become trapped. Insofar as possible, basins or buckets should be made self- 


cleaning and the bed material just downstream from the end sill should be 
stabilized or reduced in elevation to prevent bed material from being swept 
into the basin during periods of necessary unsymmetrical operation. It has 
been common practice on model tested basins to use an end sill which 
creates a ground roller and moves bed material upstream against the end 
sill to prevent undermining of the apron. The writer believes that strong 
ground rollers should be avoided by providing end sills with flatter slopes 
or dentated or notched sills which tend to neutralize the upstream motion. 
With continued motion of the bed material against the end sill, some bed 
material is bound to find its way into the basin and become the source of 
abrasion damage. 

In discussing the damage at Lucky Peak Dam, Mr. Berryhill indicates thai 
perhaps the damage was similar to that found in laboratory tests on concrete 
specimens subjected to the impingement of high velocity jets. In the case of 
the laboratory tests, the damaged area was in the form of an annular ring 
adjacent to the impingement area. The writer believes that at Lucky Peak 
Dam the damage resulted from an entirely different cause since the damage 
was located some distance from the impinging jet. It is very possible that a 
vortex formed at the upstream corners of each slide gate, similar to those 
which form behind radial gates on a spillway. At small gate openings, the F 
low pressure core of the vortex was directed toward the floor, as mentioned 
by Mr. Berryhill, and deflected inwardly by the gate slots. Cavitation then 
occurred in the vortex core near the floor. Downstream, a high pressure _ 
area collapsed the cavities and damage resulted as shown in the photographs 
The extreme damage further downstream was probably a progressive result 
of the damage which began closer to the gate. The eroded areas in Figure i 
indicate that the low pressure area originated behind the gate. The vortex Y 
theory should be investigated at the first opportunity in a hydraulic model. 
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) DIVIDED FLOW THROUGH A DIVERGENT INLET CONDUIT* 
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SYNOPSIS 
| This investigation is concerned with the determination of the characteris- 
S of divided flow phenomena. Two models have been analytically and ex- 
»rimentally investigated. One is a conduit with a constant cross-section 
aving an opening perpendicular to its length. The other is a conduit witha 
-Vergent inlet having a side opening at the end of the divergent. The analy- 
‘sis developed by means of the free stream line theory, assuming two- 
-mensional irrotational flow of an incompressible fluid. The experiments 
ere conducted under approximately two-dimensional flow conditions. The 
.eory, together with the experiments, has shown quantitatively the influence 
the velocity distribution in the inlet on the characteristics of the divided 
Ow phenomena. 


INTRODUCTION 


Efficient diversion of a portion of a main flow for auxiliary purposes is 
tally important in designing large scale hydraulic systems. Experiments 
ive shown that a unique angle exists for a branch conduit which will keep 
sses ata minimum. A conduit with a slit or hole in its side produces a 
ee jet stream at a certain natural angle to the main flow. This angle cor- 
‘sponds to the angle of least loss for the case of a branch conduit. 

The objective of this investigation was to determine the parameters affect- 
g this natural angle and other flow characteristics and to study their influ- 
ice on the jet diverted through the slit opening. 

Earlier experimental work(1,2) and analytic studies(3,4,5) have shown that 
e natural angle and other flow characteristics such as the contraction coef- 
sient, width of jet etc. depend on two factors. These are (a) the geometry of 
¢ junction and (b) the ratio (r) of the discharge diverted into the branch to 


te: Discussion open until May 1, 1958. Paper 1492 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
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the total discharge in the main conduit beyond the junction point. Previous 
scale model tests conducted by the Bureau of Reclamation(6) have shown 
large discrepancies with the results obtained in small laboratory tests and 
also have differed with prototype measurements. 

At the suggestion of Professor A. Craya, the present study was initiated to 
determine the influence of the velocity distribution in the upstream sections 
of the main conduit on the natural angle and other characteristics of the flow 
diverted through the slit. No known experimental or analytical work had 
considered this effect in the past. 

The problem of a general velocity distribution in a main conduit preceding 
the junction point is neither theoretically tractable nor experimentally feasi- 
ble. The mathematical difficulties of the problem are increased greatly be- 
cause the use of the complex variable technique is not possible in a rotational 
flow. In addition the problem of altering the velocity distribution at the inlet 
of the conduit by laboratory means is difficult.(7) However, the use of a 
divergent inlet conduit offers a means of altering the velocity distribution. 
The divergent inlet conduit is advantageous in that a small angle of divergence 
produces considerable changes in the velocity distribution, as shown by ‘ 
Nikuradge.(8) Ordinarily such velocity distribution changes would require an 
inlet length of 20 to 50 feet before the actual experimental section. The in- 
fluence of the angle of divergence itself, i.e., in the case of irrotational flow 
on the jet characteristics can be determined analytically. This influence has 
been shown to be small for small angles of divergence. 

The theoretical part of this investigation dealt with two models. One is a 
conduit with a constant cross-section and a side opening, while the other is a 
conduit with a divergent inlet having a side opening at the end of the sloping © 
part. Both cases were investigated as problems in two-dimensional irrota- 
tional flow of an ideal incompressible fluid by means of the free streamline © 
theory. The experimental part of this investigation dealt with the same two — 
cases of straight and divergent inlet conduits under approximate two- 
dimensional flow conditions. 

Comparison of theoretical and experimental results obtained for the 
straight conduit model served to determine the degree of the validity of the 
free streamline theory together with the previously stated assumptions. The 
second model with the divergent inlet conduit was used to determine the in- 
fluence of the non-uniform velocity distributions. 

If the deviation between the theory and the experiments in the case of the : 
divergent channel is greater than the differences observed between the theory 
and the experiments in the case of the straight conduit, this deviation may be 
attributed to the influence of the velocity distribution. This was assumed 3 
since all other parameters (ratio of discharges, geometry of the opening and 
property of the fluid) were kept constant and only the velocity distribution 
was changed. 


In summary, the present investigation was concerned with: 


1. the theoretical determination of the natural angle of the jet and other flow 
characteristics for divided flow through a conduit with a constant cross- 
section as well as for a divergent inlet conduit, based on the free streamline 
theory of irrotational flow; 4 
2. the experimental determination of the natural angle of the jet in both 
models. The angles used were: ¢ = 00, $= 2045'48", ¢ = 6925! 43": 

3. the magnitude of the influence of the ‘velbcity distribution at the inet on 


ASCE TSAKONAS 1492-3 


the characteristics of the jet diverted through the side opening. This also 
required measurement of velocity distributions in the experiments. 

The experimental results showed agreement with the free streamline 
theory when the velocity distribution at the inlet was uniform as in the case 
pf the straight conduit. However, when the velocity distribution was non- 
uniform as in the divergent conduit, the experiments showed a deviation from 
the theory, but they remained similar in their general aspects. 
| The experiments thus showed that the velocity distribution in a inlet con- 
Huit has a definite influence on the flow pattern of a divided flow field. 


: 
) 
) Analysis 

The derivation of a theoretical solution for the flow characteristics of the 
~roblem of a divided flow through a divergent inlet conduit will be considered 
in the following analysis. 

The inflow through the divergent channel A", A', N and P was divided into 
flow through a slit (PQ) and a straight conduit (NBQ) as shown in Figure 1. 
y successive transformations the physical plane was transformed to the 

nodograph plane (Figure 2), then to the logarithmic plane (Figure 3), and 
inally, by the Schwartz - Christoffel transformation, to the £ -plane (Figure 
) where physical properties of the flow pattern could be identified easily. 

By considering the physical plane with a source of strength (q) at the point 
of intersection of the sides (AP and AN) of the divergent channel, the result- 
ing idealized flow is such that the field of interest is bounded by streamlines 
along which the velocity is either constant in direction or constant in magni- 

de. The line AO (O is a stagnation point) is the boundary which divides the 
Flow field into an upper part discharged through the side opening (PQ) and a 
lower part conveyed through the main conduit (NBQ). The point N is another 
tagnation point of a different nature since it does not split the flow field. 

__ The hodograph plane represents the complex conjugate velocities of the 
‘low field, hence 


| 


W =u - iv, (1-1) 


where u and v are the velocity components of the point (x,y) in the z-plane. 
The logrithmic plane is related to the W-plane by 


U = log W. (1-2) 


[he origin is chosen at the point P. This requires the magnitude of the 
velocity at P to be equal to unity as well as to the points Q and C. This in- 
volves no loss of generality. 

The ¢ -plane is the result of the Schwartz-Christoffel transformation: 


= = e 24 (1-3) 
Ba ot Me cew) es p(cev) Moar (ten) 
yhere P is a complex constant and the points Q, P and N have been chosen at 
*=0,+, -1, respectively. 
_ Equation (1-3) may be written in a partial fraction form as 


A B C 
| opted) © 


j 


| 
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where 


(1-5) 


Using t= /€ and spliting the terms of Equation (1-4) into complex conjugate 
forms and then integrating them yields: 


Vt -ivo. -iJ/a B eae e Sead meas ( Vo -ivw a9 
U = —_—1log |———_——"— —log |—*~ ‘Log _ 
iva Vt +i va Per Vo taVvy iJVw Vo tivw 


The condition Up = 0 at t = + ~ has been used to eliminate the constant of 
integration. 

The constants of the transformations have been determined by relating the 
€ -plane to the U-plane, since on traversing the boundary in the U-plane the 
change of arguments of U is known. 

Using Equation (1-6) in the following form 


re) gw) el IE Sl 


in order to describe the singular points O, A and N leads to 


petoun --f.-x, alas (1-7) 
i vo i Va ivy 


However, it must be kept in mind that the constants A, B and C of the 
transformation are purely imaginary since P can be readily shown to be 
imaginary.(10) 

Combination of Equations (1-5) and (1-7) yields 


1 apie ( = (Hh - 1)(@-w) 

Ve A Ca=n)Ca-B) =a) =? (1-8 
WwW ey 

wg = {be wit wht a= : cape 


where rv = 1 is due to the transformation. 
The substitution of /a@ = aand /W = w in Equation (1-8) leads to: 
peak tata We alay har Eas (1-9 
H(a +a+1) -a° - pw - aw | 


and 


t 
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| k(w- w')( w+ a)(a- 1) =wp-D(p-a®) (1-10) 


}quation (1-10) is quadratic in w and yields the solution 


|é eii(a=1)(aew~)+(1+a~ wl = V [k(a-1)(a-w)+/( L+a” ae] “+h [eC an?) <wikc(e-1)aw+a yl 
2 {k(a-1)-w] 

(1-11) 
there the negative sign is the proper one satisfying the relation “>a, im- 
josed by the transformation. (See Figure 4) 
| The complex potential function (Z) in the ¢ -plane for a source of strength 
,at ¢ = -a and sinks of strength q' and q" at ¢ = -Band ¢ =¥Y respectively 
6 given by 


Z=- Frog t+a) +L log(t+ p) +2 rog(t-7) (1-12) 


| The complex conjugate velocity and complex potential are related by 


ES (1-13) 


} 
etting W = Oat the stagnation point (f = -w) leads to 


a-p_..Y+8 (1-14) 
) Q-W 7 ¥+W , 


here ve ue with r being the ratio of the discharge diverted into the branch 


» the discharge into the main conduit before the junction point. 

_ Substitution of Equation (1-7) for the values of A, B and C in Equation (1-6) 
nd use of Equation (1-2) gives the velocity at any point in the field by the fol- 
owing: 


| 
| 
| 
| 
s a result of this Equation (1-13) leads to 


(1-16) 


wo OP feces ede ak bag Bi 7 
sheet desde ca: Raoreas al Corie 
ere 

k k 


Mena) can) rae) 
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Equation (1-16), which relates the ¢ -plane to the physical plane through a 
differential form, becomes 


- .@9 (a-p)(a+y) C+w £ 
Wed lt Rome par EFanteptayy KPa Se 


after using Equation (1-14). At this stage of the development the parameters 
a,W, and ¥ are still undetermined. An additional relation between these 
parameters is evident from the geometry of the configuration shown in 
Figure 1. 


h! 


Be tad 2 


which, by means of Equation (1-17), yields 


(Q-8)(Q+Y) tand f M(t+w)o(¢ ) Caen? (1-18) 
a=W rede LrGNePUt-¥) 2a 
Combining the expression of the velocity at B, found from Equation (1-15) for 
¢ = -B, with 
wl a ' h! ca -r h" 
1 os aah dk Slate 
leads to 
k k 
Bi WO eve a ) VB -Vw (1-19) 
PRRs VE va 


This relates the width of the diverted jet to the height of the conduit. 
Use of Equation (1-19) together with (1-14), changes (1-18) into the form: 


- fe C - 
+Q@ + = 
MSE GSE Tes? (1-20 
k k : 
m(p -w) a (eee al a . 
ern ton ele me Va+/F 1 -VF/\VF -Va/ - 5 


This integral relation will determine the parameter Y for givena +B 5 Wis 
Using the intermediate relation 


the width of the slit is given by Equation (1-21) after the integration of 
Equation (1-17) from £ = 0 to f =+0. 


7 
vi 
hay ey Ee eee 


ee ee ee) ee 


— 
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if mip +y)(a~ B) eae i ~JE\ (af ( Ss ) at (1-21) 
Pe ™ (+0) \ va-f6/ \l +e (orale + BEY) 


| The position of the stagnation point of the flow pattern is given by: 
sg a . o (O-8)(a+Y) [Gaps (0+w)a(o jas (1-22) 
n a-wW (t+aCc* By(t-y) . 


The contraction coefficient “c” of diverted flow is defined as 


Roe fa (1-23) 


: 
: 
| 
| 
: 
: 


vhere ¢ and h" are given by Equation (1-21) and Equation (1-19), respectively. 
| The angle of the diverted jet at ~ is given by setting ¢ = Y in Equation 
1-15) 


| Mo = (Sey (vey aOR ee 
° Crrve) ra) ASE) 


. 2 a -i0 = 
Since We = val e and wat =, 


hen 


=i, 
“18 | idk tan” =V@_- ik fan eee ta) tanh YO 
VY vr y 


Taking the principal values of the multivalued functions, the angle of the 
mpinging jet is given by 


pial ira token Yoyo tan ; (1-24) 
Aaya 


The analysis developed thus far permits determination of all characteris- 
ics of the flow pattern providing the ratio of discharge and the geometry of 
1e configuration are known. 

For simplicity the following method of computation was adopted in this 
iestigation. 

Assuming values of @ > 1 and 0<W <1, uw was determined by Equation 
[-11), 8 by Equation (1-9) and ¥ by the integral relation of Equation (1-20). 
he flow characteristics and the corresponding boundaries were obtained as 
slows: 

6 by Equation (1-24), (/h by Equation (1-21) , C by Equation (1-23), 

"/h' by Equation (1-19), S by Equation (1-22) and r by Equation (1-14). 

his computation procedure has the advantage of eliminating all but one of the 
itegral relations between the parameters. The integral relation which must 
2 used is given by Equation (1-20). 


{ 
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2 - Solution for a Straight Conduit 


Using the expression evolved for the general case of the divergent inlet 
channel, the flow characteristics for the straight conduit (k = ies 0) will be 
developed. 

For the case of k = 0 Equation (1-11) gives 

L= Lor t= 6; 


The former case requires “ =a = 8 =w = 1 andis a special case which will 
be discussed later. The latter case (“4 =a) requires 8 = 1, thus allowing 
Equation (1-15) to be reduced to 


Uee0 on 


= OO 


Jt +id/w 
When ¢ = -B= - 1 and ¢ = - a, Equation (2-1) becomes, respectively, 


1=-VW 


Wa = | "elas mee 


and 


v, | ~ Va -/w 
Va + /@ 


Using the latter expressions, h and h’ can be determined as 


Hiedee trees y ee WS 
v q 
A Va =- JW 


and 


ht = ee wioldeniee Lvs 


(2-2 
Oar: | 


oe yer 


Since h = h', Equation (2-2) takes the form 


5 
(VE+va)® (1+vo)* (2-3 
at+y Lit fi 

} 

after using Equation (1-14). Once the values for « and w are chosen, the - 
value of Y may be obtained from Equation (2-3). : 
The angle of diverted jet is determined by Equation (1-24) as 7 


4 1 Vw 4 
@= 2 tan 7 3 ( 
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The corresponding ratio of discharges is given by Equation (1-14) as 


sama ene 


is el = Y+W (2-5) 
Soin e 


The width of the slit, upon integration of Equation (1-14)(10) leads to: 


efa' = > [(a+w)(  *1)loga-( ¥-w)( a-1)10g +2 VOV7(a=1)s] 
T(1+J/@ )“(a+ Y) 

) (2-6) 

| The corresponding contraction coefficient given by 

j 

; 

) ee eee eee el) A) (2-7) 


(@+W)( ¥+1)loga-( Y-w)(a -L)log ¥+2 JW JY (a -1)n 
and the ratio of the width of the jet to the height of conduit is determined by 
Equation (1-19) as 


| Hl eee) yer) 
. UF RO = See eee 2-8 
: BS (erY Gari)” ae 


Two limited cases (r = 1 anda >~ with 0< r< 1) were investigated and 
are described as follows: 


1) Case r = 1, 


2 
yi = ba =( Vo- 1) (2-9) 


then WJ = 1 and 5 
C/O) rely . 


Ifo > 1, 
then ee aS 1 
and 
9 = 90° ae -& +0 and C = 5-B— = 0,611 (2-10) 
a—_ 


This limiting case corresponds to a flow through an orifice on an infinitely 


large container. 
Iia>- 


° Ae ie re 
then Lim | = co” = 3; 6100 verti ht ee 


independently of the values of ¥ andw, 


¢ 
and lim Ci=0; 
a> 


See eo . 
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This case (r = 1,®> ~) corresponds to the flow pattern through an orifice, 
the dimensions of which may be considered infinitely large in respect to the 
height of the conduit. The jet is diverted through an angle of 60° without any 
contraction. 


2) Casea>o and O0<r<il, 


then parameters Y and w are determined, respectively, by 


2 
——e-=l-r and G+Jw )" = 1 
1+ Jw LeSy: 
and the characteristics of the flow field are determined by Equations (2-4), 
(2-6), (2-7) and (2-8). Noticing that w decreases to zero faster than ?, it fol- 
lows from the characteristic relations of the flow pattern that 6 = 0, 
ne = 0 and C = 0 independently of the values (/ h', when w= 0. 

The case &@ > © corresponds to an orifice with an infinitely large opening 
compared to the height of the conduit. The jet is diverted without any contrac- 
tion and the quantity h''/h' takes a maximum value for each ratio of discharge. 
The jet emerges through the free surface in the vicinity of the downstream 
edge of the slit which extends an infinite distance from the upstream leading 
edge. 

Results for a straight conduit solution as well as for its limiting cases, 
are plotted in Figures 5,6 and 7. Although it does not seem possible to show . 
the identity of the expressions developed in the present analysis with the cor- 
responding expressions developed by Craya(3) and McNown,\®) the plotted re- 
sults show complete agreement with Craya’s results and good agreement with 
McNown’s results. The small discrepancy in the latter case is probably due 
to the fact that, in reducing the results to the same basis, use was made of 
values obtained from the graphs plotted by McNown. 


3 - Solutions for a Small Angle of Divergence 


After the successive transformations of ¢ = - o and /O =t the integral 
relation (1-20), which determines the parameter ?, leads to 


p [OMe VP _ (tavll\( te) gg af sa) 
Ja (t ~ a) (t°-p) (t7+ y) t4Vq_ /\to tan0(@- B)(p +y 
k k 
aa = VR 
a+ JP 1 JB 


Upon integration and substitution for B by Equation (1-9) this takes the 
form 


(3-1) 


£027 50s K)o= OY, Se) (3-2) 


Having these two functions f(a, Y , w , k) and #(@, ¥,w, k) which are 
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alytic functions of k in some common region (R) with respect to all 
arameters (@, ¥,w) and the independent variable k, we actually postulate 
that only on some surface S(@, ¥ , w), which belongs to the regions of analy- 
jticity, does the equality (3-2) hold. The functions may be expressed as power 
gseries in k in the region of common analyticity: 


| 
Ca, YW, HAL a,Y,W) HALA, FWA ay WK, 


la 5 ¥,W,k)=B(a,y,w)t+Bla,y,w)k+Bla yy pw )kt.. 


[Then on the surface (S), the equality (3-2) implies that coefficients of like 
faewers of k are equal 


) AC a,Y,W)=B(a,y,w ) (a) 
| A(a,7,w)-B(a,7,w) (b) (3-3) 
) A(a,Y,;w)=B(a,7,w ) (c) 


. 
Equations a, b, c, etc. define several surfaces ina, Y,w space. If any two of 
e surfaces differ they may or may not intersect. If they do intersect, Equa- 
ion (3-2) would apply only to restricted values of one of the parameters, 
mamely those satisfying the equation of the intersection. However, this con- 
itradicts the postulate based on physical grounds that no restriction exists on 
any of the parameters in the region P of the analyticity. 
If the surfaces do not intersect, no values of the parameter exist which are 
common to both surfaces. Hence, there is no surface for which (3-2) holds. 
This also contradicts the postulate that S(a, 7, W) must exist. 
Consequently, the remaining equations in (3-3) must be identically satis- 
fied or equivalent to (3-3a) which is the only one that specifies Sia, 7, w). 
Equation (3-1), expressed in k powers with attention being concentrated on 
only the first term,(10) leads to: 


1 (va+Vq@)* 


+ B,( A, Py %yW )FkKE( A, Py %yW )tee0™ 
(a= p)(a+7) 


(3-4) 


(Vp + va)? 0(25 B,Y,W) + kO(a, By Y,W )*.. 
(Q- B)(p +7) 


Fie 


On the other hand were Equation (1-9) written as 


PE CES Ness wien w) eg woe (3-5) 


p(a‘ta+1) - a -pw -aw 


Equation (3-4) then obtains its final form: 

: (Jie Je) = +F(0, 7,0 Jour, (A, YW )+...8 
a 
k 


yO O(a, Y,W )+kd,( 2, te eee 


(3-6) 
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Therefore by the argument stated above, the surface S(a, ¥,w) to which 
Equation (3-1) holds is given by 


(Va+J/a)® | (14+Vu)® (3-7) 
OSes y ike A ‘i 


This is the same relationship as Equation (2-3) found in the analysis for the 
case of k = 0, without using the integral relation (1-20). 

It can be readily found that the integral relation should include the case of 
k = 0 by means of Equation (3-5), since 


k 
lim At BS) — + 21 
k--0 1- /B 


p-1 


and the right hand side of Equation (3-1) does not vanish. 

The discussion in the preceding paragraphs eliminates the need for evalu- 
ating the integral in Equation (3-1), which was attempted both by a closed- 
form method and by numerical computations. The closed-form method was 
not tractable, while numerical computations were arduous, difficult and did 
not maintain the required accuracy. 

By means of Equations (1-9), (1-11) and 3-7) the values of 8, u and ¥ can 
be determined for assumed values of @> 1 and 0< W <1. The flow charac- 
teristics then can be computed in close-form except for the ratio of the slit 
opening to the height of conduit ((/h' ) and the contraction coefficient (c), 

The investigation of some limiting cases leads to: 


1) when u> ~, 
Equations (1-11), (1-9) and 3-7) yield, respectively, 


2 
a- (+2) : 


a a+ yew (8-8) 


ae (b) 


and 
y- Lva+ Va) Vo] Vo 
NN (0) 
Va+1 +2VW 


For the case H >~ and assumed values of 0<w <1, the values ofa, B, 
and Y can be computed by Equation (3-8 a, b, c) afterwhich the flow charac- 
teristics 6, r, h'/h' can be determined by the general expressions given in 
section 1. This case physically corresponds to a conduit with an extremely — 
large opening in its side. The width of the diverted jet is smaller than that 
attained in the case of straight conduit. Furthermore, it is noticeable that, 


as the size of the slit opening increases and tends to infinity, the width of the 
jet diverted through the opening approaches zero. 


2) When U >= and a= 1, 
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| 
| Equations (3-8) leads to: w >oB> 0.3333 and Y >0 (¥ does not approach 
Pro as quickly as w) 

jor this case @ = 0, h"/h' = O andr =0. The entire discharge is conveyed 
hrough the main conduit. 


) Whena > 1, 

Equation (1-10) leads to & = a = 1, independent of the value of w. Then 
jquation (1-9) gives B the value of 1 for any value of w, while Equation (3-7) 
w (3-8c) yields Y= VW . However, Equation (1-14) has meaning only if 
»= 1, which leads to: 


: 
| 
| epted, {O--'s and h"/h' — 0 


this case corresponds to the physical case of an infinitely large container 
jith a small side opening. The jet is diverted through an angle of 90° and 
me angle of divergence has no influence on the jet’s characteristics. 

Up to this stage of the analysis no approximation has been made. The flow 
naracteristics r, 9, h''/h' have been determined in closed-form eral 


nd the results are plotted in Figure 7 for the cases of oe 64 and @= ae 


| In order to reach a conclusion on the influence of the angle of divergence 


4 the case of an irrotational solution, the study of the physical characteris- 
ics (/h' must be considered. At this stage it is ees that k is small, 
e., the angle of divergence ($) is smaller than a = 6926'. Equation (1-21) 


be written in an expanded form, where only the two first terms are taken 
to account(10) as 


Lint « of™ £( jab + i2ck ie £(6 OC f ao 4... , or (3-9) 
Apa = cI 


1 + i2ckI, 


| CD 
| iO) s coaGEe CEE. 
' o( Ss) = tan Te - sani oe 


k 
eye eral B) ae cea 
n( JB +Va)* Va - JE + VB 


The first term (I,) of Equation (3-9) is the value of (/n' for k = 0 and its 
fficient (c) reduces to 


al 
1 on( +Y )(1+/a)* : 
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The integrand of the second term Ig is a combination of the integrand of Ij 
and the function $(¢) which is bounded. Its maximum value is attained as 
é=VarewAsia result 


max.0(¢) = fan ss 2 tan \ Va 


Therefore, the following inequality holds: 


|2cx £(6)0(E Jab | < 2c,k0(VE) fey We £($)at |~ 20,k0( VE) | 1, | =2", 


where 


k k 
Naa) ee) 


The actual value of the second term of Equation (3-9) is much smaller than 
I'5. It can be shown by calculation that I'y = (1.5 = 3%) Ij. Since I, <1, the 
use of the second term (Ig) in Equation (3- %) with the first term (Ij) a alter 
the values of (/h' to the second or third decimal point. 

Therefore, for the case of a small angle of divergence and for a < 9 which 
covers the range of practical applications, it can be assumed that the lines of 
(/h' remain almost invariable as k varies from 0 to 45 35° 


Using this fact the theoretical results obtained from 


h’ 
(Figure 7) can be transferred to @ versus r plane (Figures 10 and 11) which 
has been used for comparison with experimental work. 


versus @plane 


4 - Experimental Investigation 


The main features of the equipment and instrumentations as well as the 
experimental procedure are discussed briefly in this section. The entire ex- 
perimental work was conducted in the Fluid Laboratory at Columbia Univer- 
sity. The fluid used in the experimental work was water. The open channel 
in which the experiments were conducted was 21 feet long, 15 inches wide an 
37 inches deep. It is part of a circuit formed by a reservoir, pump, constant 
head tank, entrance section, two exit tanks with measuring weirs and a re- — 
turning channel leading to the reservoir or to a volumetric measuring tank. 

The actual test section was made by inserting two carriers (Figure 8) a 
a partially movable bottom into the channel so that it was divided into four 
parts: (1) an entrance chamber, (2) a main conduit, (3) an upper chamber 
and (4) the remainder of the channel: The lifting mechanism under the par- 
tially movable bottom provided a maximum angle of convergence and diver- 
gence of some Seven degrees in the inlet section of the main conduit. : 

The measuring equipment consisted of a pitot-tube with a traversing 
mechanism, standard point gages, a protractor, a micromanometer, a set of 
tubing to feed sufficient quantity of dye into the test section and a pair of tri- 
angular weirs to measure the two discharges. 


| 
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| The pitot-tube consisted of a streamlined stem and three removable tips 
jo permit reading of the differential pressures at three different stations of 
jhe inlet conduit. 
| The protractor consisted of a quadrant which was divided into one degree 
jivisions, a carriage which bore the pointer and a micrometer head used for 
small adjustments. The protractor was used for measuring the angle of the 
iverted jet. This was done by centering the pointer in the jet flow and by 
aking it parallel to the jet. 

The main problem in the experimental work was to regulate and control 

e discharge of water passing through the slit and the main conduit, to read 
ihe angle of the diverting jets and to determine the fluid’s velocity at different 
ections in the inlet conduit. 

Regulation of the discharge was achieved by calibrated weirs of known dis- 
‘-harge coefficients. Main readings were taken of the two discharges, the 
wngle of the side jet, and the velocity distribution in the upstream inlet sec- 
tion. In order to read the angle of the jet, a sufficient quantity of dye was in- 
ii to make the jet visible. A picture also was taken of each test. An- 


: 


ther important feature of this investigation was the velocity distributions 
used for comparing the results of straight and divergent inlet conduits. Dif- 
serential pressure between the pitot-tube and hydrostatic wall tap was mea- 
sured by means of the pitot-tube and micromanometer. The velocity distri- 
Fgsons was obtained by 


| ve = 10.0 A p ft.“/sece® 

| p 

vhere A p = differential Sig A in inches of water 
ind p = density in slug / ft.3 

| Velocity distributions were ae at three different stations, 2 7/16, 6 3/8 
ind 12 5/16 inches from the leading edge of the opening. The results were 
{ptted in a non-dimensional form y/d versus u/um (see as an example 

igure 9) where 

y = distance from the central line of the conduit 

d = 1/2 of the corresponding height of the conduit 

u = velocity at the point under consideration 


Uj, = Maximum velocity at the section under consideration 


| Table I gives values of hj, dj, 1;, and the corresponding area (A;) for the 
‘arious cases under consideration. 


: 
: 
: 
: 


| Experimental work was carried out for a straight se dees as well as for a 
| a inlet channel with angles of divergence p= a ithe 2°45' 48" and 


oo = 6925'40". Several slit openings and various discharge ratios were 


4 _ each geometry. 

In the straight conduit case the velocity distribution was almost uniform 
<cept at the boundary later. When an angle of divergence was used, the de- 
iation of the velocity distribution from uniform became more and more pro- 

ounced as the angle was increased. In most cases the velocity distribution 
‘as nearly symmetric, the exceptions being cases in which separation was 


1minent. ; 
| & plot of the experimental results on charts representing the theoretical 


; 


Ss 
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Table I 


Station 3 


Station 1 
£, = 12 5/16" 


(,= 2 /é" 


Angle of Scheme 


divergence 


tation 2 
e = 6 3/8" 


h, =1.905" h, = 1.905" 


1 
d.= 0.9525" 


0.198 ft." 


d, = 0.9525" d, = 0.9525" 


O.198rt.? | A, = 198hft.* 


1.9852" 


1.5918" 1.3001" 


b = 2°58" 0.8926" |4, = 0.796" 


Q 
i] 


0. 6506" 


041859ft.- 0.1658ft.* | A, = 0.135Ufte* 


1.6303" 1.1867" h, = 0.5177 
b = 6°25) 3" O.8152" | d, = 065931" a, = 0.2588" 
0.1698ft.* 0.1236fte> | A. = 0.0539fte.’ 


solution shows a close agreement in the case of a straight conduit with uni- 
form velocity distribution (Figure 5), but depicts a definite deviation in the 
case of a divergent conduit with non-uniform velocity distribution. As the 
velocity distribution becomes more and more non-uniform the deviation of 


the theoretical and experimental points becomes more and more pronounced 
(Figures 10 and 11). 


CONCLUSIONS 


The free-streamline theory is the only one which offers an analytic solu- 
tion to the divided flow phenomena, but the stumbling-block for the generali- 
zation of its principles was believed to be attributed to the velocity distribu- 
tion in the inlet. The main goal of this investigation was to determine 
qualitatively the influence of the velocity distribution on the solution obtaine: 
by the free-streamline theory. 

Both the experimental work and the analysis established the fact that the 
upstream velocity distribution is another essential factor for the determina 
tion of the divided flow pattern by the free streamline theory. ~ 

This conclusion has been established by comparison of the theoretical 
solution, obtained by the free streamline, with corresponding experiments fi 
several different angles of divergence, which are equivalent to considering 
several different velocity distributions. ‘ | 

As long as the upstream distribution remains essentially uniform, the fr 
streamline solution is verified fully by the experiments. However, when the 
velocity distribution is not uniform, the characteristics of the diverted jet 
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phows a systematic difference, but remain similar in general aspects. 
| An angle of divergence of ¢= 2945' 48", which alters the velocity distribu- 
tion slightly from uniformity, lowers the angle of diverting jet by 1 1/2 to 4 
Hegrees. An angle of divergence of @= 6925' 43"', which alters the velocity 
fecross., considerably, lowers the angle of the jet inclination by 2 1/2 to 
degrees. 
| The similarity in results found between analysis and experiments suggests 
An attempt to incorporate the velocity factor to the principle of the free 
| ei theory. This may lead to a generalization of the free streamline 
rheory for rotational flow, which at the present is relatively intractable from 
theoretical stand-point. 
The analysis used in this investigation leads to a close-form solution for 
the flow characteristics in a straight conduit case as well as in the case of a 


divergent inlet conduit with a small angle of divergence (b<35)- However, 
for any angle of divergence <0, close-form expressions are obtained for 


the natural angle, the ratio of discharges and the ratio of the width of the jet 
to the height of the conduit. 

Figures 5 through 10 show the variation of the flow characteristics for 
different slit openings and discharge ratios. These variations are discussed 
in Reference 11 as a function of the pressure variation. 

The analysis developed herein, based on the free streamline theory and 
applied to a small angle of divergence of the inlet conduit, is expected to be 
useful for design purposes or/and as a guide for the interpretation of experi- 
mental data for similar geometrical configurations. 
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